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ABSTRACT 


This investigation is concerned with the optimization of 
arch structures. The DOT optimization code is used to minimize 
the volume of arch structures which are constrained by limits 
on stress, design geometry, and section dimensions. Modeling 
the arch structure by a series of bar-beam elements, the 
finite element method is used to compute element stresses. The 
DOT optimization code selects section dimensions to prevent 
failure due to element stresses exceeding the material yield 
Stress. Specifically, through coordinate transformations 
between oC asl element coordinates and global system 
coordinates the element stiffness matrices transform into the 
global stiffness matrix. The resulting system matrix equations 
are then solved for the system degrees of freedom, that is, 
displacements and slopes. The system degrees of freedom, in 
turn, are transformed back to the element level to compute the 
internal forces and moments and hence, the stresses. Results 
are presented for a number of caseS with regard to 


optimization scheme and stress analysis. 
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I. INTRODUCTION 


A. BACKGROUND 

Over 5,00 years ago, evolution of Ene post and lintel 
structures of the stone age gave rise to the arch. Highly 
regarded for its graceful shape and design suitability, the 
Simple arch structure has been applied to engineering and 
architectural designs ever since. The ancient Roman Coliseum 
and agqueducts, great cathedrals of the Middle Ages, and 
railway bridges of modern history are just a few of the many 
examples of structures comprised of arches standing today 
eevomre 1.1). Throughout its history, engineers and architects 
have labored to improve the design of the arch in order to 
enhance the overall design structure. This desire for 
perfection has led engineers to devise a rational, directed 
design procedure and hence, the concept of optimization was 
created. 

The advent of the computer era has lead to 20 years of 
extensive development in the use of numerical optimization 
techniques. These techniques offer a logical approach to 
design decisions where intuition and experience previously 
prevailed. Coupled with trends toward material and cost 
efficiency, numerical optimization has prompted considerable 


research in the field of automated design [Ref. 1]. As 


(a) 
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Figure 1.1 (a) The Coliseum 
(b) Firth of Forth Railway Bridge 
(c) Isernia, Italy, railway bridge 
Photos by Mrs. P. Menzies and CDR D.C. Warner 


a step in design optimization of structures, the arch has been 
the subject of numerous optimality studies to enhance 
applicability in engineering and architectural designs. 

One such study was performed by Farshad in 1976 
[Ref. 2]. Using calculus Ot variations, he derived 
optimality conditions for nonlinear partial differential 
equations for hinged-hinged arches. The total potential energy 
of the system, augmented with several objective functions via 
Lagrange multipliers, was minimized with respect to design and 
state variables to achieve equilibrium and optimality. The 
nonlinear systems of equations for optimal thrust, minimum 
length of the arch, and minimum volume were presented but not 
solved. 

In 1980, Rozvany et al. [Ref. 3] used the Prager- 
Shield criteria to optimize statically determinate arches. His 
‘arch’ consisted of two inclined funicular frame beams 
ridgedly interconnected with a concentrated load applied at 
the joints. In the optimal ‘arch’ only bending or axial forces 
develop depending on the ratio of 4L/D, where L is the span of 
the structure and D is the depth of the cross section. Ratios 
greater than eight to one produced axial forces only and the 
optimal shape has a height of half the span. Ratios smaller 
than eight to one develop only bending and the optimum 
Structure is a straight beam. In each case, the width of the 
beam segments for the optimal ‘arch’ varied linearly from the 


hinged support to the axis of symmetry. 


That same year, Lipson et al. (Ref. 4] used the 
‘complex’ method to optimize parabolic arches subject to 
uniform loading. His ‘arch’ was comprised of equal length 
Straight beam sections of thin walled rectangular tubes. 
Maintaining constant depth and width for each segment, the 
vertical and horizontal wall thicknesses determined the arch 
shape which was optimized for minimum total weight. An arch 
with a rise of 0.342 times the span length proved to be the 
optimum. 

Tne cies Ang et al. [Rec 5) solved the arch 
optimization problem by parametricing the unspecified arch 
axis using spline functions and employing a smoothing function 
to approximate the non-smooth objective function. The ‘arch’ 
was considered to be a ‘plastic’ design of rectangular cross 
section subject to bending and axial compression. Three types 
of boundary conditions were imposed, simply supported-simply 
supported, clamped-clamped, and simply supported-clamped. The 
optimum shape of the arch is claimed to be a parabola with a 
rise of 0.433 times the span length. Apparently, there is some 
disagreement between these results and those previously noted. 

In addition to arch optimization studies, Ding and Esping 
[Ref. 6] solved the minimum weight design problem for 
frame structures when stress and displacement constraints are 
considered. Using dual numerical methods, seven cross- 
sectional shapes were treated by approximating the stresses 


with pseudo and virtual load techniques. Results were 


presented for a beam clamped at both ends, a portal frame, a 
Pee oegriliage, and a helicopter tail boom structure. Although 
Ding and Esping’s investigation does not specifically solve 
for arch structures, the approximations used are completely 
detailed with convincing results. 

In December of 1990, Charles Scott McDavid of the Naval 
Postgraduate School presented his thesis, "Weight Optimum Arch 
eeemeaures,™” which optimized circular arches subject to 
various loadings and end conditions. Specifically, he 
optimized arches segmented into rectangular boxes that varied 
in width only. Through his research he concluded that a 
bar/beam element model is a viable technique for the 
approximation of arch structures, and that an arch structure 
that 1s more statically indeterminate is more efficient under 
identical loading. Additionally, he proposed possibilities for 
future researcn which includes varying both the height and 


Pieri dimension, the major thrust of this investigation. 


B. PROBLEM DEFINITION 

In order to provide an in depth study, each of the cited 
investigations began with a problem definition and specific 
assumptions about the type of arch to be considered. For this 
investigation, the arch is defined as a structure of constant 
Curvature (i.e., circular arches) which when supported at both 
ends and loaded laterally develops perpendicular reactions. 


This is intended to eliminate thick walled curved beams and 


Straight beams which develop virtually no perpendicular 
reactions when loaded laterally. Additionally, the cross- 
section dimensions are small relative to the radius of 
Curvature and therefore the centroidal and neutral axes are 
assumed to coincide. Without the thin depth assumption, 
complications arise in the calculations of the displacements 
and the slopes because the arch no longer behaves as predicted 


by the beam equilibrium equation: 
Cl DV ee ee) (1.1) 


and the bar equilibrium equation: 
(AEu’)! =-P,(s) (1.2) 


where the prime superscript notation denotes differentiation 
with respect to the independent variable, s, and 


= Young’s Modulus 

= Cross-sectional Moment of Inertia 
= Lateral Displacement 

= Lateral Loading 

Cross-sectional Area 

= Axial Displacement 

= Axial Loading 

= the Independent Variables 


a 


* 
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In order to facilitate the development of a finite element 
code to approximate the local displacements, the arch is 
approximated by a series of straight segments. From the local 
displacements, the virtual load techniques, as described in 
the Ding and Esping paper, are applied to determine the 
internal psuedostresses. Once the stress distribution is 


determined, the arch volume is minimized to a structure that 


Maintains the developed stresses below the predefined maximum 
allowable stress. 

The thrust of this investigation is to minimize the total 
weight of a linearly elastic, isotropic, and homogeneous arch 
under a variety of loadings and end conditions. Optimization 
in this investigation refers to the variance of the cross- 
sectional dimensions (that is, the design variables) to obtain 
optimum least weight structures. Design Optimization Tool 
CDOT ) software (Ref. 7] is used ES perform the 
optimization subject to prescribed constraints on the design 
variables as well as on the stress limitations. The objective 
is to minimize the total volume of the arch while maintaining 
stresses below the yield strength of the arch material. The 
intent of this study is to provide direction and guidance on 
which further research for weight optimization may be 


developed. 


II. PROBLEM FORMULATION 


A. PROBLEM STATEMENT AND ASSUMPTIONS 

As noted in the introduction, the purpose SoOGmiea> 
investigation iS t0 optimize arch Structures to itera 
foundation upon which further research can be based. These 
arch structures, subject to Specified loadings and end 
conditions, vary in cross sectional geometry to minimize the 
weight. In order to limit the scope of this 'stua7 


approximations and specific assumptions are made as follows: 


¢ The arch maintains a conStant radius Cf Cucvatumaer 


* The arch iS approximated by a series of Straight segments 
of a solid rectangular cross sectional geometry. 


* Cross section design is restricted to ensure the 
applicability of beam and bar equilibrium equations (1.1 
and 182) 

° To prevent failure the internal stresses developed due to 
the loading must not exceed the yield strength of the 
material. 

° The arch structure is composed of a linearly elastic, 
isotropic, homogeneous material. 

To begin the design optimization process, the arch 
Structure iS approximated by contiguous Straight fame 
segments. Each segment is modeled by a bar-beam structure 


connects to the adjacent segment at a point defined as the 


nodal point. At each nodal point, the cross section base and 





height dimensions are selected as the design variables. From 
this model, the optimization problem can be formulated into 
objective and constraint functions which are functions of 


these design variables. 


B. MATHEMATICAL MODEL 

Due to the complex nature of this problem, the constant 
radius arch Structure is modeled by a series of straight 
contiguous elements where the arch radius of curvature, R, and 
the number of elements used to approximate the arch, NEL, is 
specified. (Figure 2.1) For simplicity, the length of each 


element is constant such that: 
1, = ORE Nee 


where @ represents the subtended arc of the arch. 





Figure 2.1 Arch Structural Model 


At each nodal point, there exists a base and height 
dimension such that the cross sectional dimensions from one 
element to the adjacent element maintains smooth piecewise 
continuity. (Figure 2.2) The resultant element shape is that 
of a three dimensional trapezoid whereby the volume is 
calculated by multiplying the average base and height with the 
length of the element. In mathematical terms, the volume of 


the ith element is calculated as follows: 


Volume (2) = Bae) Seen (2 ees (2.1) 
where — (B (i) + B (it1))/2 (2.2) 
(H(i) FRG) 72 (2.3) 


the Nodal Base Dimension 
the Nodal Height Dimension 
the Element Length 

the ith Element 
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Figure 2.2 Arch Elements 


10 


Defined in the problem statement, the optimal arch is 
achieved by varying the cross sectional dimensions, the base 
and height, in order to minimize the weight. Thus the nodal 
base and height dimensions are the design variables for which 


the objective function is defined. 


C. OPTIMIZATION PROBLEM 

The objective of this study is to minimize the weight of 
an arch structure while maintaining a stress distribution 
which does not exceed the yield strength of the material. 
Additionally, other constraints on the design variables are 
imposed. Since the arch is composed of a homogenous material, 
the weight of the arch is directly proportional to the volume 
of the arch. Thus, the objective of this investigation is 
satisfied by minimizing the total arch volume. The total arch 
volume, Vtot, is the sum of the elemental volume, v(i). Thus 
in mathematical form, the objective function is as follows: 


NEL 
Objective = MIN (V..,,) sap» via) (2.4) 


i=) 


where the elemental volumes, v(i), calculated by Equation 
(2.1), 1S summed for all elements to compute the total arch 
volume. 

In keeping with the assumptions made in the problem 
Statement, the objective function is constrained in order to 


impose practical and important physical restriction on the 


Pi 


problem. Properly defined, the constraints are used to avoid 
undesirable behavior such as yielding, to ensure validity of 
the governing equilibrium equations, and to provide a 
realistic design. For this study, the constraints falieaimeo 
three categories, strength criteria, geometric limitations, 
and Side constraints. 

First, for specified loadings and end conditions, the 
optimized arch must not ‘fail by yielding.’ Assuming the arch 
material to be linearly elastic, the applied loading must not 
cause the structure to exceed the elastic limit of the 
Selected material. Therefore, the internal stresses developed 
must remain below the yield strength of the material. 


Mathematically, the strength criteria is as follows: 
O (2) S757 
Or in normalized storm. 


(6 (i) (S77) 20 = 0.0 (2.5) 


where O(1i) is the maximum stress developed at the ith nodal 
point of the arch and Sy is the yield strength of the arch 
material selected by the designer. Unfortunately, the stress 
distribution, in terms of the design variables is not readily 
available. However, using the beam and bar equilibrium 
equations (1.1 and 1.2), a finite element scheme based on the 
model can be developed to determine the arch’s displacements 


and slopes due to a given loading. Knowing how the 


12 


displacements and slopes change throughout the arch, the 
stresses at the nodal points can be calculated. 

Secondly, limits must be imposed on the cross sectional 
geometry in order to ensure applicability of the bar and beam 
equilibrium equations. Limiting the cross section base and 
height dimensions relative to one another prevents the 
Structure from becoming either a shell-like or deep curved 
beam structure. To maintain the geometry of the arch, the 


following conditions are imposed: 

B(i) -3.0*H(i) $0.0 (2.6) 
and 

H(i) -10.0*B(i) $0.0 (2.7) 


Finally, the side constraints are imposed to ensure a 
Bealaistic solution. The arch is a physical object that must 
have a realistic finite cross sectional area; however, these 
section dimensions must also remain small relative to the 
radius of curvature by definition of the arch. Thus, the side 


constraints for the base and height dimensions are as follows: 


O-0sen. <= Bajisi6. an. (2.8) 


OS OBIE GSES Jel OE) = eRe Osi e)- (2.9) 


In the future, additional constraints should be considered 


Swen as Global buckling and local crippling. 
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III. OPTIMIZATION ANALYSIS 


To perform the computer optimization, the Design 
Optimization Tools (DOT) software package is used due to its 
availability, user friendliness, and reputation. DOT, a 
FORTRAN 77 optimization software package available from VMA 
Engineering, uses numerical search methods to seek a minimum 
value of one function, the objective, subject to the limits of 
others, the constraints [Ref. 7]. DOT has two methods for 
iteratively solving constrained optimization problems, the 
Modified Method of Feasible Directions and the Sequential 


Linear Programming Method. 


A. MODIFIED METHOD OF FEASIBLE DIRECTIONS 

Modified Method of Feasible Directions is a numerical 
method that deals directly with nonlinear problems. For this 
method, a search direction vector, S, is first found. The 
design point is then moved in this direction to update the 


design variable vector, X, according to the equation: 
xX =X ._+a°s (3.1) 
ox q 


where the scaler quantity a* defines the distance moved in the 


S direction, and q represents the iteration number. 
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For an initial design, say X,, the design is moved in the 
direction of the steepest descent until a constraint is 


encountered. 


wl Q =O 











) = — TRVE OPTINUM 
re) x? 


Figure 3.1 Modified Method of Feasible Directions: The 
Search Direction 


Having encountered the constraint boundary, a new search 
@erection is found by solving the subproblem: 
Maximize: 


MAXIMIZE: pry (3.2) 


les) 


Sub JeCES toe 


Ays0 (3.3) 
va ys 1 (3.4) 
where 
Sy 
S> 
yr: (3.5) 
>, 
LB 
0 
; (3.6) 
oS |: 
0 
1 | 
"g, (X) 
V’g, (x) 
oS : (3 37) 
Vg, (x) 
| V7F (x) | 


The search direction, S, will follow the constraint yet allow 
the design to leave a constraint boundary if the objective 
will reduce farther. In general, the form for inequality 


constraint problems is: 
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Maximize: 


-VF(x)° S (3.78) 

Stlesect "to: 
Vg ,(x) ‘sya JET (3.9) 
Se eS Al (3.10) 


When the search direction is away from a currently active 
constraint and the scaler product of the gradient of each 
Critical constraint with the S vector is less than zero, the 
constraint is omitted from the set of active constraints. If 
S is the null vector or numerically small, the optimization 
process is terminated because the Kuhn-Tucker conditions for 


optimality have been met. 


B. SEQUENTIAL LINEAR PROGRAMMING 

The second numerical method, Sequential Linear Programming 
(SLP), linearizes nonlinear objective and constraint functions 
and then obtains a solution using linear programming methods. 
Once the approximate solution is found, the functions are 
linearized about the new design point and the a linear 
programming problem approximated and solved. By repeatedly 
linearizing and solving the resulting problem, a precise 
solution is achieved. 

In general format, the nonlinear functions are linearized 


via a first-order Taylor series expansion as follows: 


ay) 


Minimize: 


F(x) = (gave x ) = OX (3.11) 
SUD JeCtmcon 
g,(X) © 9,(X,) +V9,(X,) > OXS0 j=l1,m (3.12) 
where 
OxX=X-X (3.13) 


and the zero subscript identifies the point about which this 
Taylor series expansion is performed. At the initial design, 
X,, the objective and constraints are linearized to give 
straight line representations of the functions. 

Typically, this method converges to the optimum solution 
with fewer iterations than the previous method mentioned. 
However, as seen in Figure 3.2, the optimum of the 
approximated linear problem is infeasible (1.e., a design that 
violates some or all of the constraints). Additionally, 
certain linearizations produce unbounded linear problems. 
However, imposing move limits on the linear approximation 


helps ensure that the optimum will eventually be reached. 


C. DOT PROGRAM PARAMETERS 

For both numerical methods, there are several parameters 
that can be adjusted within DOT in order to ‘fine tune’ the 
program for a specific problem. Fine tuning is a process in 


which the program parameters are internally adjusted to 
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a9 


optimize the optimizer performance. With proper tuning, the 
optimization process can be designed to remain within 
specified tolerances and operate more efficiently. A complete 
listing of all the DOT parameters is contained in Appendix A. 
However, for the purpose of this investigation, only the 
constraint boundaries, auto scaling, and termination tolerance 
parameters were tuned to enhance optimization performance. 
fRef. 7] 

First, for constrained optimization, the constreacm: 
boundary must be established. Mathematically defined, the 
constraint is considered active if its numerical value is 
between the value of CT and CTMIN, and violated if its 
numerical value is greater than CTMIN. By using a narrow band 
to approximate the constraint function, the optimizer is less 
likely to exceed convergence criteria without achieving an 
optimal design. In the realm of design, CTMIN is of particular 
concern. Principally, it is a small positive number 9am 
controls how far the design can deviate from the constraint 
boundaries and still be considered a feasible design. In 
theory, CTMIN can be reduce to zero to avoid any constraint 
violations, however, it is not practical due to the large 
number of iterations and computer expense required. 

In addition, it is normally considered good engineering 
practice to normalize design variables and nondimensionalize 
basic parameters [Ref. 1]. For optimization, variables are 


Scaled to affect normalizing by evaluating the diagonals of 
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the Hessian matrix of the objective and constraint functions. 
As the optimization proceeds, reevaluation is sometimes 
necessary to rescale the variables. The DOT parameter ISCAL 
may be selected to rescale the design variables over an 
interval or eliminate the scaling function all together. 
Unfortunately, the DOT manual indicates that there is no 
established theory for scaling. Scaling is therefore a 
meee ron Of trial and error. 

Last, the termination criteria also has a major effect on 
the efficiency and reliability of the optimization process. 
Termination criteria is established so that the design process 
is stopped when the number of iterations exceeds a specified 
limit. DOT parameters ITMAX and JTMAX specify the maximum 
number of iterations allowed for the Modified Method of 
Feasible Directions and the SLP method respectively. This 
ensures that the program will not iterate indefinitely. 
Furthermore, the progress of the optimization is checked for 
convergence. Design convergence is achieved when the change in 
the value of the objective function from one iteration to the 
next approaches zero. The DABOBJ parameter is a specified 
tolerance for which the maximum absolute change in the 
objective function between iteration is numerically small. 
Additionally, ITRMOR and ITRMST are parameters which specify 
the number of consecutive iterations for which the design 
change is less than DABOBJ for Modified Method of Feasible 


Directions and the SLP method respectively. 
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IV. STRESS ANALYSIS 


The objective of this investigation is to minimize the 
total weight (volume) of a load bearing arch subject to 
specified constants. To obtain an optimal structure, DOT is 
interfaced with an analysis program which computes the values 
of the objective and constraint functions in terms of the 
design variables, specifically the cross sectional dimensions. 
Since the strength constraint requires that the stresses at 
any point do not exceed the yield strength of the arch 
material, the stress distribution over the domain of the arch 
must be known. However, as indicated in Chapter II, the stress 
distribution is not readily available in terms of the cross 
sectional dimensions. Therefore the following stress 


development is pursued for optimization. 


A. STRESS DEVELOPMENT 

For this study, the strength constraint requires that the 
applied load will not cause the arch to fail by yielding. 
Therefore, the internal stresses developed must remain below 
the yield strength of the material. For this study, the 
stresses considered are composed of normal stresses due to 
bending moments and axial forces where the total normal stress 


is the algebraic sum of these components expressed as follows: 


LE: 


Oo, =6,+6, (4.1) 


foOtelememumal Stress 
Oo the normal stress due to bending 
0, the normal stress due to axial force 


where oO, 





Figure 4.1 Normal Stresses Due to Bending Moments and Axial 
Forces 


Shear stresses may also develop within the arch from shearing 


forces; however, the side constraints limit the geometry such 
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that these stresses are negligible. (See Appendix B for the 
complete justification for Shear stress omission.) 

To compute the two normal stress components, the arch is 
sectioned and approximated by straight frame elements. Thus, 
the stresses can be determined for each element endpoint (or 
nodal point) in order to establish the stress distribution. 
Each element is considered to behave as both a tapered beam, 
to calculate the stresses due to bending, and a tapered bar, 
to calculate stresses due to axial forces. 

First, for a straight beam segment, the maximum normal 
stress due to bending, hereafter referred to as bending 


stresses, is defined by the following equation: 


6, = (4.2) 


where c is the distance from the neutral axis to the point 
furthest from the neutral axis. The moment, M, at a section is 


calculated by: 
M=EIv"/ (4.3) 


resulting from the beam equilibrium equation (1.1)  aiaen 


substitution and simplification, Equation (4.2) becomes: 
6, = ECv!! (4.4) 


In the same manner, the normal stress due to axial 
behavior is determined. For a bar element, the normal stress 
due to axial forces, hereafter referred to as axial stresses, 


is defined by the equation: 
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F 
= 4. 
= (4.5) 


where A is the cross section area and the axial force, F, is 


calculated by: 
B= Abu’ (4.6) 


Milsaing from the bar equilibrium equation (1.2). Again, 


substituting and simplifying, Equation (4.5) becomes: 
GO =Eu’ (4.7) 


meal Substitution into Equation (4.1) results in an 


equation for total normal stress as follows: 
Cea Pia! +11.) (4.8) 


where Young’s Modulus of elasticity, E, is a function of 
material selection, the distance from the neutral axis to 
extreme fiber, c, is a function of cross section height, and 
u’ and v" are the first and second derivatives of axial and 
lateral displacements respectively. Using the Galerkin Finite 
Rlement Method, approximate values for the axial and lateral 
displacements can be determined at element endpoints. From 
these values, the stress distribution is computed and the 


optimization process can proceed. 


B. THE FINITE ELEMENT BEAM EQUATION DEVELOPMENT 
The Galerkin Finite Element Method (FEM) is an 


approximation method which transforms a linear differential 
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equation into a system of linear algebraic equations. Using 
the beam equilibrium equation (1.1), approximate lateral 
displacements for the arch can be determined at the system 
nodal points. For this method, a family of hermite cubic shape 
function which possess the Kronecker Delta property, are 
introduced in order to maintain the necessary function and 
Slope continuity for the fourth order beam equation. An 
approximate solution, v, for displacement, V, is formed as 


follows. 
Viet = Ou (4.9) 


where v is the exact solution of the beam equation in 
continuous space, ¥ is the approximate solution in discrete 
space, Q* is the transpose of a column vector of the cubic 
shape functions, and V is the vector of lateral displacements 
and slopes. 

After the approximation is formulated, the next step in 
the Galerkin method is to form the residual, R, in the 


following Formac; 
R=Q(7) = Os) (4.10) 


where p, is the lateral excitation force and @ denotes the 
differential operator which in the case of the beam 


equilibrium equation is defined by: 


(vy) = [EI(v!’) ]” (4.11) 


With substitution, the residual becomes: 
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eae One ale = py (S) 


From the residual, the Galerkin Equations are formed: 
forryas=0 (4.13) 
D 


meen as the null vector. Further substitution for R into 


the Galerkin vector equation results in 


fo ter(gty)”’ 1"ds - fop,(s) ds = 0 (4.14) 
D D 


To solve the Galerkin Equation, integration by parts is 


performed twice which yields: 


QL EI(Q*v)""1’|,-QEI(Q’v)""|, 


(4.15) 
/ Tos = = 
“fo pie ((@"v.) ‘ds Jae, (s) ds 0 


where |, denotes evaluation of these vectors at the boundary 
points of the structure. Recognizing that the lateral 


displacement and slope vector is constant, Equation 4.15 is 


rewritten as: 
(Ol aa (OR Ys \0 se ed B28 (Sel Ae 


(4.16) 
/ aR) = = 
fe Bio Cas. [o0,(s) ds 0 


From the beam equilibrium Equation (1.1), the shear, V, is 


defined by: 
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Y= oi (4.17) 
and Moment, M, by: 
M=EIv"’ (4.18) 
Thus, the boundary term load vectors are defined by: 
Vea eel) v1 (4.19a) 
and 
M=QEI(Q")"‘vl, (4.19b) 


Additionally, for convenience a system stiffness Mat rixjeas 


is defined by: 
Ke = [9 ET(Q")""as (4.19c) 
D 
and a system Force vector, E”, by: 


D 


Substitution of Equations (4.19 a through d) into BQuiaeaueg 
(4.16) results in the following system of linear algebraic 


equations: 
V| Ml a+K’v-P=0 (4.20) 


Further simplification is possible by defining F° as the load 


vector of internal and external applied lateral loads by: 
gees" +M| ad B (4.21) 


Thus, Equation (4.20) reduces to: 
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Kev=r* (4.22) 


where the global or system bending stiffness matrix, K*, is 
constructed from the union of all the elemental bending 
emeenecs matrices k* and the global bending force vector, 
F®, is constructed from the union of all the elemental bending 


merece vectors, Ff’, 


C. THE FINITE ELEMENT BAR EQUATION DEVELOPMENT 

In a similar manner to the beam equation, the Galerkin 
Finite Element Method is applied to the bar equilibrium 
equation (1.2) to approximate the axial displacements at the 
endpoints of a bar element. However, the bar equation is only 
a second order linear differential equation. Therefore, a 
family of linear shape functions which posses the Kronecker 
Delta property, are used in order to maintain the necessary 
eemte~eron continuity only. An approximate solution, a, for 


axial displacement, u, is formed as follows: 
uz=0=GTu (4.23) 


where u is the exact solution of the bar equation in 
continuous space, 0 is the approximate solution in discrete 
space, G’ is the transpose of a column vector of the linear 
shape functions, and u is the vector of axial displacements. 

After the approximation is formulated, the next step in 
the Galerkin method is to form the residual, R, in the 


hollowing format: 
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R=@(0) +p,(s) (4.24) 


where p, is the axial excitation force and @ denotes the 
differential operator which in the case of the bar equilibrium 


equation is defined by: 
Y(u) = [AE(u)’]! (4.25) 
With substitution, the residual becomes: 
R= [AE(G"u)’}’+p,(s) (4.26) 


From the residual, the Galerkin Equation is formed: 


for) eee (4.27) 
D 


where Q represents the null vector. Further substitutiomaneo 


the residual equation results in 
fetaz atu) ')/ds+ (ep, (s) ds=0 (4.28) 
D D 


Unlike the beam equation development, only single 
integration by parts is performed to solve the Galerkin 


EqQudiwon,. This resuisec seem: 
AEG (G"u)’| ,- {G’ [AB(G"u)Ids+ [Gp,(s) =9 (4.29) 
D D 


where |, represents evaluation at the boundaries of the 
Structure. Recognizing that the axial displacement vector is 


constant, Equation (4.29) is rewritten as: 
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G(AEG")'ul ,- [G! [AE(G™]dsu+ (Gp, (s) =0 (4.30) 
D 


meemecie bar equilibrium Equation (1.2), the axial force, F, 


Pemeaerined by: 
F=AEu’ (4.31) 
Thus, the boundary term load vectors are defined by: 
P=AEG(G")’ul, (4.32a) 


Additionally, for convenience a system stiffness Rake tesre al amgas = 


is defined by: 


Ki= G [AE(G"] ds (4.32b) 


D 


eemeetmoystem Force vector, F*, by: 


F?=|Gp,(s) (4.32c) 


D 


Substitution of Equations (4.32 a through d) into equation 
(4.30) results in the following system of linear algebraic 


equations: 
Bau + B=. (4.33) 


meeener Simplification is possible by defining F* as the load 


vector of internal and external applied lateral loads by: 
Pee (4.34) 


Thus, Equation (4.33) reduces to: 


Sel 


K*u=F4 (4.35) 


here the global or system axial stiffness matrix, K*, is 
constructed from the union of all the elemental axial 
stiffness matrices ce and the global axial force vectOu aa 
is constructed from the union of all the elemental axial force 


VECCOES peraa. 


D. THE ELEMENTAL STIFFNESS MATRIX 

The global Galerkin FEM Equations (4.22 and 4.35) are 
constructed from the union of elemental axial and bending 
stiffness matrices, Bs cline isa and axial and lateral force 
VECEOrs, tee and ries For the beam element, the elemental 
degrees of freedom in which the elemental forces act are shown 


in Figume”(422)2 


\ 2! 
LOCAL 
Nope 1 


Figure 4.2 Beam Element - Degrees of Freedom 
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iiieyethe stiffness matrix, k* for bending results in a 4X 4 
matrix of the form: 
ae ie fa 
hee Pee Peres are (4.36) 
| ks kay’ kaa Kae 
[kat kez kes kay 
For the bar element, the elemental degrees of freedom in 


which the elemental forces act are shown in Figure (4.3). 


Figure 4.3 Bar Element - Degrees of Freedom 


Thus the stiffness it~ U eb lb ea Sala for axial force results ina 


Men 2 Matrix of the form: 


5). 


pao ai 
ee Kr (4.37) 


= Ks kos 
To simplify, the elemental degrees of freedom are 


redefined for bar-beam elements as depicted in Figure (4.4). 


§ 
ow 


3! 


Figure 4.4 Bar-Beam Element - Degrees of Freedom 


This results in a combined 6 X 6 stiffness matrix, ee of the 


form: 
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Peer Oe ks 0 0 

O kif kiz 0 kag kis 
elec kxi k2p 0) kas kz (4.38) 
oreo Oke, C0 «0 

QO ky k3y O ka3 koe 

QO ka Kas 0 Kaz Kaa 


The elemental displacements and forces follow suit and are 


defined as follows: 


The elemental displacements vector, 6%, becomes: 


/ sll af i! f ef rw 
(07)7=< 61, pa gHO, 7 Os , One (4.39) 


where for the ith element 


6,7 = the axial displacement at local node 1 
oo = the lateral displacement at local node 1 
ena = the beam slope at local node 1 

- = the axial displacement at local node 2 
6,7 = the lateral displacement at local node 2 
a = the beam slope at local node 2 


The elemental force vector, f*, becomes: 


(fe =< fe lets be > (4.40) 


where for the ith element 


icp = the axial force at local node 1 
ea = the lateral force at local node 1 
ia = the moment at local node 1 

ical = the axial force at local node 2 
a = the lateral force at local node 2 
ie = the moment at local node 2 


Thus, the combination of the Galerkin Beam and Bar Equations 


for each element simplifies to: 


S5 


Porras (4.41) 
where the elemental stiffness matrix, k*, in terms of known 


quantities becomes: 


AE/®, 0 0 ~AE/t, 0 0 
O 12E1/ti3 6EI/¢ 0 -12EI/t 6E1/¢ 
pene O G6EI/ti? 4EI/t, O ~6EI/G 2EI/t, (4.42) 
- ~AE/€, 0 0 AE/¢; 0 0 
0 -12EI/t -6EI/G 0 12EI/t -6EI/G 
0 6EI/(G, 2EI/t; 0 ~-6EI/G 4EI/t, 


It should be noted that the bar and beam have uncoupled 
behavior. 
E. COORDINATE TRANSFORMATION OF THE ELEMENTAL SYSTEM 

OF EQUATIONS 

For curved structures such as the arch, each element has 
a unique orientation with respect to the global x and y axes. 
Therefore, to solve the global system of equations, the 
elemental Galerkin Equation (4.41) is transformed from local 
to global coordinates. The horizontal and vertical axes of the 
arch are chosen for a global reference coordinate system. 
Figure (4.5) depicts the angle the i'® element makes with the 
horizontal x-axis as @,, and the compliment angle, B,, as the 
angle the it" element makes with the vertical y-axis. 

From these definitions, the local displacements and 
forces, marked by a prime to indicate element degree of 
freedom are defined in terms of the reference coordinates axes 


as follows: 
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Figure 4.5 
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of 
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i™ ELEMENT 


Freedom 


and 


Accordingly, a transformation matrix, I°, for the ithvegemege 


becomes: 


which reduces the notation of Equations (4.45) and (4.46) to: 


and 


ll 


6,7 cos (G,) Gee scos (b,) 
=O, Cos {BJ eaemmces (7) 
8.3 

6, cos (G,) 480. ge ocamioL) 


=0,’ COs (0, JmeenO meCos micoe 


On 


ff," cOosmia,) Sue, Cocumpe 
-f, cosim( Jeet, "eosmcr 
tee 

fF, COS (Oi) tt. Cocmaioy, 


=f, cos (b;) 7 62 aeos 


ae: 


COS (Aco ci(in aay 

-cos(B,) Ges(a ame 
0 0 1 
0 0 Oe Coste.) 
0 0 Omecos (fp) 
0 0 Q 0 


§ 2! = 1'8! 
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cos (f,) 
eos (7) 
0 


Sola S&S S&S 


(4.45) 


(4.46) 


(4.47) 


(4.48) 


EE 





se Bl ve (4.49) 


where 
(87) 7=< 8:, 83, 53, 83, 85, 5; > (4.50) 


(Eo ya Ena usar ealt aa nae ree (4.51) 
Thus, the transformed elemental stiffness equation becomes: 
S10, Ge (4.51) 


Pyemoe@osticuting Equations (4.48) and (4.49) into Equation 


(4.41). 

Pyemultiplying both sides of Equation (4.52) with the 
inverse of the transformation matrix, lie, an orthogonal matrix 
meee) *=17), yields: 

(I) hae (T*) 8? a ae (4 753) 


where the elemental stiffness matrix, k’, in terms of the 


Global x and y coordinates is defined by: 
Hoe (nue (nt) (4.54) 


F. SOLUTION 


Recall from the Beam and Bar FEM development that the 
Global system of equations result from the union of the 


elemental stiffness matrices and force vectors such that: 


oo 


KA=F (4.55) 


where the global or system force vector, F, is the union of 
the transformed local force vectors, f°, and the Gl@baime 
system stiffness matrix, K, is the union of transformed ee 
stiffness matrices, xk’. Thus, Equation (4.55) is solved=ime 
the global displacement vector, A. 

These global horizontal, vertical, and rotational degrees 
of freedom are transformed back to local axial, lateral, and 
rotational displacements by the same transformation 
:elationships of section E (Equations 4.45 and 4.46) baaew 
these local displacements, the virtual loads at the element 
endpoints are computed from Equation (4.41): 

a0 = fi (4.56) 
where the elemental stiffness matrix, k’, is defined by 
Equation (4.42). 

The node point virtual loads, f*, equate to the virtual 
axial and lateral forces, and bending moments located at the 
endpoints of each element. From Equation (4.2) and (4.5), 
bending and axial stresses are calculated. For continuity, the 
stresses of internal global nodal points are averagedmeum-. 
physically, local nodal point 2 of the i*® element is the me 
point as local nodal point 1 of the it? + 1 element. Therefore, 
using Equation (4.1), the normal stresses can be determined 


for each global nodal pointe 
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V. PROGRAM DESCRIPTION AND VALIDATION 


mmome the development of Chapters II and IV, a VAX Fortran 
77 Code for FEM analysis of an arch was written to interface 
Regmene DOT SOftware package. The main program, ARCH OPT.FOR, 
and associated common program, ARCH COM.FOR, are contained in 
meeemearx C. Briefly, ARCH OPT.FOR opens and reads an input 
meee ncH IN.DAT, before it 1s divided into several 
Subroutines that perform the FEM analysis. Table 5.1 lists the 
Baemercata fields required of ARCH IN.DAT along with a brief 


@eeeription of each. 


TABLE 5.1 ARCH IN.DAT FIELD PARAMETERS 


Input File ANGLE, RADIUS, YOUNG, YIELD, 
form: NEL, METHOD, IPRINT, 
DV1BG, DV1LO, DV1UP, 
DV2BG, DV2LO, DV2UP, 
CLAN, FX, FY, FM, FA, 
OPTDCS, ITERATE, PRCSN, 
BX1, BY1, BM1, BX2, BY2, BM2, 
LABEL 


Parameter Description 


ANGLE The angle from 0 to 359 degrees subtended by the arch 
Seruceure. 


RADIUS The length of the arch radius of curvature. (The dimension 
is arbitrary, however all remaining inputs must be 
consistent.) 


Young’ Modulus of Elasticity for the arch material. 


The yield strength of the arch material. If a factor of 
safety is desired, it should be accounted for prior to 
aljpyohbie, c 


An integer number of elements, from 1 to 32, used to 
approximate the arch structure. 
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METHOD 


The optimizer method ten eesused. 


METHOD = 0 or 1: Modified Method of Feasible Directions 
METHOD = 2: Sequential Linear Programming 
IPRINT On screen print control parameter. Integers from 0 to 5 


DV_BG 


DV_UP 


CLAN 


FX 


FY 


FA 


OPTDCS 
OPTDCS 


OPTDCS 


ITERATE 


PRCSN 


PRCSN 


PRCSN 


indicate increasing screen printout. 


The best guess for design variable 1, the base dimension, 
or 2, the height dimension. Nodal point dimensions are 
initialized to the best guess value, thus establishes the 
optimizationmstarting point. 


The lower limit or side constraint for design variable l, 
the base dimension, or 2, the height dimension. 


The upper limit or side constraint for design variable l, 
the base dimension, or 2, the height dimension. 


An integer from 1 to NEL + 1 that indicates the node at 
which the concentrated load is to be applied. 


The magnitude of the concentrated load in the horizontal 
direction applied at node CLAN. 


The magnitude of the concentrated load in the vertical 
direction applied at node CLAN. 


The magnitude of the concentrated moment applied at node 
CLAN. 


The magnitude of the uniformly distributed load in the 
radial direction which spans the entire length of the arch. 


Optimization optiaen 


ae Optimize the dimensions of the problem. 


“a: Do not optimize the problem. Based on the initial 
design, calculate the stress distribution only. 


The number of iterations performed. The resulting optimized 
variables are re-entered into DOT and the optimization 
performed ITERATE times to effect an iteration. 


Computer precision used by the equation solver. 
= 2: Single precision 


= 2: double precision 
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BX_ Boundary conditions for horizontal displacement at 1, the 
Piast emocde @Lmeaerarech, mode 1,,0r 2, the last node of the 
aech, node NEL + L. 


BX_ = 0: Mneenodemsefree LO move horizontally. 
BX = 1: The node is not free to move horizontally. 
BY_ Boundary conditions for vertical displacement at 1, the 


first node of the arch, node il, or 2, the last node of the 
arch, node NEL + l. 


BY = 0: phe Nodemuce tree toemove vertically: 
Boe = 1: The node is not free to move vertically. 
BM _ Boundary conditions for the beam slope at 1, the first node 
QithewaLcn, mode 1, or 2, the last node of the arch, node 
Neo + 1. 
BM = 0: The node is free to rotate. 
BM = 1: The node is not free to rotate. 
LABEL A character string used to identify the output. 


Poeeuelined in Figure 5.1, the main program, ARCH OPT.FOR 
is divided into subroutines. In general, subroutine Geometry 
1s called in order to generate the x and y coordinates of the 
global nodal points and determine the orientation of each 
meememe. Following Geometry, subroutine Optimization tool 
establishes the DOT parameters prior to the first call of the 
DOT program. The first call serves only to record the DOT 
parameters selected in DOT’s internal arrays. After DOT is 
Meeeea, the Optimization tool subroutine, calls Eval to 
Memeiece the objective function and constraint functions 
Originally outlined in the problem formulation of Chapter II. 

As detailed in Chapter IV, the constraint functions are 


made functions of the design variables through Finite Element 
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DATA FILES MAIN PROGRAM OPTIMIZATION LOOP = EXTERNAL SOFTWARE 


| ARCH_OPT.FOR =| 
| 
| 
| 

ee 

| | 





Figure 5.1 Arch Opt Program Structure 


Mebnoa analysis. Subroutines Form and Force vector develop the 
global stiffness matrix and force vector, which are modified 
by subroutine Bndary for the appropriate boundary conditions. 
The equation solver, L2ARG, from the IMSL library is called to 
solve for the global displacements, which in turn are used to 
calculate the nodal stresses. Once the constraints are 
evaluated for the initial design, the problem is returned to 
DOT where the move direction is computed and an updated design 
Moume, chosen. The objective and constraint functions are 
reevaluated for the updated design point before returning to 
Beeeror further iteration. 

Once termination criteria for optimization are reached, 
Biemiealim program creates the output file, ARCH OUT.DAT. This 
file contains the problem parameters, optimized design 
variables, and the resulting objective function value along 
Seema variety of additional information. Summarizing, for a 
Given geometry, loading, and set of end conditions, the 
jm@eegram iS capable of finding the optimum cross section 
dimensions of each nodal point along the length of the arch. 

To validate the FEM analysis, several non-optimum straight 
beam and arch problems with known analytical solutions were 
solved. A straight cantilever structure, subject to a 
concentrated lateral end load, axial load, and end moment; and 
a quarter cantilever arch, subject to a lateral end load were 
analyzed. These test problems established the program error 


for stress and displacement calculations. Additionally, the 
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quarter cantilever arcn and a hinged-hinged semi-circular ag 
structure, subject to a lateral load on the axis of symmetm@ 
establish trends in a relationship between the number 

elements used to approximate the arch and accuracy. T 
remainder of this chapter is a summary of the results and t 
conclusions drawn from each validation problem studied. T. 


complete solution of each problem is contained in Appendix | 


A. VALIDATION I: CANTILEVER BEAM 
A cantilever beam is subject to a concentrated end load | 


Shown in Figure 5.2. (Ref. 8] 


5 xa 0.0 
L=45.00 inches 


4 1 25° B=1.50 inches 
H=3.00 inches 





3 
22 5" ta = 3.375 inches! 
P=1,000.01bf 
a 
37.75 Bee rai 
ae 
mF 
» 0 eS 
45, £ 


Figure 5.2 Validation Case #1 —_ 


ARCH OPT.FOR was run for this beam stsucturecesicee 


angle of 45.0 X 107° radians and a radius of 10° inches J 
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approximate a straight beam of 45 incnes. The four element FEM 
Peet Om 1s COMpared £6 the analytical solution in Table 5.2. 


TABLE 5.2 


= ae | 
STRESS STRESS 

| 4 | _sooo.0] 4999.9] 0.0020 | 

‘Saeed 0.0 0 


0.0000 | 
















wnere the percent error is defined as: 


$ Error = ( Theory - FEM Analysis )/ Theory * 100 


B. VALIDATION If: PRISMATIC BAR 
Similarly, a prismatic bar 1s subject to an axial load as 


Seow In Figure 5.3. {Rer. 8} 


ZL=45.00 unches 
B=1.50 inches 
H=3.00 inches 

A=B: H=4.50 inches? 
P=1,000.0 ibf 


B= 2S Wee wea 
a 


a 





Figure 5.3 Validation Case #2 


Input values for angle and radius remained the same to 
approximate the straight bar. The four element FEM solucwenmae.e 


compared “iielac temo. 





ee a ABLE 5.5 
sang la FEM ANALYTICAL | %* ERROR 
= ae ae aia 
a eee Z eT eee Z | 9.0000. 0000 






+ 222.2 222.2| 0.0000 | 
a3 aoe 22 D2 222. 25| OF 00On 
meee SODED 222.2| 0.0000 





C. VALIDATION III: CANTILEVER BEAM 
The cantilever beam is subject to a concentrated moment at 


the free end as shown in Figure 5.4. [Ref. 8] 


M L=45.00 inenes 
5 ) 45." B=1.50 inches 
H=3.00 inches 


BH? 
2 


E=30 x10° pen 
M=10,000lbf 


I= =3.375 inches 


Mx 
Slope S Ey 





Displacement S= Mx* 
2EL 


Figure 5.4 Validation Case #3 
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The four element FEM solution for both slope and 
displacement is compared in Table 5.4. 


TABLE 5.4 


2 THEORETICAL | FEM ANALYTICAL | % ERROR | 
SLOPE SLOPE | 


[9.000000] 0.00000000| 0.0000 | 
[2 [0.00111111 | 000321109] 0.0019 | 
[4 [000333333] 0. 00333328] 0.0016 















ODES el 
DISPLACEMENT DISPLACEMENT 
om a Oa0e S50 00m 0.02499940 0.0024 


| 4 | 0.05625000 0.05624880 0.0021 
0.10000000 0.09999790 0.0021 





D. DATION IV: CANTILEVER QUARTER ARCH 


A cantilever quarter arch is subject to a lateral load as 


Seeeeeowm in Figure 5.5¢°(Ref. 8] 
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Pp L= 45.00 inches 
B=1.50 inches 
H=3.00 inches 


3 
| I= 5 =3.375 inches‘ 


x E=30x10° psi 
P= 1700070 Lbs 


5 = PR? 
\ Mosse 


Figure 5.5 Validation Case #4 





ARCH OPT.FOR was run for this Structure using an aueomen 
90.0 degrees and a radius of 45 inches. To approximate the 
arch a four, six, eight, ten and 12 element FEM solution is 
solved and compared to the analytical solution presented in 
Table 3.55 


TABLE 5.5 





THEORETICAL | FEM ANALYTICAL | % ERROR 
Ox Ox 


ees | 

a] 0.450000 0.446951] 0.677556 

a | 0.450000 0.448382} 0.359556 
i 






0.450000 0.448790 0.254667 | 
0.450000 0.448790 DE 26 cc 








0.450000 0.4491009"  On200000 
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E. VALIDATION V: HINGED-HINGED SEMI-CIRCULAR ARCH 
A hinged-hinged semi-circular arch structure is subject to 


a lateral load along the axis of symmetry as shown in Figure 


Seo. (Ref. 8] 


R=32.0inches 
B=1.50 inches 
H=3.00 inches 
2) GO) 6) Boke 


_ PR. _ _ PR 
M = (1 cos@) — (8m0) 


Figure 5.6 Validation Case #5 


Results are tabulated in Table 5.6 for comparison of the 
four, Six, eight, ten, 12, 14, and 16 element FEM solutions to 
jee analytical solution. It should be noted that by using 
symmetry, the arch structure iS approximated by twice the 


number of elements shown in the calculations. 


ae 


TABLE 5.6 


THEORETICAL | FEM ANALYTICAL % ERROR 
STRESS STRESS 
| 4 





0 eh) ear se 
22.5 = | 119 a i123 
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8 oof lw OL O00) 


pant 

p.ai.25 [| 7465.5] 7509.3] 0.586658 
| 22.5 
33.75 
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F. CONCLUSIONS 

The four element approximation for a straight cantilever 
structure produced an error mo greater than 0.0l@2eae 
cantilever quarter arch produced an error less than 0.70% for 
the four element model which reduced to less than 0.20% with 
12 elements. The results of the hinged-hinged arch indicate, 
as expected, that the more elements used the better the 
solution. Considering only meaningful stresses, stresses in 
excess of 10,000 psi, the error is less than 2% for eight 
elements and less than 0.8% for 12 elements. 

In general, the percent error recorded for the first three 
validation cases seemed Tig es aCe pgiaL te Aer QUE.» Four element 
approximations sufficed to solve the stresses, slopes, and 
displacements for straight structures. Therefore, it was 
concluded that the program was producing accurate results for 
analysis of straight beams. 

Unfortunately, for the arch structures, the error mae 
four element model was significant (greater than 45%). 
However, the error reduced significantly when more elements 
were used to approximate the structure. Grid independence, (2% 
error), was not achieved for the hinged-hinged arch until at 
least eight elements are used to approximate the structure. 
This indicates that an element cannot be used to span more 
than 11.25 degrees of arch. The resulting trend, as expected, 
confirms that the more elements used, the better the model. 


However, computer time and computer error increase with 
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Mmiiverease in the number of elements and models of more than 


eight elements were not used. 


SS) 


VI. CASE STUDIES 


Results are presented for a number of cases with regard to 
optimization scheme and stress analysis. The case studies 
range from the simple cantilever beam to complex arch struc- 
tures. In addition, for many cases one parameter of the same 
structure was modified and the problem was reoptimized to 
establish a comparison. The straight beam is examined first, 
followed by five cases studying the quarter cantilever arch 
with varied loadings. Cases #7 and #8 are symmetric 
semicircular arches comparing simply supported arch structures 
with fixed-end arch structures. The remaining cases are 
asymmetric semicircular arch structures. Cases #9 through #11 
investigate various end conditions and Cases #12 through #14 
various combined loadings. The cases conclude with Case #15 
which combines a concentrated lateral load, applied moment, 
and distributed load across the arch structure. 

For each case, interpretations of the results 
accompanied by a schematic drawing of the structure modeled, 
a plot of the cross section dimensions and area as functions 


of nodal points, and a plot of the axial and bending stresses 
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as functions of nodal points. Eight elements were selected to 
model the arch structures and the material properties were 
selected such that the yield strength was imputed as 52,000 
psi and Young’s Modulus as 30,000,000 psi. For reference, the 
Modified Method of Feasible Directions will be referred to as 
Method 1 and the Sequential Linear Programming Method will be 
referred to as Method 2. Additionally, each endpoint, unless 
geometrically restricted by imposed boundary conditions, can 
have three ‘means of displacement,’ MOD. An endpoint can 
meme amout the z-axis, displace in the x direction, and 
@mempmace 1n the y direction. For reference, an endpoint will 
be described by a number from zero to three reflecting the 
means of diplacement. As an example, a fixed end is considered 
to have zero means of displacement because it cannot rotate or 
displace in either the x or y direction. A free end which can 
rotate and displace in both the x and y direction is 
considered to have three means of displacement. A hinge which 
can only rotate has one MOD. The complete computer data 


printout is presented in Appendix D. 


Dal 


A. CASE #1: CANTILEVER BEAM WITH LATERAL LOADING 

The cantilever beam was optimized first in order to 
provide guidance for adjusting the various parameters 
discussed in Chapter III. Satisfactory results were produced 
by turning the auto scaling function off, reducing @C iaaw. 
CTMIN, and establishing the termination criteria. Using the 
Modified Method of Feasible Directions, henceforth referred to 
as Method 1, the cross section dimensions and stresses were 
plotted. As expected, the dimensions form a parabolic function 
over the length of the beam. Furthermore, the beam exhibits 
only stress due to bending moments. The normal stresses are 


virtually nonexistent which likewise is as expected. 


Loads 


2000* Lateral = 2,060 Qs 
Axial = 0 ibs 
Moment = C in-lbs 
End conditions 

Node 1 0 MOD 

Noce 9 3 MOD 
Dimensions 

f 


y 32 Radius 


Theta 


not applicable 
not applicable 


Total volume 


a Volume = 33.13 Sine 
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Case #1 





(ui) UoisueWwig 


Node 


—+— Height —— Area (in* 2) 


; ~#- Base 





5 6 
Node 
—=a— Axial —+— Bending 





Sage Chit clei: 


r 7 ; 


9189S DO7- (isd) ssens 
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B. CASE #1A: CANTILEVER BEAM WITH LATERAL LOADING 

For comparison with Case #1, the same cantilever beam was 
optimized using the Sequential Linear Programming Method, 
henceforth referred to as Method 2. The results are quite 
Similar. In total structure volume, the difference is less 
than 0.07%. The only Significant difference appears at nodal 
point 9, the free end. In theory, the free end of a beam can 
support no bending stresses. For this case, nodal point 9 has 
no stresses unlike the previous case which had relatively 
small bending stresses at nodal point 9. However, from this 
result alone it iS not concluSive that Method 2 iS Superior to 


Method 1. 


Loads 
Lateral = 2,000 lbs 
2000 # Netlog zs 0 lbs 
Moment = 0 in-lbs 
End conditions 
Node 1 0 MOD 
y Noae 9 3 MOD 
Dimensions 
Radius - not applicable 
32 Theta = not applicable 
Total volume 
Volume = 33.156 
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Dimension (in) 


Case #1a 


a 


Stress (psi) - LOG Scale 









2 3 4 5 6 7 8 
Node 
| -™#- Base —— Height —— Area (in~2) 
a 


ae ee | 


10000: 
i ; <—e 
| | 
1000: |. oe 
5 - 4 
: : a7 
100-1 
eee Ss eae se 6 87 gs 
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C. CASE #2: CANTILEVER ARCH WITH LATERAL LOADING 

Since Case #1 did not strongly suggest a preferential 
method, the cantilever arch was optimized with Method 1. At 
most nodal points, the total stresses were well below the 
yield stress which indicates that this design is far from an 
optimum Strierune.” Addictions the height and base 
dimensions hovered around the initial starting pow wee 
inches by 2 inches and produced a Structure only 7.42% less in 
volume than that of the initial Structure. It appears that the 


optimizer failed to achieve an optimum solution using this 


Method. 
Loads 
Pareral = 2,000 lbs 
Rxval = 0 lbs 
2000 End conditions 
x Node 1 C MOD 
Node 9 1 MOD 
Dimensions 
: Racius = 32 in 
preos = 90 degrees 
Total volume 
32” Volume = 186.15 ein. 
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Case #2 
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D. CASE #2A: CANTILEVER ARCH WITH LATERAL LOADING 

For comparison, the same arch structure was reoptimized 
with Method 2. Each element of the structure now supports 
Stresses equal to the yield stress producing an efficient 
structure. The total volume was reduced from the initial 
starting point by 61.32%. For this structure, Metheameae. 
produced results with fewer iterations than Method 1. With 
these observations in mind, Method 2 was selected as the 
preferred method for quarter arches. Additionally, it is 
interesting to note that the axial stresses only remotely 
effect the stress total for the first 5 nodal pointsjpee@ene- 
the first 45 degrees of arch. After node 5, the height reduces 


Significantly, however the area remains roughly the same. 


Loads 
Lateral = 2, 000 Ibs 
Axial = 0 lbs 
Moment ~ 0 in-lbs 
End conditions 
Node 1 Q MOD 
Node 9 3 MoD 
Dimensions 
Radius = 32 an 
Theta = 90 degrees 
Total volume 
Volume = 77.78 in? 
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Case #2a 
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E. CASE #3: CANTILEVER ARCH WITH AXIAL LOADING 

This case presents a quarter arch structure subject to an 
axial load vice the lateral load of Case #2. Unlike Case #2 
the axial stresses increase significantly along the length of 
the arch and the bending stresses decrease. The net result is 
an arch structure of 27.15% less material. This seems to 
indicate a dominant relationship between area and bending 
stress. Additionally, this case exemplifies the difficulty 
experienced by approximating an arch of 90 degrees with eight 
Straight segments. The plots appear very disjointed, hence the 
data points seem circumspect. However, the effect can be 


minimized as presented in Case #3a. 


Loads 
Lateral = 0 lbs 
y Axial = 2,000 lbs 
Moment = 0 in-lbs 
End conditions 
2000” Node 1 0 MOD 
xX Node 9 3 MOD 
8 Dimensions 
Radius = ae in 
Theta = 90 degrees 
Total volume 
ae Volume = 56.66 in 
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Case #3 





(u!) UoisueWIGg 


Node 


| 
| 


—*— Area (in~ 2) 


—+— Height 


| -#- Base 





{ { we 
: ae ar a 
ce nee ‘ + ‘_™ co 
7 | a ' 
: | aN 
ne 6 2 
' I 
ie AD wo 
@® 
bhe Ww 5 
za 
eye wv 
op] 
N 
ae 
Tew feito Vee UN atae Tee Sine TTT a oe Oh ek cee aa 





a i. a: 


Qj29S DO1- (isd) ssens 


67 


F. CASE #3A: CANTILEVER ARCH WITH LATERAL LOADING 

Thus far, each case has started with an initial design of 
2 inches by 2 inches. For comparison, the arch structure of 
Case #3 was optimized a second timeusing the results of Case 
#3 as the initial design. Reoptimizing had the desired effect 
of smoothing the results and in graphical form, both the area 
and stress curves take on a fairer shape. In terms of total 
structure volume, the reoptimized arch was 27.03% smaller than 
that of Case #3. In all subsequent cases this strategy of 


reoptimization will be referred to as a two-stage optimization 


SEGaAECOy. 
Loads 
Lateral = O Tee 
| Axial = 2,000 lbs 
Moment = O in-lbs 
- End conditions 
x 2000 Node 1 O MOD 
— ~ Node 9 3 MOD 
Dimensions 
} Radius = 32 an 
Theta = 90 degrees 
Total volume 
32" Volume - 41.34) ine 
ee 
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G. CASE #4: CANTILEVER ARCH WITH LATERAL LOADING AND MOMENT 

For this case, a lateral load and concentrated moment were 
applied at nodal point 9. The shape of the dimension plot 
curves are very Similar to those of the cantilever beam, 
parabolic in form. In comparison with the same structure 
subject only to the lateral load, Case #2, the total structure 
volume is reduced by 18.67%. The concentrated end moment 
negates the effect of lateral load on the extreme fibers by 
producing compressive stresses on the outer fibers and tensile 
stresses on the inner fibers of the arch. Thus, the cross 
sectional dimensions necessary to withstand the total normal 


stress 1s reduced thereby reducing the total structure volume. 


Lateral = 2, 000=lss 
Ps Axial = O lbs 
iy 24000 Moment = 1,000 in-lbs 
End conditicns 
Noce 1 Oi MSD 
Noace 9 3 MOD 
Dimensions 
Radius = a in 
Theta = 90 degrees 


Total volume 


Volume = 63.25 ame 
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H. CASE #4A: CANTILEVER ARCH WITH LATERAL LOADING AND MOMENT 

To further emphasize the effect of the concentrated end 
moment, the structure of Case #4 was subject to the same 
lateral load while the moment at the end point was increased 
by a factor of 10. By increasing the applied moment, the 
effect of the lateral load on the extreme fibers is negated 
further which reduces the cross sectional area necessary to 
withstand the total stresses. Expectedly, the volume reduced 
from Case #4 by 15.88% for a total reduction from Case #2 of 
31.60%. It iS interesting to note that the shape of the 


dimension curves still remain parabolic in form. 


Loads 


Lateral = 2,000 lbs 
Axial == O es 
Moment = 10,000 in-lbs 
End conditions 
Node 1 C MOD 
Noage 9 3 MOD 
Dimensions 
Radius = 322i 
Theta = 90 degrees 
Total volume 
Volume = 5 2 iin 
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Dimension (in) 


Stress (psi) - LOG Scale 
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I. CASE #5: CANTILEVER ARCH WITH DISTRIBUTED LOADING 

This case is presented to display some of the versatility 
of the program. A load acting radially inward is distributed 
along the length of the arch. The cross section dimensions and 
area curves appear to be almost linear and the bending 
stresses dominate the total stresses. Since the bending stress 
is a function of height Squared, the optimizer “trveduwae 
maximize the height dimension until the geometric constraint 
was violated. At each nodal point, the height is 10 times the 
Size of the base except at the end point for which both 
dimensions reach the minimum side constraint. Had the arch 
Structure not been optimized, the volume necessary to Support 


the distributed load would increase by 225%. 


Loads 


Lateral = 0 lbs 
4 Axial = 0 lbs 
Moment = 0 in-lbs 
Distrib. = 1,000 Ibs/in 
rm 100 LBS/jn, = | 
Ena conditions 
Node 1] 0 MOD 
Node 2 3 MOD 
Dimensions 
Radius = 32 
Theta - 90 degrees 


Total volume 


Volume = 55.70 oe 
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Case #5 





(uj) UOjSUBWIG 


—— Height —— Area (in* 2) 


, -—#-— Base 


net Renee eH ce) Sie SB pa 
H ! i | 
t, va i! ‘ ’ 
eter tj ————T4 
bes ay respect ‘ ‘ eee 
=. Sn 
aS tet 
i ; joe 
alae ae ntte 2 
tt 
Vy er 
YM ae 
ae jf ee 
i : ; 
i ial 
aa a 
ties oe oe 
es | 
od J 
' ioe 
Ti t ORS 
i ‘| 3 
i 
wie 
F Il - 
\ ee 





tok Met ot 
vv" 


aledS 50 1- (isd) ssans 





fas 


J. CASE #6: CANTILEVER ARCH WITH LATERAL AND 
DISTRIBUTED LOADING 


To build on Case #5, a lateral load was applied at the end 
point in addition to the distributed load. In comparison gene 
volume required to withstand the lateral load only is 77.78 


: (Case uma ame The volume required to withstand the 


in 
distributed load only is 55.70 in*. Yet the wvodumeuae 
withstand both the lateral load and the distributed load 
presented in this case is 97.47 in’. By combining loads which 
produce opposing bending moments, the volume of the resultant 
optimized arch is not equal to the sum of the volume of arches 
optimized subject to the individual loads. Therefore, it is 


possible to achieve a more efficient structure through 


resourceful combination loadings. 


Loads 


2,000 tbs 


Lateral 


Axial = 0 -Ibs 
Moment = O in-lbs 
Diserie. ze 100 1bs7an 
End conditions 
Node l O MOD 
Node 9 3 MOD 
Dimensions 
Radius = 32 328 
Theta = 90 degrees 
Total volume 
Volume = 97.47 in? 
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K. CASE #6A: CANTILEVER ARCH WITH LATERAL AND 
DISTRIBUTED LOADING 


For comparison, the same structure (Case #6) was optimized 
with the DOT auto scaling function switched on. Changing this 
parameter seemed to have little effect on the overall volume 
indicated by an increase by only 4.40%. However, the 
computation effort judged by total computer time nearly 
doubled and both the dimension and stress curves have 
unexpected behavior near the endpoint. This comparison 
confirmed that better results were achieved by switching the 


auto scaling function off for these structures. 


Loads 


Lateral . 2,000" lbs 
Axial = QO Ibs 
Moment = Q in-lbs 
Distrib: = 100 lbs/in 
End conditions 
Node 1 Q MOD 
Node 9 3a Mos 
Dimensions 
Radius = 327i 
Tneta = 90 degrees 
Total volume 
Volume = 101.76 in? 
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L. CASE #7: HINGED-HINGED ARCH WITH LATERAL LOADING 

For the remaining cases, it was observed that the only 
reliable and consistent results were obtained by using Method 
1 for optimization. It 1S theorized that reStiaseeeeem 
displacements at both endpoints may have caused Method 2 to 
become mathematically unstable and therefore unsuitable to 
solve such problem. For thiS particular case, meee 
interesting to note that at the base, node 1, and 56.25 
degrees from the base, node 6, the axial stress completely 
dominates the total stresses because there is virtually no 
binding force. At these points, the dimensions of the cross 
section, dictated strictly by the axial stress, form a square 


to produce the minimum area. 


Lateral = 16,000 lbs 

Axial = O ibs 

Moment = 0 in-lbs 
End conditions 


Node ]1 1 MOD 
Node 9 1 MOD 


Dimensions 


Radius = S32. ag 

Theta = 180 degrees 
Total volume 

Volume ~ 129.12 aa 
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M. CASE #7A: HINGE-HINGED ARCH WITH LATERAL LOADING 

Similar to Case #3a, the arch structure of Case #7 was 
reoptimized using the results achieved as the initial design 
in order to apply the two-stage optimization strategy. Again, 
reoptimizing had the desired effect of smoothing the results, 
however, this effect was not as dramatic for Method 1 as for 
Method 2. The two-stage optimization strategy only reduced the 
volume by 4.47% using Method 1 as opposed to the 27.03% 
reduction using Method 2. Additionally, at node 6, §eieueeea 
stresses exceeded the yield stress by 2.38%. Fortunately, this 
occurrence did not repeat in any other cases due to 
reoptimization. Therefore, the two-stage optimization strategy 


was applied for the remaining cases. 


Loads 


16,000 lbs 
O lbs 
O in-lbs 


Lateral 
Axial 
Moment 


End conditions 


Node 1 1 MOD 
Node 9 1 MOD 
Dimensions 
Radius = 32.328 
Theta = 180 degrees 


Total volume 


Volume = 123.35 eee 
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Case #7a 
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N. CASE #8: FIXED-FIXED ARCH WITH LATERAL LOADING 

For this case, the same loading of Case #7 was applied to 
a semicircular arch with fixed end points. This produces a 
Statically indeterminate structure with zero means of 
displacement at both the boundaries. As a result, the peaks of 
the axial stress curve are dampened and shifted towards the 
center by approximately 15 degrees. A larger bending moment is 
produced at the base since it 1S no longer free to rotate. 
However, the net results is that the total structure volume of 
Case #7 is reduced by 14.08% by changing the end conditions 
from Simply-Supported to fixed. From this, aS expecmecuaea 
structure more statically indeterminate results in a more 


efficient structure. 


Loacs 


16,000506s 
O° ibs 
O in-lbs 


bateral 
Axial 
Moment 


End conditions 


Noae 1 O MCD 
Node 9 O MOD 


Dimensions 


s2 in 
180 degrees 


Radius 
Theta 


Total volume 


Volume = 105.98 in? 
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©. CASE #8A: FIXED-FIXED ARCH WITH LATERAL LOADING 

In Chapter I, thesis research performed by Scott McDavid 
was mentioned as the predecessor for this investigation. His 
results optimized a fixed-fixed arch subject to the same 
lateral load with respect to the base dimension only. By 
holding the height dimension constant, the structure must have 
twice the volume in order to withstand the loading. For most 
arch structures, the bending stress is usually the more 
dominate stress. Therefore, the height dimension has more 
effect on the total volume than the base dimension because the 
bending stress is a function of base times height squared. 
When the bending stress dominates, the optimizer will seek to 


maximize the height. 


Lateral = 
Axial = O lbs 
Moment = 


End conditions 


Node 1 O MOD 

Node 9 O MOD 
Dimensions 

Radius 32 in 


Theta 180 degrees 
Total volume 


Volume = 218.77 in’ 
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P. CASE #9: FIXED-HINGED ROLLER ARCH WITH AXIAL LOADING 
Unable to invoke symmetry on the remaining cases, the 
elements used to model an asymmetric semicircular arch must 
span the full 180 degrees. The largest number of elements used 
to model the structure and produce consistent results remained 
only eight. It iS again suggested that restricting 
displacements at both endpoints cause the problem to become 
mathematically unstable. For this particular case, the arch 
has zero means of displacement at node 1 and two means of 
displacement at node 9. As expected, the arch is quite large 
at node 1 to support the resultant moment. At node 3, 45 
degrees up from node 1, the axial stress dominates the total 
stress and the size decreases. For reference, this arch is 


more than twice the volume of the arches in Case #7 and #8. 


Loads 


Lateral = 0 lbs 
Axial = 16,000 lbs 
Y Moment = 0 in-lbs 
End conditions 
Node 1 0 MOD 
Node 9 2 MOD 
Dimensions 
Radius = 32 en 
Theta = 180 degrees 
Total volume 
Volume = 287.15 eine 
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Case #9 
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QO. CASE #9A: HINGED-HINGED ROLLER ARCH WITH AXIAL LOADING 
To emphasis the conclusion drawn from case #9 about the 
endpoints, the arch structure and loading studied for Case #9 
was modified by adding an additional means of displacement at 
node 1. Allowing nodal point 1 to rotate freely, the dimension 
and stress curves alter drastically. The total structure 
volume increased by 19.92%, yet the structure cannot withstand 
the stresses. The total stresses exceed the yield stresses by 
54.81% resulting in an infeasible design. It appears that the 
optimizer failed to achieve an optimal solution for this arch 


Structure due to the additional means of displacement. 


Loads 


Lateral = 2,000 lbs 
Axial = 16,000 tbs 
y Moment = C in-lbs 
End conditions 
Node l 1 MOD 
Node 9 2 MOD 
Dimensions 
Radius = 32.07 
Theta = 180 degrees 
Total volume 
Volume = 344.34 ane 
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R. CASE #10: FIXED-FIXED ROLLER ARCH WITH AXIAL LOADING 

For this case, the same arch structure and loading studied 
in Case #9 was modified by reducing one means of displacement 
at node point 9. This produces a more redundant structure with 
a resultant decrease in total volume of 10.64%. In comparison, 
Case #8 with one less degree of freedom than Case #7 at both 
node 1 and node 9, had a reduction in total volume of 14.08%. 
Again, it is suggested that a structure more statically 
indeterminate results in a more efficient structure. 
Additionally, it is noted that when the axial stress dominates 
the total stresses, the area is reduced significantly and the 


cross section dimensions reduce to form a square. 


Loads 


Lateral = O lbs 
y Axial = 16,000 lbs 
Moment = O in-lbs 
End conditions 
x Node l 0 MOD 
Node 9 1 MOD 
i6000 * Dimensions 


32 on 
180 degrees 


Radius 
Theta 


Total volume 


Volume = 256.6102 n 
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S. CASE #11: HINGED-FIXED ARCH WITH LATERAL LOADING 
AND MOMENT 


To investigate the possibility that dominant axial 
Stresses result in volume reduction, a semicircular arch with 
one degree of freedom at node 1 and zero degrees of freedom at 
node 9 was subjected to a lateral load and applied bending 
moment. From this, it appears that the cross sectional area is 
inversely proportional to the axial stresses. Additionally, it 
appears that the dimension and stress curves of the left half 
of the structure behaves exactly as those of Case #7 which has 
the identical end conditions. Similarly, the curves of the 
right half of the arch behaves exactly as those of Case #8. 
This suggests that the boundary conditions do not effect the 


Structure past Jehewn apo. 


Loads 


Lateral 12,000 lbs 


Axial = O lbs 
y Moment = 1,000 in-lbs 
End conditions 

Node 1] 1 MOD 

Node 9 O MOD 
Dimensions 

Radius = 32 in 

Theta = 180 degrees 
Total volume 

Volume = 153.08 im: 
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T. CASE #11A: HINGED-FIXED ARCH WITH LATERAL LOADING 
AND MOMENT 


In order to test the possibility that dominant werwes 
stresses might reduce the cross section area and hence reduce 
the total structure volume, the structure of Case #11 was 
subject to the same bending moment while the lateral load was 
increased by a factor of 2. As a result, the axial stresses 
increased overall. Again, it appears that the cross sectional 
area 1S inversely proportional to the axial stresses. The 
dimension and stress curves displayed the same shape as noted 
before but the total volume increased by 57.95%. Therefore, it 
was concluded that increasing axial stresses may reduce the 
cross sectional area at specific nodes but the overall 


structure volume 1S neopereduced= 


Loads 


24,000 lbs 
0 lbs 
1,000 in-lbs 


Lateral 
y Axial 
Moment 


End conditions 


Node 1 1 MOD 

Noae 9 0 MOD 
Dimensions 

Radius 32 in 


Theta 180 degrees 
Total volume 


Volume = 241.78 in? 
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U. CASE #12: FIXED-FIXED ARCH WITH MULTIPLE LOADING 

To demonstrate further versatility of this program, a 
fixed-fixed arch was subjected to a combination load applied 
at an angle 45 degrees up from node 9. The load consisted of 
a concentrated lateral and axial load and an applied bending 
moment. As anticipated, there is a jump in the dimension 
curves at node 7 were the load was applied. Interestingly, 
22.5 degrees from each endpoint, the axial stress dominates 
and according a, the cross sectional area reduces 
significantly. Additionally, at node 5, the midpoine yeoman 
arch structure, the cross sectional area is significantly 


smaller as a result of an increase in the axial Stress. 


Loads 


Lateral = 17,000 ss 
Axial = 9,000 lbs 
4 Moment = 1,000 in-lbs 
End conditions 
x Node 1 0 MOD 
wk Node 9 O MOD 
Dimensions 
Radius = 32 ae 
Theta = 180 degrees 
Total volume 
Volume = 156.5500 n- 
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V. CASE #13: FIXED-FIXED ROLLER ARCH WITH AXIAL LOADING 
AND MOMENT 


For this case, the same structure of Case #10 is subjected 
to an equivalent axial load with an additional bending moment 
applied at nodal point 6, 67.5 degrees up from node 9. 
Shifting the load by 22.5 degrees and adding the applied 
bending moment appeared to have little effect on the overall 
design. In fact, the volume is increased from Case #10 by only 
3.65% and the dimension curves exhibit very similar 
characteristics. However, the dip observed previously in the 
axial stress curve at node 6 of Case #10 is not present in the 


axial Stress CULVeouotee eee cee 


Loads 


0 lbs 
16,000 lbs 
1,000 in-lbs 


Lateral 


4 Axial 
Moment 


tor ol 


Ena conditions 


Node i 0 MOD 
Node 9 1 MOD 
Dimensions 
Radius = 32s 
Theta = 180 degrees 


Total volume 


Volume ~ 265.96 in? 
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W. CASE #13A: FIXED-FIXED ROLLER ARCH WITH AXIAL LOADING 
AND MOMENT 


For comparison, the same structure and loading were 
reoptimized starting from a different initial design point. 
Previously, each optimization began from an initial design of 
2 inches by 2 inches at each node. For this case, the base 
dimension at each nodal point was 0.5 inches, and the height 
dimension at each nodal point was 3.5 inches. Incredibly, the 
volume of the resultant Structure is 33.96% less than the 
volume of the structure optimized in Case #13. Obviously, 
optimization can be a function of the Starting@ieonae 
Fortunately in the previous cases, various initial design 


points were tested and this occurence did not repeat itself. 


Loads 
Lateral = O ies 
4 Axial = 16,000 lbs 
Moment = 1,000 in-lbs 
End conditions 
Node 1 0 MOD 
Node 9 1 MOD 
Dimensions 
Radius = 32-18 
Theta = 180 degrees 
Total volume 
Volume = 175.65 in. 
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X. CASE #14: FIXED-HINGED ARCH WITH MULTIPLE LOADING 

In contrast to Case #13, this asymmetric arch was loaded 
at an angle on the side of the arch with zero means of 
displacement at the endpoint. The behavior exhibited by the 
dimension and stress curves was similar to that of Case #12 
which had zero means of displacement at both endpoints. 
Allowing for the difference in the magnitude and direction of 
the load, the only significant difference between Case #12 and 
this case appears at node 8 and 9. It is presumed that the 
added means of displacement at node 9 caused such a 
difference. To ensure that a true optimum had been reached, 
attempts were made to optimize this structure for several 
different initial starting points. Consistent results were not 


obtained. 


Loads 


5,000 lbs 
9, 000 tbs 
1, OO0S iis 


Lateral 
Axial 
4 Moment 


End conditions 


Node 1 O MOD 
Node 9 1 MOD 
Dimensions 
Radius = 32 in 
Theta = 180 degrees 


Total volume 


Volume = 121.28 22. 
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Y. CASE #15: FIXED-HINGED ROLLER ARCH WITH MULTIPLE LOADING 

For the last case studied, a concentrated lateral load and 
bending moment are applied at the midpoint in combination with 
a load acting radially outward distributed along the length of 
the arch. The cross sectional area behaved as anticipated from 
Case #11, inversely proportional to the axial stress. Again, 
to ensure that a true optimum had been reached, attempts were 
made to optimize this structure for several different initial 


starting points. Consistent results were not obtained. 


Loads 
Lateral = 16,000 lbs 
Axial = 0 lbs 
y Moment = 1,000 in-lbs 
Dast rub: = 100 > Tes/in- 
End conditions 
Node 1 O MOD 
Node 9 2MOD 
Dimensions 

Radius = 3z in 
Theta = 180 degrees 

Total volume 
Volume = 285.75 in? 
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Z. CASE #15A: FIXED-HINGED ROLLER ARCH WITH MULTIPLE LOADING 

For comparison, the structure of Case #15 was subject to 
the same bending moment and distributed load while the lateral 
load was doubled in value. As demonstrated by Case #lla, an 
overall increase in the axial stresses results does not effect 
the shape of the dimension and stress curves. However, the 
volume from Case #15 is increased by 79.65%. In comparison, 
doubling the lateral load for Case #11 resulted in an increase 
in volume of 57.95%. Of interest, it appears that thesmanogmum 
of the volume increase is centered around the midpoint were 


the increased load was applied. 


Loads 
Lateral = 32, 000 ibe 
W Axial = O lbs 
Moment = 1,000 in-lbs 
Disoriie == 190 lbs/in 


Ene Conditions 


Noae 1 O MOD 
Noae 9 2 MOD 


Dimensions 


Radius = a2 el 
Theta = 180 degrees 


Total volume 


Volume = 516.58 in? 
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VII. CONCLUSIONS 


The conclusions of this study are as follows: 


The bar-beam model for stress analysis yielded results 
which deviated from known analytical solutions with an 
error of less than 2%. Therefore, the technique of 
modeling arch structures with bar-beam elements is deemed 
a viable approximation. (Chapter V) 


From the specific cases studied, the Sequential Linear 
Programming method, (Method 2), best performed the 
optimization for cantilever arch structures. The Modified 
Method of Feasible Directions (Method 1) best performed 
the optimization for arch structures with reseriveea 
boundary conditions at both endpoints. (Case #2, 2a, 7) 


Reoptimization of an optimal solution has the effect of 
Smoothing the results and reducing the volume of the 
Structure. The effect of this two-stage optimization 
Strategy was more significant for Method 2 than Method 1. 
(Case #3a, 7a) 


The DOT auto scaling function inhibited the optimizer 
performance. (Case #6a) 


Applying moments that produce stresses that oppose the 
stresses produced by a concentrated load reduce the total 
Structure volume required to withstand the combined load. 
Through prestressing one-way loaded structures, more 
efficient structures can be achieved. (Case #4, 4a, 6) 


The cross sectional shape is dependant on the type of 
stress experienced. When bending stresses dominate, the 
optimal cross section forms a tall rectangle limited only 
by the geometric constraint. When axial stresses dominate, 
the optimal cross section dimensions form a Square. (Case 
#5, dy 150) 


Structures which are more statically indeterminate are 
more efficient under identical loading than less redundant 
Structures. (Case #8, 10) 


Asymmetric structures are more likely to produce erroneous 


results due to the limit of the number of elements used to 
Obtain results. (Case #9) 
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e The boundary conditions act as an excitation which follow 
Mme St. Venant principle. The information from the 
boundary condition diminishes such that for a semicircular 
arch, the boundary conditions do not effect the cross 
sectional shape past the arch midpoint. (Case #11, 14) 


Meee zatilon is a function of the initial design starting 
point. (Case #13a) 

From this investigation, the following is Suggested as a 
possibility for future research in the realm of weight optimum 
arch structures: 

e Continue to record results for a comprehensive study of 
all combinations of parameters, loadings, and end condi- 


Eons 


Meeeeimize the arch structure using varied cross sections 
emerm as a C, L, .or I beam, a box beam, or a circular beam. 


¢e Remove the assumption that the arch maintains a constant 
radius of curvature and optimize the arch shape. 


meeeply additional constraints such.as global buckling in 
Order to present amore accurate model. 


le 


APPENDIX A 
DOT PROGRAM PARAMETERS 


The information in the following tables is taken from 


(Reise, | 


SCALAR PARAMETERS STORED IN RPRM 


RPRM (1) Port ey a aoa 


RERMaG2 } Cia OF00s 


Sue 
ei) -0c001=a5ste 
Sal or 
RPRM (5) 

Sat eae 


RPRM ( OF 1 


RPRM (8) 2 a 0.24, AKIX(1)] 


0.0003 
Exner wwe __fe.4____ 


| RPRM(12) DABSTR POmIOON = Aion O eOOCrim 


RPRM (14) -RPRM(20) RESERVED FOR INTERNAL USE 





NOTE: FO = The value of the objective function at the start of 
optimization (for the initial values of X). 


ele 


DEFINITIONS OF PARAMETERS CONTAINED 
IN THE RPRM ARRAY 


DEFINITION 
oll A constraint iS active if its numerical value is 
more poSitive than CT. CT is a small negative number 
2 CTMIN A constraint is violated if its numerical value is 
more positive than CTMIN 
DABOBJ 
ITRMOP consecutive iterations to indicate 
convergence in optimization 


i Ooms 
i 


apace me ~ 7 — ee - 
a s . e 4 5 
6 ' 
° 
s ¢ 8 4 * 


Se 





Maximum absolute change in the objective between 






Relative change in the objective function attempted 


DELOBJ Maximum relative change in the objective between 
ITRMOP consecutive iterations to indicate 
convergence in optimization 
on the first optimization iteration. Used to 
estimate initial move in the one-dimensional search. 


5 DOBJ1 
Updated as the optimization progresses. 
DOBJ2 Absolute change in the objective function attempted 
on the first optimization iteration 
7 DX1 
to estimate the initial move in the one-dimensional 
search. Updated as the optimization progresses 
DX2 Maximum absolute change in a design variable 
attempted on the first optimization iteration. Used 
to estimate the initial move inthe one-dimensional 
search. Updated as the optimization progresses. 
DCH Relative finite difference step when calculating 
gradients 
LO FDCHM Minimum absolute value of the finite difference step 
when calculating gradients. This prevents too small 
a step when X(1) is near zero | 


ah RMVLMZ Maximum relative change in design variable during 
the first approximate subproblem in the Sequential 
Linear Programming Method. This is, each design 
variable is initially allowed to change by +40%. 
This move limit is reduced as the optimization 
progresses. 

2 DABSTR Maximum absolute change in the objective between 
itrmst consecutive iterations of the Sequential 
Linear Programming method to indicate convergence to 
the optimum 

Blas 


DELSTR Maximum relative change in the objective between | 
ITRMST consecutive iterations of the Sequental 
Linear Programming method to indicate convergence to 
the optimum 












Maximum relative change in a design variable 
attempted on the first optimization iteration. Used 













ales 


PARAMETERS IN THE IPRM ARRAY 


WitOG Ag LON NAME DEFAULT VALUE 


PNAME 

IPRM (1) : 
IPRM (2) NDV 
TPRM (3) 

4) 2 
IPRM (5) 
iP RM GG) NCON+NDV, but at least 2*NDV and 

not gmonc Lon DO Nia, 


ee RMI IGMAX 


( 
PERM 


) 
saat zt 
cane ; 
[terRM(@0) /JeRINT@ =| Oe | eae eo) 
PrpRM(11) | IPRNT?, {Oe 
|tprM(i2) _[/iprnt2_ 0) ee 
LiprM(i3__.| owns) |ou 0 ek J 


T PRM 4) TER ie) REO Giga EOIN” al Ua oes 


) 
IPRM (19) NEWITR INTERNALLY DEFINED 
IPRM (20) INTERNALLY DEFINED 
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DEFINITIONS OF PARAMETERS CONTAINED IN THE IPRM ARRAY 







™TAr 

dod \ud Nr « 68 

1 IGRAD Specifies wnetner the gradients are calculated by DOT 
(IGRAD=0) or by the user (IGRAD=1) 

2 ISCAL Design variables are rescaled every ISCAL iterations. Set 
oer ero. ern or. scaling 

2 1MAX Maximum number of iterations allowed at the optimize level 
. The number of consecutive iterations for which the absolute 






or relative convergence criteria must be met to indicate 
convergence at the optimizer level 














Number of columns in constr 












If IGMAX=C, only gradients of active and violated 
GOnsesaints are Calucaved, Jf IGMAX>0, up to NCOLA 

Gradienzs are calculated; including active, violated, and 
near active constraints 














Maximum number of iterations allowed for the Sequential 
Linear Programming method. This is the number of linearized 
Suppo ems SOlVedr 








Impex of consecutive iterations for which the absolute 
ive convergence eee Mus. be mec to indicate 
in the Sequential Linear Programming method 


) 

















Sequential Linear Programming subproblem print. 
JPRINT>O, IPRINT is turned on during approximate linear 
SUPpLOD IEMs his isetor debugging only 


~ 










print scaling factors for the X vector 


IPRNT2=1, print miscellaneous search information. If 
turn on print during one-dimensional search 
Tiss 2s for debuccing only 





PVleenurmber LOGWHGemiteration history information tomethis 
1S USG€r:uUL fOr USIMG pOStprocessing program to plot the 
iteration process. This is only used if JWRITE>O 
















NEWITR No laren: =—-1|, Set “nr ac Ene Stare Of a New iteration, where 
pers MPEneCmmuUroe. Of sohewiLeravilon~j just completed, If 
METHOD<O, 1, this is after each one-dimensional seaarch. If 
f|THOD=2, this is after each approximate optimization If 
JWRITE>0O, the optimization information will have just been 
Welkeen cOmtmies fClule. It YOUrwleh to stop alter each 

ioe nae LOnmOr llLer dupaltlcllar iteration) and then re- 
Start later, NEWITR is a flag to do this. NEWITR is defined 


incernally by Doe 




















Tne number of constraint gradients needed. Ii the user 
Supelies oracient > Tombol, sthiscewilileoe needed. The 
GOnstrvaint numbersmremaywaich Gradients are needed are 
contained in positiooon 1-NGT of the IWK array. NGT is 
defined internally by DOT 
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APPENDIX B 
JUSTIFICATION FOR OMITTING SHEAR STRESSES 


(The following Appendix is taken from [Ref. 9]) 


The shear stress distribution through a beam of 
rectangular cross-section Nas a parabolic diStributioneauen. 
the height of the member. The maximum shear stress, located at 
the neutral axis of the beam, is 

Go = aloes (B.1) 


where T is the maximum shear stress, V is the shear force, 


and A is the cross-sectional area of the beam. [Ref. 8] 
The normal stress due to bending is given by the equation 
O62 4=oMe/ a (B.2) 
where 6, is the maximum normal stress, M is the bending 
moment, and I is the cross-sectional moment of inertia which 
for this case is bh’/12 where b and h are the width and height 
respectively of the cross-section. 
Redefining the normal stress in terms of the cross- 
sectional dimensions yields 
6, = M(h/2) / (bh?/12) 
or 
O, = 6M/nhA (B.3) 
The ratio of the maximum shear stress to the normal stress 
due to bending, is denoted by r and given by the expression: 


L = Trax/Op (ee 


IG 


Meeereuting Equations (B.1) and (B.3) into Equation (B.4) 
yields 


pee ()] ov / A) yom / nA) 


ee — ny (B.5) 
For the cases investigated in this study, the maximum value r 
G@weeetetain 1s when the loading is that of a uniformly 
Geer oured load, p,. Then, where: 
V = pL (B.6) 
M = p,L?/2 (B.7) 
Wee upOn substitution into Equation (B.8) yields 
io eeomin) noe lie 2) 
which simplifies to 
r = h/2L (B.8) 
The use of the beam equation requires the length of the 
beam to be at a minimum ten times the height, that is: 
fe On (B.9) 
TO maximize the value of r, let L equal 10h, the minimum 
allowable length. Substitutiing this value of L into Equation 
(B.8) yields 
Taye Ol) 
Or simply 
sag 8 (B.10) 
Hence, the maximum shear stress accounts for less than 5% of 
the bending stress developed in the structure. Five percent is 


high considering this analysis over-assumed the value of the 


Pi? 


shear stress by assigning the maximum shear stress to the 
entire cross-section of the beam. Moreover, at the outermost 
fibers where 6, iS a maximum, the shear stress is zero. 
Therefore, under the circumstances of this study, the addition 


of shear stresses was deemed to be unwarranted. 


ies 
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APPENDIX C 
ARCH OPTIMIZATION COMPUTER CODE 


ARCH OPTIMIZATION 


TREES ESESESERSSESESLSSSESO ES SO SESE RECESSES ESLLOSSSSLSLCLCCCCCTCCCCLCCCLCCLLLL 


® 


* 


* 


ARCH OPTIMIZATION ANALYSIS CODE * 


® 


RRA AAAARAHARKEHARARAAAKHAARHARHRAMARAHRAAAAAAAARARRARARARARAAHARARARKHARKRKRAARKRAKKAH 


* 


ALPHA....TRANSFORMATION ANGLE OF ELEMENT (ANGLE TO X-AXIS) 
ANGLE....TOTAL ANGLE OF ARCH (IN DEGREES) 
BAVE.....THE AVERAGE BASE DIMENSION ACROSS AN ELEMENT 
BASE.....DOT ARRAY CONTAINING THE ELEMENTAL BASE DIMENSIONS 
BASEL....DOT ARRAY CONTAINING THE ELEMENTAL BASE DIMENSIONS LOWER 
SIDE CONSTRAINT 
BASEU....DOT ARRAY CONTAINING THE ELEMENTAL BASE DIMENSIONS UPPER 
SIDE CONSTRAINT 
BETA ....TRANSFORMATION ANGLE OF ELEMENT (ANGLE TO Y-AXIS) 
B1......BOUNDARY TERMS APPLIED AT END "1" 
B 2......-BOUNDARY TERMS APPLIED AT END "2" 
cI,..,C5.CONSTANTS RELATED TO ELEMENT STIFFNESS COEFFICIENTS 
CLAN.....CONCENTRATED LOAD APPLICATION NODE (THE NODE FX,FY,FM ARE 
APPLIED 
COUNT....COUNTS THE NUMBER OF ITERATIONS COMPLETED 2 
DOF......DEGREE OF FREEDOMS (UNKNOWN DISPLACEMENTS & SLOPES) _ 
DSN......DESIGN VARIABLE FOR EACH ELEMENT 
DESIGN...DOT ARRAY CONTAINING THE ELEMENTAL BASE AND HEIGHT DIMENSIONS 
DESIGNL..DOT ARRAY CONTAINING THE ELEMENTAL BASE AND HEIGHT DIMENSIONS 
LOWER SIDE CONTRAINT 
DESIGNU..DOT ARRAY CONTAINING THE ELEMENTAL BASE AND HEIGHT DIMENSIONS 
UPPER SIDE CONSTRAINT 
DV1BG....DESIGN VARIABLE #1 (BASE DIMENSION) INITIAL ESTIMATE 
DV1LO....DESIGN VARIABLE #1 (BASE DIMENSION) LOWER SIDE CONSTRAINT 
DV1UP....DESIGN VARIABLE #1 (BASE DIMENSION) UPPER SIDE CONSTRAINT 
DV2BG....DESIGN VARIABLE #2 (HEIGHT DIMENSION) INITIAL ESTIMATE 
DV2LO....DESIGN VARIABLE #2 (HEIGHT DIMENSION) LOWER SIDE CONSTRAINT 
DV2UP....DESIGN VARIABLE #2 (HEIGHT DIMENSION) UPPER SIDE CONSTRAINT 
ERK... sees 6X6 ELEMENT STIFFNESS MATRIX IN LOCAL X,Y COORDINATES 
EKPR.....6X6 ELEMENT STIFFNESS MATRIX IN ELEMENT LOCAL COORDINATES 
ELEN..... LENGTH OF ELEMENT 
eee. PT ORCE VECTOR. OF SYSTEM 
FA.......-CONSTANT DISTRIBUTED LOAD OUTWARD FROM END TO END 
FM......-CONCENTRATED MOMENT AT FREE END 
FX.......CONCENTRATED LOAD IN X DIRECTION AT FREE END 
FY......e.CONCENTRATED LOAD IN Y DIRECTION AT FREE END 
G........THE ARRAY OF CONSTRAINT FUNCTIONS 
GAMMA....6X6 ELEMENT TRANSFORMATION MATRIX 
GK.......-(NDOF)X(NDOF) GLOBAL STIFFNESS MATRIX 
HAVE.....THE AVERAGE HEIGHT DIM. ACROSS THE ELEMENT 
HGT......DOT ARRAY CONTAINING THE ELEMENTAL HEIGHT DIMENSIONS 
HGTL.....DOT ARRAY CONTAINING THE ELEMENTAL HEIGHT DIMENSIONS 
LOWER SIDE CONSTRAINT 
BGTU..... DOT ARRAY CONATINING THE ELEMENTAL HEIGHT DIMENSIONS 
UPPER SIDE CONSTAINT 
INDSN....INITIAL (UNIFORM) DESIGN DIMENSION 
INFO.....DOT PARAMETER USED TO SIGNAL THAT THE OPT IS COMPLETE 
IPRINT...DOT PARAMETER USED SELECT THE DATA OUTPUT FORMAT 


IPRM.....-DOT SELECTABLE INTEGER PARAMETERS 
ITERATE. .THE NUMBER OF TIMES DOT IS TO BE RELOADED WITH THE 


PRECEEDING DATA 


.DOT INTERNAL WORK SPACE ARRAY 


ee 


METHOD...DOT PARAMETER USED TO DEFINE THE OPTIMIZATION METHOD 
MINMAX...DOT PARAMETER USED TO MINIMIZE/MAXIMIZE THE PROBLEM 
NCON.....NUMBER OF DESIGN CONSTRAINTS 
NDO Po. NUMBER OF DEGREES OF FREEDOM 
NDVeyas « =. NUMBER OF DESIGN VARIABLES 
NEL......NUMBER OF ELEMENTS 
NRIWK....DOT INTERNAL WORK SPACE ARRAY DIMENSION 
NRWK..... DOT INTERNAL WORK SPACE ARRAY DIMENSION 
NSNP saree NUMBER OF SYSTEM NODAL POINTS 
OBI. eer THE OBJECTIVE FUNCTION OF THE OPTIMIZATION 
OPTDCS...OPTIMIZATION DECISION TO OPTIMIZE THE PROBLEM OR NOT 
Pl1...P5..PARAMETER DIMENSION CORRESPONDING TO THE NEL, NSNP, NCON, 
NDOF, AND NDV RESPECTIVELY 
PH eee SUBTENDED ELELENT ANGLE (ALSO, PHIANG IN DEGREES) 
PRCSN....THE PRECISION DESIRED TO SOLVE THE FEM SYSTEM OF EQUATIONS 
RADIUS...ARCH RADIUS 
.....DOT SELECTABLE REAL PARAMETERS 
SIGMA_B..THE ELEMENTAL NORMAL STRESS DUE TO BENDING 
SIGMA_N..THE ELEMENTAL NORMAL STRESS DUE TO AXIAL FORCES 
SIGMA_T..THE MAXIMUM TOTAL STRESS IN EACH ELEMENT 
Sewanee THE “DISPLACEMENT” VECTOR OF THE SYSTEM OF LINEAR EQUATIONS 
WK.......DOT INTERNAL WORK AREA 
X........GLOBAL HORIZONTAL COORDINATE 
b ororon -GLOBAL VERTICAL COORDINATE 
YIELD....YIELD STRENGTH OF THE ARCH MATERIAL 
YOUNG....YOUNG’S MODULUS OF THE ARCH MATERIAL 
ReARARM RRA aRRRAERRAERKARKKKKRKAKKEKRAHKEKRKKRKKKKAKRARKAKKAHRHRKKRKKRERRERKRRRERKAKKKKRa HK 
-..-declare the. variables 32... ee 0 6 5 ee ei eeiels ois selene nn 
INCLUDE ’ARCH_COM. FOR’ 


_ 
- 


ona orrrnrneneee te + Hee HH OO UU UU UU UU 
ys] 
Y 
yt] 
4 


cee read the input parameters eoeoeeeeee#eeeeee#eee oeoeeeeosqo5oxefe#eeoeeeeoeoeee#eeeseeoesses8 @ 
OPEN(8, FILE=’ARCH_IN.DAT’, STATUS=’OLD’ ) 


READ(8,*) ANGLE, RADIUS, YOUNG, YIELD,NEL,METHOD,IPRINT,DVI1BG, 


& DV1LO,DV1UP,DV2BG,DV2L0,DV2UP,CLAN,FX,FY,FM,FA,OPTDCS, 
& ITERATE, PRCSN,BX1,BY1,8M1,BX2,BY2,BM2, LABEL 
e 
e -.. -Gefine COnstanmtes. 2h. eRe < 0 ec eisseee AP mer ore 0 « So ee te eum 
NSNP = NEL + 1 
NDOF = 3*NSNP 
NCON = 3*NSNP 
NDV = 2*NSNP 
C 
C ...-.-determine the system nodal coord and element orientation.... 
CALL GEOMETRY (NEL,NSNP,ANGLE,RADIUS,X,Y,ALPHA, BETA, ELEN) 
©: 
G -.. Gefine thes size of the work acrays for e9bOT.°3. (>)... eee 
NRWK = 38800 
NRIWK = 1000 
C 
C -.--Optimize! the problems... eS eas See ee eee Pe 
CALL OPTIMIZATION TOOL 
C 
(e ----compile and format the GUCpie. «ce cr meeet etter eee - 
CALL ARCH OUTPUT 
C 
END 
RecA RK MRAM RKAKRR RAK KRAMER ARHRRKKKNRKKERKKKKARKEKKKKAKKHRKASK 
* * 


SUBROUTINE GEOMETRY (NEL,NSNP,ANGLE, RADIUS,X,Y,ALPHA, BETA, ELEN) 
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| This routine is used by main ARCH OPTIMIZATION to generate | 
| the x-, y-coordinates of each system node, to determine | 
| the orientation of each element, and to calculate the | 
| length of each element. | 
meseaeclace the variables... ..csccsvccssecs Meer ee Tonish sisi) cateneais 
INTEGER NEL,NSNP,P1,P2 


AANANAANAN 


PARAMETER (P1=32,P2=33) 


REAL ANGLE,RADIUS,ELEN,X(P2),Y(P2),ALPHA(P1),BETA(P1), 
& PI,PHI,ANG, YNUM,XDEN 


PARAMETER (PI#3.141593) 


C merecerermine the geometric constants......csssccesee i cee 
PHI = (ANGLE/NEL)*(PI/180.0) 


X(1) = 0.0 
men) = 0.0 


ANG = 0.0 


DO 100 isl, NEL 
ANG = ANG + PHI 
X(i+1) = RADIUS * (1.0 = COS(ANG) ) 
Y(i+l1) RADIUS * SIN(ANG) 
YNUM eed le) eect 1) 
XDEN CX i+ l= eG) ) 
ALPHA( i) ATAN 2 ( YNUM, XDEN ) 
BETA(1i) (PI/2.0) - ALPHA(i) 
100 CONTINUE 


c meraetermine the length of each element... ......ccssereccccee 
ELEN = SQRT(X(2)**2.0 + Y¥(2)**2.0) 


RETURN 
END 


RHR BKM RRRRRAKRRRKRKRRARRKEKARRRRRARRAAARRAEERARARKRKRKRKKHKEKRARAKRARKRAAKAE 
* * 
SUBROUTINE OPTIMIZATION TOOL 

This subroutine directs the program flow optimization decision 
i.e., optimize the problem or not. It also serves to set up & 
execute the DOT optimization software. 

~--- Geclare the eT JEAN) ee Gieecs. che cree er teen. coesyeiets «cits ceeéececeeee 
INCLUDE ’ARCH COM.FOR’ 


AANAAN 


INTEGER i 


aa 


fe were out the RPRM and IPRA @t rays. 2... wee eee eee 
DO 100 i#1,20 
RPRM(i) = 0.0 
IPRM(i) = 0 
100 CONTINUE 


c emo ialine “COUNT... ... re ree oO 6 OMY Se S665 5 Ce 
COUNT = 1 


ya 


aqa-n 


gaa AN 


an 


qgaaQq Aan 


200 


250 


300 


ec s Define the Constraint ce uacrence se |e cee « OR 
RPRM(2) = 0.0001 
RPRM(3) = 0.001 


oe turn Off DOT s auto scalingeeees. or. os oo ee eee «(0 5 Om nee 
IPRM(2) = -1l 


wee increase DOT’s default number of iterations........... eee 
IPRM(3) = 1000 
IPRM(8) = 1000 


...-increase DOT’s number of consecutive convergence criteria. 
IPRM(4) = 3 
IPRM(9) = 3 


»...define MINMAX=-1 to minimize the objective function....... 
MINMAX = -1] 


...-imitialize the design variable limits and best guess...... 
DO 200 i=#=1,NSNP 


BASE(i) = DV1BG 
BASEL(i) = DV1LO 4 
BASEU(i) = DV1UP . 5 
HGT(i) = DV2BG 
HGTL(i) = DV2LO 
HGTU(i) = DV2UP 
CONTINUE 


....combine base and HGT arrays into design array......... 
DO 250 i#1,NSNP 
jJ=NSNP+1 
DESIGN(i) = BASE(i) 


DESIGNL(i) = BASEL(1) 
DESIGNU(1) # BASEU(1) 
DESIGN(j) = HGT(i) 
DESIGNL(j}) = HGTL(i) 
DESIGNU(j) = HGTU(i) 
CONTINUE 
~..make optimization decision............. 2 see. eee + ee 


IF (OPTDCS .NE. 1) THEN 
CALL EVAL 
RETURN 

ENDIF 


.... ready to OptiMiZe@.. .. «css 6 cas & 0 6 ese eenememenenene @eeee#?@@ @ee47e#eese?¢ @ 
INFO = 0 


CALL DOT (INFO,METHOD, IPRINT,NDV,NCON, DESIGN, DESIGNL, DESIGNU, 
OBJ ,MINMAX,G,RPRM, IPRM,WK,NRWK, IWK, NRIWK) 


...-evaluate the objective function and constraints........... 
IF (INFO .GT. 0) THEN 

CALL EVAL 

GOTO 300 
ENDIF 


AZ 


mee refinemstiemsolucion Vector bY GE0ptimi zing. sccsasuseoes-.. 
fe) (COUNT .LYT. ITERATE) THEN 
INFO = Q 
COUNT = COUNT+1 
GOTO 300 
ENDIF 


RETURN 
END 


RR RRR RARE ERE RR ERRARRKRERR RRR R RRR RERRRERRRRRREREAKARE 


® 


ANANAANANN 


an 


50 


100 


200 


210 


® 
SUBROUTINE EVAL 
, oe se & £ 8 8 6 8 EF EE EE EE BR EG UE EE UB EG UG UE UE UE US U8 8 8 UE UE UB UE UU UE UU UE UF UE EF EG gg eg FF ee EF Fe ee eg oe og oe OP ee se Ue Cy 
This subroutine is used to evaluate the Objective function, 
constraint functions, and side constraints of the optimization 


problem. 
Pr CI@eLAre Che. Var  ableS <ccc6 4.4 6 oe cies Ss. ee os Bs 5 5 


INCLUDE ’ARCH COM.FOR’ 
INTEGER i, 3 


....separate the design array into base and HGT arrays 
DO 50 i#=1,NSNP 
J=eNSNP+i 
BASE(i) = DESIGN(i) 
HGT(i) = DESIGN(j) 
CONTINUE 


...-- calculate the objective pa ghee Vol lay We Ug ee 
OBJ = 0.0 


DO 100 i#1,NEL 
BAVE(i) = (BASE(i)+BASE(i+1))/2.0 
HAVE(i) = (HGT(i)+HGT(i+1))/2.0 
OBJ = OBJ + BAVE(1)*HAVE(i)*ELEN 
CONTINUE 


mm lierarice che Ge@sign Constraint: VeECtOr. 20... etc eee cee es 
DO 200 i#=1,NCON 

G(i) = 0.0 
CONTINUE 


...-Getermine the design constraints..... Oe eS ee ecinoete 
CALL ARCH_STRESS 


DO 210 i=#1,NSNP 
j= i+NSNP 
kei+(2*NSNP) 
G(i) = (SIGMA T(i)/YIELD) - 1.0 
G(j) = BASE(i)-(3.0*HGT(i)) 
G(k) = HGT(i)-(10.0*BASE(i) ) 
CONTINUE 


RETURN 
END 


eM MR RRR RRR RRA EARAHREREAKAREE AKER ERRER EKER 


® 


* 


SUBROUTINE ARCH STRESS 


F235 


ANANAN 


‘@) fo) 2) NAA AANA AN AN 


This subroutine is used to perform the Finite Element analysis 
of the stresses developed in an arch or beam for a given load- 
ing. 

...-declare the variables........ eee re re er er ee i 
INCLUDE ’ARCH_COM.FOR’ 


INTEGER IPVT(99) 
REAL F(P4) 


REAL*8 BK(P4,P4),BF(P4),BU(P4),FAC(9801),WORK(99) 


....form the element and system matrices........... i: eager ere ene 
CALL FORM 
form the FOrce Vector, ©... ee éSevedeeenene 


CALL FORCE VECTOR (NEL,NDOF,ELEN, ALPHA, BETA, FA,F) 


.set the boundary conditiona and loads. ne 
CALL BNDARY (NDOF,GK,CLAN,FX,FY,FM,F,BX1, BY1, /BM1, /BX2, BYZ .BM2) 


... BOlVve Ehe System Of s€équations. .. 3. . sneme EP eer oe 
IF (PRCSN .EQ. 2) THEN 
....-Cchange GK and F arrays to double precision...... ee 
CALL UPSCALE (NDOF,GK,F,BK,BF) 
oe SOLVE the system oOfequatlonic. . crs. 0 neueeeees ya 
CALL DL2ARG (NDOF,BK,P4,BF,1,BU, FAC, IPVT, WORK ) 
.... change BU arrayetousingle presi cicn au. oe 
CALL DOWNSCALE (NDOF,BU,U) 
ELSE 
.-solve the system of equations. ROM Gok ec 
CALL L2ZARG (NDOF,GK,P4,F,1,U,FAC, IPVT, WORK ) 
ENDIF 


~oe-Getermine the stress disStribution.......... 2 ig 
CALL STRESS 


RETURN 
END 


REAEKKKRRERRHRREKEKEKEKEKAREAREKREKKRRRKEKRRKAEKRHKRAEKRKEKARKRKKKKKEKAEKRERKKERKERKEKERE 


® 


AANANA 


® 


SUBROUTINE FORM 


This subroutine is used to construct the global stiffness mat- 
rix for the arch problen. 


....aeclare the variables........ oeee#eeee oee#eskeeeeee+eet3eeesteeee#e#eee#e 
INCLUDE ‘ARCH COM. FOR’ 


INTEGER IEL,1,J,%,11,JJ,K, [11,009 
REAL Ci,c2,C3,C4,C5,CA,CB,Ek(P1,6,6) ,GAMMA(6,6), ENXGA( 6,678 
GAEKGA(6,6),BH,BH3 


os2e0Ge@fine the constants ¢x....-- sie: 6. 6 eon. Cee ee er i a 
C1 = YOUNG/ELEN 


124 


ZnO 
205 
220 


Za0 
Z5p 
240 


Z50 
260 
SAG 


140 
150 


280 
290 
300 


Boe220 I = 1,6 
DO 215-8 
DO 210° K 
EKGA ( 
CONTINUE 
CONTINUE 
CONTINUE 


..determine the 
DO 240 I = 1,6 
DO. 235 J = 
DO 230 K 
GAEKG 
CONTINUE 
CONTINUE 
CONTINUE 


....-copy the GAEK 
DO 260 I = 1,6 
DO 250 J = 
ER( LEG 
CONTINUE 
CONTINUE 
CONTINUE 


me. 2nicialize th 
DO 150 I = 1, NDO 
po 140 J #1, 
GK(I,J) = 0 
CONTINUE 
CONTINUE 


mee.constrcuct the 
DO 300 IEL = 1, N 
II = 3*(I 
BO.290 J = 1, 
JJ = II 
DO 280 K 
KK 
GK(JJ 
CONTINUE 
CONTINUE 
CONTINUE 


RETURN 
END 


as 
= 1,6 
I,J) = EKGA(I,J) + EKPR(IEL,1I,K) *GAMMA(K,J) 


GAEKGA array.......- DE oo 6 eo 1 6 6 6 6a as eee 


reo 
= 1,6 
A(I,J) = GAEKGA(I,J)+GAMMA(K,I)*EKGA(K,J) 


GA array into the EK array........ 5. ee 


6 
,J) = GAEKGA(I,J) 


See Corte LO Vos: «ements memento. 5's i slemeasl ste tetene igtetels e+ « « 
F 

NDOF 

20 


GR miele x... ere, 6 + oe edice ne 6% AE ee 
EL 

EL-1) 

6 
+ J 

= 1, 6 

= II + K 

»KK) = GR(JJ,KK)+EK(IEL,J,#) 


Raentkkkkkknrkhkakekkakkkkkakkkkkkkakkakkkakkakkkaekkekakkkakkenkke ek kkkkeakaekahkakank 


® 


qa A ANNAN 


* 


SUBROUTINE FORCE VECTOR (NEL,NDOF,ELEN,ALPHA, BETA, FA,F) 


This subroutine i 
FEM problem speci 


PARAMETER( P1=32,P 


s used to construct the force vector for the 
fied. 


4=99 ) 


REAL ELEN,ALPHA(P1),BETA(P1),FA,F(P4) 


Ze 


90 
100 


C2 = (1.0/ELEN) ** 2ya0 

C3. = (12.0)7(2 0 “EEN? 

C4 = (1.07/80 

€5 = (1707) c.0 

ose eb Nit PalizesGner wOUK at bavis en ars 
DO 120 IEE si, NEE 


DO 10071 = 176 
DO 90 J= 1,6 
EKPR(IEL,I 
GAMMA(I,J) 
EKGA(I,J) 
GAEKGA(I,J) 
BK{ TEL, 173) 
CONTINUE 
CONTINUE 


,J) 

= Q. 
= 0. 
= 0. 
= Q. 


...-Calculate the area and 
BH = BAVE(IEL) *HAVE(IEL) 


inertia terms 


BH3 = BAVE(IEL) *(HAVE(IEL) **3.0) 


...-Getermine 
EKPR(IEL,1,1) 
EKPR(IEL,1,4) 
EKPR(IEL,2,2) 
EKPR(IEL,2,3) 
EKPR(IEL,2,5) 
EKPR(IEL,2,6) 
EKPR(IEL,3,2) 
EKPR(IEL,3,3) 
EKPR(IEL,3,5) 
EKPR(IEL,3,6) 
EKPR(IEL, 4,1) 
EKPR(IEL,4,4) 
EKPR(IEL,5,2) 
EKPR(IEL,5,3) 
EKPR(IEL,5,5) 
EKPR(IEL,5,6) 
EKPR(IEL,6,2) 
EKPR(IEL,6,3) 
EKPR(IEL,6,5) 
EKPR(IEL,6,6) 


~e. Getermine 


the ERPR matrixeee....:..<2. 


C1*8H 
-C1*BH 
C1*C2*B8H3 
ei es *8u3 
-C1*C2*BH3 
Cl*C3*BH3 
Giwies <5 as 
C1*C4*BH3 
-C1*C3*BH3 
elves *BH3 
-C1*BH 

ei? BH 
-C1*C2*BH3 
=-Cl*c53. hes 
€1*€2*BH3 
=Cl*C3" Bus 
eles "Sus 
C1*C5*BH3 
-C1*C3*BH3 
C1*C4*BH3 


thie GAMMA Mmatrix....+... 
CA = COS(ALPHA(IEL) ) 
CB = COS(BETA(IEL) ) 


GAMMA(1,1) = CA 
GAMMA(1,2) = CB 
GAMMA(2,1) = -CB 
GAMMA(2,2) = CA 
GAMMA(3,3) = 1.0 
GAMMA(4,4) = CA 
GAMMA(4,5) = CB 
GAMMA(5,4) = -CB 
GAMMA(5,5) = CA 
GAMMA(6,6) = 1.0 
...-determine the EKGA array............ 


126 


oe eoeeeee#eeseee#e#eeeee##%#### @ @¢ @ @ 


oe ef ¢@ @ @ @cGHhUhM OCU OCUCUCOrCUCUcCOrhUCUCOrChMUcChOrhUCUchOhrmhUlUcPOrhUhCUchOhrhUCUchhmhUCUcChhUhhU Fh} 


eoeet?eee8e#e#e@e@6mUMOMmUCUCOrhUCUC HhrhUCUchmhUhHhUhOrCUhLhChUhOhUCUcCOrmhUCUC OrhUCc HhCU HhlUhO 


oeoeeeeeee#srkesee##eeeeeeeeeeesese#s#eese @ 


F(l) = (ELEN/2.0) * (-COS(BETA(1))) 
F(2) = (ELEN/2.0) * (COS(ALPHA(1))) 
F(3) 0.0 


DO 100 i#2,NEL 
Il] =» (i-1)*3 + 1 
I2 = (i-1)*3 + 2 
I3 = (i-1)*3 + 3 


E( 12) ose (ELEN72..0)*({(-COS( BETA(1))) 
& +(ELEN/2.0)*(-COS(BETA(i-1))) 
F(I2) = (ELEN/2.0)*(COS(ALPHA(i))) 
& +(ELEN/2.0)*(COS(ALPHA(i-1)) ) 
F(I3) = 0.0 
100 CONTINUE 


F(NDOF-2) = (ELEN/2.0)*(-COS(BETA(NEL) ) ) 
F(NDOF-1) = (ELEN/2.0)*(COS(ALPHA(NEL) ) ) 
F ( NDOF ) = 0.0 


Cc momescalesthe F-vector by FA......%. “carnrien IEA Oo a 
DO 200 i#1,NDOF - 
F(i) = FA*F(i) = 
200 CONTINUE . 


RETURN 
END 


MRR HRRHHRRR ARR RRRRR eR RRRHRRRRRRE RK eee RRR ARERR KRKAERAE KK 
x ® 
SUBROUTINE BNDARY (NDOF,GK,CLAN,FX,FY,FM,F,BX1,BY1,BM1,BX2, 

& BY2,BM2) 
This subroutine is used to impose the boundary conditions upon 
the global stiffness matrix and force vector. 
meoec lace thes varlaOleg... 6.2565 « «ns MR ee cass bc ce te 
INTEGER NDOF,BX1,BY1,BM1,BX2,B8Y2,BM2,CLAN,i,N,I1,12,13,P4 


AAANANA 


PARAMETER( P4=#99) 
REAL GK(P4,P4),FX,FY,FM,F(P4) 


G invoke, ther wseential DOouN@arY CONGICIONS. ......seesc esse wes 
TP (BXl .EQ. 1) THEN 
CALL IMPOSE BC (NDOF,GK,1,F) 
ENDIF 


fF (eyl’.EG. 1) THEN 
CALL IMPOSE BC (NDOF,GK,2,F) 
ENDIF 


IF (BMl .EQ. 1) THEN 
CALL IMPOSE BC (NDOF,GK,3,F) 
ENDIF 


IF (BX2 .EQ. 1) THEN 

N=NDOF-2 

CALL IMPOSE BC (NDOF,GK,N,F) 
ENDIF 


LZ 


IF (BY2 s5ow 1) THEN 
N=NDOF-1 
CALL IMPOSE_BC (NDOF,GK,N,F) 
ENDIF 


IF (BM2 .EQ. 1) THEN 
CALL IMPOSE BC (NDOF,GK,NDOF,F) 


ENDIF 


G ....add the concentrated load to the force vector............ 
Il=#(CLAN-1)*3+41 
I2=(CLAN-1) *3+2 
I3=(CLAN-1) *3+3 


F(I1l)#F(I1)+FX 
F(I2)eF(I2)+FY 
F(I3)=2F(I3)+FM 


RETURN 
END 


See eee EERE EER RERR RRR RRR RRR RRR RRR RRR RRR RRR RR RRR RR RRR RRR RRR ES ESS ES 
® zi * 
SUBROUTINE IMPOSE BC (NDOF,GK,N,F) = 

This subroutine is used to do the redundant leg work of impos- 
ing the boundary conditions. 

-.. Geclare the variables... vc ctesececececectvseceeeec ee 
INTEGER NDOF,N,1,P4 


AaAANAANAN 


PARAMETER(P4#=#99) 
REAL GK(P4,P4),F(P4) 


E ...e-impose the boundary condition on ate GK and F arrays...... 
DO 100 i#1,NDOF 
GK(N,i) = 0.0 
100 CONTINUE 


RETURN 
END 


kkkkrkkrkaknrkkkkkkkkkkkkkkakkhkkkkhkkkkkkkrkkkkkkakhkkkkkkkhkakkkekaekeaekakaekenenk 
* ® 
SUBROUTINE UPSCALE(NDOF,GK,F,BK,BF) 

This subroutine is used to change the stiffness matrix & force 
vector from single precision to double precision in order to 
solve the linear system of equations in double precision. 
-+--Geclare the variables ...<s.5e - eee on ee errr 5 
INTEGER NDOF,i,j,P4 


AAANANANAN 


PARAMETER (P4=299) 
REAL GK(P4,P4),F(P4) 


REAL*8 BK(P4,P4),BF(P4) 


12s 


»...generate the doubleprecision compliments of GK and F...... 
DO 110 i#1,NDOF 
DO 100 j=#=1,NDOF 
BK(i,j) = GK(i,j) 
CONTINUE 
BF(i) = F(i) 
CONTINUE 


RETURN 
END 


RRMA ERR RRAMARRAAARREMRAAAARREAARHERRAKRRARREARAARARREAAKRKRRKAKRAKKRKRAREE 


® 


ANAANAN 


100 


® 
SUBROUTINE DOWNSCALE(NDOF,BU,U) 
This subroutine is used to do down scale the double precision 
solution of the linear system of equations back to single pre- 
cision. DOT could have problems with double precision numbers! 
merece clace the Variables... .cccrccvcces oi’ «Sere eretemene ehenécé-is wee ee 
INTEGER NDOF,1,P4 


PARAMETER (P4#=#299) u 
REAL U(P4} 
REAL*8 BU(P4) 
...-generate the doubleprecision compliments of GK and F...... 
DO 100 i=1,NDOF 
U(i) = BU(i) 
CONTINUE 


RETURN 
END 


-c eee cece eee ee eee eee eee eee eee eee eee eee eee eee ERR RRR RRR ERR REE Eee ee ee 


* 


aaNAaAN 


® 
SUBROUTINE STRESS 
This subroutine computes the stress at each nodal point. 
Meee OeClarations<.:eec.sccccc es eee # eoee#e?e¢#eee8 @®#8 @ @ o@e@@e@ee¢eee#?8¢#¢e08#¢#¢ ¢#2 ¢ @ 
INCLUDE ‘ARCH _COM.FOR’ 


INTEGER tll 2 374, 15516 
REAL CA1,CB1,K1,K2,FPR(P4,6),UPR(6),NORM1,NORM2, 
BEND1,BEND2 


...-determine local forces from stiffness and displacement.... 


DO 100 i=l, NEL 
Tle(i-1)*341 
T2=(i1-1)*3+2 
I3m(i-1)*3+3 
I4ma(i)*341 
T5=(i)*3+2 
T6=(i)*3+3 


CBl= COS(BETA(i)) 


Ze 


CAl= COS(ALPHA(1) ) 


C 
UPR(1)= U(I1)*CAl + U(I2)*CBl 
UPR(2)= -U(I1)*CBl + U(I2)*CA1 
UPR(3)= U(I3) 
UPR(4)=—= U(I4)*CAl + U(I5)*CBl 
UPR(S)= -U(I4)*CBl + U(1IS)*CAl 
UPR(6)= U(I6) 
C 
DO 250 L#1,6 
FPR(i,L)= 0.0 
250 CONTINUE 
S 
DO 300 Jel,6 
DO 350 K=1,6 
FPR(i,J)= FPR(i,J) + EKPR(i,J,K)*UPR(K) 
350 CONTINUE 
300 CONTINUE 
100 CONTINUE 
G ... determine the bending and normal stresses...........22000. 
SIGMA_N(1) = ABS(FPR(1,1)*(1.0/(BASE(1)*HGT(1)))) 2 
SIGMA _B(1) # ABS(FPR(1,3)*(6.0/(BASE(1)*(HGT(1)**2.0))))- 
SIGMA_T(1) = SIGMA_B(1) + SIGMA_N(1) = 
DO 400 i=#2,NEL 
Kl = 1.0/(BASE(i)*HGT(i) ) 
K2 = 6.0/( BASE(i)*(HGT(i)**2.0) ) 
NORM1 = ABS(FPR(i,1)*K1) 
NORM2 = ABS(FPR(i-1,4)*K1) 
BEND1 = ABS(FPR(i,3)*K2) 
BEND2 = ABS(FPR(i-1,6)*K2) 
SIGMA_N(i) = (NORM1+NORM2)/2.0 
SIGMA_B(1i) = (BEND1+BEND2)/2.0 
SIGMA_T(i) = SIGMA_B(i) + SIGMA_N(i) 
400 CONTINUE 
SIGMA_N(NSNP) = ABS(FPR(NEL,4)*(1.0/(BASE(NSNP) *HGT(NSNP) ) ) ) 
SIGMA _B(NSNP) = ABS(FPR(NEL,6)* 
& (6.0/(BASE(NSNP) *(HGT(NSNP) **2.0)))) 
SIGMA_T(NSNP) = SIGMA_B(NSNP) + SIGMA_N(NSNP) 
C 
RETURN 
END 
RRARARKRKKRRAKRAKRRAARARRAAARAKRAARARRARAKARAAARARKEAKRAARARRAARKRARAAARKAAARRARRAKREAN 
x ® 
SUBROUTINE ARCH OUTPUT 
BREESE RBSSBEPSASSESEESISSSBTESTSVsESTESCECASsBSwSESswlsISSESPPSSISESSESEESISESPSsSASSSEPSBPSESBEA_AASsA 
Cc This subroutine formats the final results and output of the 
Cc optimization problem and stores it in a file named ARCH_OUT.DAT 
Cc BSR SSSSaBBTeBFFSBFPIsI SF SPPSSPSISS SPSS SSFP SISQPISASBSSASASVTSRAPE SAA SPSSISAASBSsIsASA S222 BS 
c .-. declare vaciabl @sg cs. ca cc ccc o's clemenememeneeeeeienens, «eu oe 


INCLUDE ‘ARCH COM. FOR’ 


ISO 


REAL VOL, VOLUME 


Ren OUtpULC El lemand writer header. ..0 2. 6. cc wc cee ee twee 
OPEN(9, FILE=’ARCH OUT.DAT’, STATUS#’ UNKNOWN’ ) 


WRITE(9,100) LABEL 
WRITE(9,100) ’ OPTIMIZATION SOLUTION’ 
WRITE(9,105) 0 ww ww nm wn we we ow eo ne eo ee ee ee ee 


100 FORMAT(/5X,A) 
10S FORMAT (5X,A) 


SDM NCRINCON IRN UIT ES sO ITE es) So choles 0 6 6 6 6-656 6 oe cs cso Glee eae eee es oe 
WRITE(9,100) ’ A) Problem Parameters: ’ 

WRITE(9,110) ’ Arch Angle :’, ANGLE, ’ Youngs Modulus:’,YOUNG 
Weete(9,110) * Arch Radius:’, RADIUS, ’ Yield Strength:’ ,YIELD 


WRITE(9,115) ’ No of Design Var:’, NDV, ’ No of Elements:’,NEL 
NEAL, PORMAT(8X,A,F12.3,7T38,A,F12.1) 
1s FORMAT(8X,A,I7,T38,A,I10) 


RAST GGT RAE Song SESS Goi eee ee aerate: 
WRITE(9,100) ’° B) Derived Constants:’ “ 
WRITE(9,120) ’ No of System Nodal Points...’,NSNP ~ 
WRITE(9,120): ’ No of Degrees of Freedom....’,NDOF 
WRITE(9,125) ’ Length per Element.......... ‘, ELEN 
WRITE(9,125) ’ Phi Angle per Element....... ‘, PHIANG 
Mmeere(9,120) " Number of Iterations........ ‘, ITERATE 

120 FORMAT(8X,A,1I6) 

125 FORMAT(8X,A,F12.4) 
(oe sp SOS UE “CM oa So bh 6 Se. Seer ee en 
WRITE(9,100) ’ C) Structure Loading:’ 
“of ca a(S) aS) aut 13> <a eee (Ex 
LS Sle (CS) Gale) I a aes eS rs ee er ae Cpe yt 
Premmteiee SC EM 5.66 score os 06 os 6 #6 6 8 8 9 eo OpEM 
mero lees) PAG ss okt oh cde ea ewe w eee es aS ee 
ESOC De 5 65s env esis vie sucleisia ses 6 6 6 $84 6 H6 HO See EO oe ee we ee 


WRITE(9,100) ’° D) Elemental Dimensions and Stress Distribution: ’ 
WRITE(9,210) ’Node’,’Height’,’Base’,’Length’,’Area’ 


210 FORMAT(8X,A,T21,A,T36,A,T49,A,T62,A) 
220 Pommar(S8X,14,T17,F10.5,T32,F10.5,T48,F8.5,T60,F8.5) 
VOLUME = 0.0 


DO 300 i#1,NSNP 
AREA = BASE(i)*HGT(i) 
WRITE(9,220) i,HGT(i),BASE(i),ELEN,AREA 
300 CONTINUE 


ee Ss COLEG) OL ig Bie rem aM ce, eee ee ee Go he Voravie tots 660 sie 6 4 0 8 0 eee 
WRITE(9,100) ° E) Obj@ctive Function: ’ 


WRITE(9,310) ’ Total structure Volume:’ ,OBJ 
310 FORMAT(/12X,A,F12.6/) 


WRITE(9,330) ’Node’,’Normal Stress’,’Bending Stress’,’Total’ 
DO 320 i#1,NSNP 


i 


400 
410 
420 
430 


WRITE(9,340) i,SIGMA_N(i),SIGMA_B(i),SIGMA_T(i) 
CONTINUE = 
FORMAT(8X,A,T18,A,T35,A,T57,A) 
FORMAT(8X,14,T15,F14.1,1732,F14.1,T49,F14.1) 


oe ee  SOCCION UF 6. sk ce ee 6s Cee ee. 
WRITE(9,100) ’ F) Boundary Conditions:’ 

WRITE(9,410) ‘Node’,’X-Displ’,’Y-Displ’,’Slope’ 

WRITE(9,430) 1,BX1,BY1,BM1 

WRITE(9,430) NEL+1,BX2,BY2,BM2 


oo e SECKION "GN. ob sik cbse. ses eee 
WRITE(9,100) ’ G) Solution Vector:’ 
WRITE(9,410) ‘'Node’,’X-Displ’,’Y-Displ’,’Slope’ 
DO 400 i=z1,NSNP 

Tl=#(i-1)*3+1 

I2=(i-1)*3+2 

IT3=(i-1)*3+43 

WRITE(9,420) i,U(I1),U(I2),U(13) 

CONTINUE 
FORMAT(T9,A,T17,A,T31,A,T46,A) 
FORMAT(7X,1I5,3E14.6) 
FORMAT(7X,15,T720,14,T734,14,T48,14) 


RETURN 
END 


2 
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2) 


DM M PH Mm Mm M 


MMMM MM M M 


ARCH COMMON 


> abel eT ee ee ee Ce eAet eee Karcen ec) i ieee ae 
Pl. ....The maximum number of elements 
P2.....The maximum number of global nodal points 
P3.....The maximum number of design constraints 
P4.....The maximum number of degrees of freedom 
pS.....The maximum number of design variables 


Nerelemeacewerine. VALLAD LOS. 66. 5 0 etc tt 8c he ee hc ecw ewe a webs 
INTEGER NEL,NCON,NSNP,NDOF,NDV, METHOD, MINMAX, INFO, IPRINT, 
IWK(1000), NRWK,NRIWK,IPRM(20),COUNT, OPTDCS, ITERATE, 
PRCSN, CLAN, BX1,BY1,BM1,BX2,BY2,BM2,P1,P2,P3,P4,P5 


PARAMETER( Pl#32,P2=33,P3=96,P4=#99,P5=64) 


REAL ANGLE, RADIUS, ELEN, X(P2),Y¥(P2),ALPHA(P1),BETA(P1), 
YOUNG, YIELD,WK(38800),RPRM(20),0BJ,G(P3), 
DV1BG,DV1LO,DV1UP,BASE(P1),BASEL(P1),BASEU(P1), 
DV2BG,DV2L0,DV2UP,HGT(P1),HGTL(P1),HGTU(P1), 
DESIGN(P5) ,DESIGNL(P5),DESIGNU(PS), 
FA,FX,FY,FM,U(P4),SIGMA_T(P4),SIGMA_N(P4), 

SIGMA B(P4),BAVE(P1), 
HAVE(P1),GK(P4,P4),EKPR(P4,6,6) = 


CHARACTER®* 30 LABEL 


Pn A We mmemeeCOMMON . .. oc Nees 6 sc 0c os ec ct Sim cc ee wow eee reece eee 
COMMON NEL,NCON,NSNP,NDOF, NDV, METHOD, “MINMAX, INFO, IPRINT,IWK, 
NRWK,NRIWK, I PRM, COUNT, OPTDCS, ITERATE, PRCSN, CLAN, 
BX1,BY1,BM1,BX2,BY2,BM2, 
ANGLE, RADIUS, ELEN,X,Y,ALPHA, BETA, YOUNG, YIELD, 
WK, RPRM,OBJ,G,DV1BG,DV1LO,DV1UP, BASE, BASEL, BASEU, 
DV2BG,DV2LO,DV2UP,HGT,HGTL,HGTU, 
DESIGN,DESIGNL,DESIGNU, 
FA,FX,FY,FM,U,SIGMA T,LABEL,SIGMA _N,SIGMA B, 
BAVE,HAVE,GK,EKPR 


1S3 


APPENDIX D 
VALIDATION CASES 


Validation #1 


OPTIMIZATION SOLUTION 


A) Problem Parameters: 


Arch Angle : O2003 Youngs Modulus: 30000000.0 

Arch Radius: 1000000.000 Yield Strength: 52000.0 

No of Design Var: ie No of Elements: 4 
B) Derived Constants: 

No of System Nodal Points... 5 

No of Degrees of Freedon.... 15 

Length per Element.......... ll.Z500 

Number of Iterations........ af 


C) Structure Loading: 


Pe a aera ae eG Renews Ret S 1000.0000 

FY .. chs se subeelaehateie erent eee 0.0000 

2 PP a ee 0.0000 .. 

oP PP eee rs cc 0.0000 : 

D) Elemental Dimensions and Stress Distribution: 

Node Height Base Length Area 
i 3.00000 1.50000 11.24996 4.50000 
Zz 3.00000 1.50000 11.24996 4.50000 
3 3.00000 1.50000 11.24996 4.50000 
4 3.00000 1.50000 11.24996 4.50000 
5 3.00000 1.50000 11.24996 4.50000 


E) Objective Function: 
Total structure Volume: 202.499207 


Node Normal Stress Bending Stress Total 
0.0 19999.7 19999.7 
2 0.0 14999.7 14999.7 
3 07:0 9999.8 9999.8 
4 0.0 4999.9 4999.9 
= 0.0 Onc 020 
F) Boundary Conditions: 

Node X-Displ Y-Displ Slope 
1 1 a 1 
5 0 0 0 


G) Solution Vector: 
Node X-Displ Y-Displ Slope 


1 0.000000E+00 0.000000E+00 0.000000E+00 
2 0.257807E-01 0.112327E-08 -0.437492E-02 
3 0.937478E-01 0.409056E-08 -0.74998S5E-02 
4 0.189839E+00 0.828721E-08 -0.937481E-02 
5 0.299993E+00 0.130985E-07 -0.999979E-02 
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Validation #2 


OPTIMIZATION SOLUTION 


A) Problem Parameters: 


Arch Angle : 0.003 Youngs Modulus: 30000000.0 
Arch Radius: 1000000.000 Yield Strength: 52000 <0 
No of Design Var: 10 No of Elements: 4 
B) Derived Constants: 
No of System Nodal Points... 5 
No of Degrees of Freedon.... 15 
Length per Element.......... ein 2 5100 
Number of Iterations........ a 
Smeeructure’ Loading: 
oo 4 Arce ere 0.0000 
6 4 er 1000.0000 
lo 6 er 0.0000 
4 6 A ire 0.0000 a 
D) Elemental Dimensions and Stress Distribution: - 
Node Height Base Length Area 
Ht 3.00000 1250000 11.24996 4.50000 
2 3.00000 1.50000 11.24996 4.50000 
3 3.00000 1.50000 11.24996 4.50000 
4 3.00000 2. 50000 11.24996 4.50000 
5 3.00000 1.50000 11.24996 4.50000 
E) Objective Function: 
Total structure Volume: 202.499207 
Node Normal Stress Bending Stress Total 
il 729) Fa 0.0 22242 
2 7A IA Sane 0.0 C22 ae 
3 we. 2 0.0 toed 
4 222.2 0.0 Zel.2 
5 222.2 0.0 Dei 2 
F) Boundary Conditions: 
Node X-Displ Y-Displ Slope 
a 1 i a 
5 0 0 0 
G) Solution Vector: 
Node X-Displ Y-Displ Slope 
1 Q0.000000E+00 0.000000E+00 0.000000E+00 
2 0.112327E-08 0.833330E-04 -0.191234E-09 
3 0.409056E-08 0.166666E-03 -0.327829E-09 
4 0.828721E-08 0.249999E-03 -0.409786E-09 
5 0.130985E-07 0.333332E-03 -0.437105E-09 


SS) 


Validation #3 


OPTIMIZATION SOLUTION 


A) Problem Parameters: 


Arch Angle : 0.003 Youngs Modulus: 30000000.0 
Arch Radius: 1000000.000 Yield Strength: 520000 
No of Design Var: 10 No of Elements: 4 
B) Derived Constants: 
No of System Nodal Points... 5 
No of Degrees of Freedom.... Ls 
Length per Element.......... dl 00 
Number Of Iterations... .. 2... 
C) Structure Loading: 
pr rer teres Oe 0.0000 
PY oo. 50s aieiece suse ee eecs meee 0.0000 
PM 55 ileus ee sae valance ae mememeeencu sera 10000.0000 
Pre ne ee 0.0000 7 
D) Elemental Dimensions and Stress Distribution: : = 
Node Height Base Length Area 
iH 3.00000 1.50000 11.24996 4.50000 
2 3.00000 1250 00:0 11.24996 4.50000 
3 3.00000 1.50000 11.24996 4.50000 
4 3.00000 Teri oie} u) 11.24996 4.50000 
5 3.00000 1.50000 11.24996 4.50000 
E) Objective Function: 
Total structure Volume: 202.499207 
Node Normal Stress Bending Stress Total 
af 02:0 4444.4 4444.4 
Ps 07.0 4444.4 4444.4 
3 0.0 4444.4 4444.4 
z| 0.0 4444.4 4444.4 
=] 0.0 4444.4 4444.4 
F) Boundary Conditions: 
Node X-Displ Y-Displ Slope 
1 i iE 1 
5 0 0 0 
G) Solution Vector: 
Node X-Displ Y-Displ Slope 
1 0.000000E+00 0.000000E+00 0.000000E+00 
2 -0.624985E-02 -0.273190E-09 0.111109E-02 
3 -0.249994E-01 -0.109276E-08 0.222218E-02 
4 -0.562488E-01 -0.245871E-08 0.333328E-02 
5 -0.999979E-01 -0.437105E-08 0.444438E-02 
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Validation #4 


OPTIMIZATION SOLUTION 


Arch Angle 90.000 Youngs Modulus: 30000000.0 
Arch Radius: 45.000 Yield Strength: 5§2000.0 
No of Design Var: 10 No of Elements: q 
B) Derived Constants: 
No of System Nodal Points... 2) 
No of Degrees of Freedonm.... i) 
Mengch p@r Element.......... ic) 9 DOL 
Number of Iterations........ 
C) Structure Loading: 
he 0.0000 
ne so 6 1000.0000 
a 5 A 0.0000 m 
PE) ics oo 0. 0 49) 008) 6 Suse eres 0.0000 
D) Elemental Dimensions and Stress Distribution: 
Node Height Base Length Area 
al 3.00000 1.50000 ool 3 4.50000 
Z 3.00000 1.50000 75613 4.50000 
3 3200000 i 50000 e561 3 4.50000 
4 30.0000 1.50000 io. 3 4.50000 
5 B00000 50000 iio 2o1 3 4.50000 
E) Objective Function: 
Total structure Volume: 316.046356 
Node Normal Stress Bending Stress Total 
1 2179 19996.8 20214.8 
2 20r. 4 WS 47 5) 5 18676.8 
3 754.1 14141.1 14295 .2 
4 83.4 Hoss .6 pea? .0 
5 a3... 3 OF. 0 43.3 
F) Boundary Conditions: 
Node X-Displ Y-Displ Slope 
1 i ah 1 
5 0 0 0 
G) Solution Vector: 
Node X-Displ Y-Displ Slope 
1 0.000000E+00 0.000000E+00 0.000000E+00 
2 -0.654528E-01 0.131494E-01 0.750557E-02 
3 -0.223473E+00 0.118865E+00 0.138688E-01 
4 -0.381495E+00 0.355493£+00 0.181207E-01 
5 -0.446951E+00 0.684692E+00 0.196139E-01 
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Validation #4 


OPTIMIZATION SOLUTION 


Arch Angle 90.000 Youngs Modulus: 30000000.0 
Arch Radius: 45.000 Yield Strength: 52000.0 
No of Design Var: 14 No of Elements: 6 
B) Derived Constants: 
No of System Nodal Points... zu 
No of Degrees of Freedom.... 20 
Length per Element...... suoueee 11.7474 
Number of Iterations........ i 
C) Structure Loading: 
9 Ro 5 caer 0.0000 
|; (rrr ei ee ErermietAn sy 6 Go okt ¢ 1000.0000 
BOM oo: 6. « oibusre 4 os lee secepaeer annus eee ee 0.0000 a 
PA 6.0. 2 eerste ye ae 6 aurea eee 0.0000 : 
D) Elemental Dimensions and Stress Distribution: 
Node Height Base Length Area 
a 3.00000 1.50000 11.74736 4.50000 
iz 3.00000 1.50000 74730 4.50000 
3 3.00000 1.50000 1.74736 4.50000 
4 3.00000 1.50000 11.74736 4.50000 
5 3.00000 1.50000 1) .74736 4.50000 
6 3.00000 1.50000 Li. 747 30 4.50000 
% 3.00000 1.50000 11.74736 4.50000 
E) Objective Function: 
Total structure Volume: 317.178711 
Node Normal Stress Bending Stress Total 
ut 220.3 19988.4 202087 
2 Zz.,.6 19310.1 19522.9 
3 190.9 Lass 3:2 17506.0 
4 529 14139.8 14295.7 
5 POS 9999.2 10109.4 
6 Sia 5176.4 5233.4 
u 29.0 Oz 2902 
F) Boundary Conditions: 
Node X-Displ Y-Displ Slope 
1 1 1 at 
” 0 0 0 
G) Solution Vector: 
Node X-Displ Y-Displ Slope 
1 Q.000000E+00 0.000000E+00 0.000000E+00 
2 -0.300317E-01 0.404077E-02 0.512948E-02 
3 -0.112085E+00 0.381154E-01 0.991003E-02 
4 -0.224178E+00 0.124215E+00 0.140157E-01 
5 -0.336278E+00 0.270393E+00 0.171665E-01 
6 -0.418343E+00 0.468604E+00 0.191473E-01 
7 -0.448382E+00 0.696858E+00 0.198230E-01 
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Validation #4 


OPTIMIZATION SOLUTION 





A) Problem Parameters: 
Arch Angle 90.000 Youngs Modulus: 30000000.0 
Arch Radius: 45.000 Yield Strength: 52000.0 
No of Design Var: 18 No of Elements: 8 
B) Derived Constants: 
No cof System Nodal Points... 9 
No of Degrees of Freedom.... 27 
Length per Element.......... 8.8215 
Nuembér of Iterations........ 1 
C) Structure Loading: 
-.... 004 Sr 0.0000 
5 o 4 oa are 1000.0000 
ere 0.0000 _ 
5 6 ee 0.0000 : 
D) Elemental Dimensions and Stress Distribution: 
Node Height Base Length Area 
ne 3.00000 1.50000 8.82154 4.50000 
2 3.00000 1.50000 8.82154 4.50000 
5 3.00000 1.50000 8.82154 4.50000 
4 2.00000 1.50000 8.82154 4.50000 
= 3.00000 1.50000 8.82154 4.50000 
6 3.00000 1.50000 8.82154 4.50000 
7 3.00000 7.590000 8.82154 4.50000 
8 3.00000 1. 500100 8.82154 4.50000 
9 3.00000 1.50000 8.82154 4.50000 
E) Objective Function: 
Total structure Volume: 317.575592 
Node Normal Stress Bending Stress Total 
i! Cok. e 19983.1 20204.3 
2 217.0 19602.4 19819.4 
3 204.4 18467.7 10072. . 
4 184.0 166220 16806.6 
5 1561.06 Manso 14294.6 
6 123.0 Li109 23 1i232..4 
“f 84.7 AOS. 0 TiSd 7 
8 43.2 3901.8 3944.9 
9 21.9 Ore 22.0 
F) Boundary Conditions: 
Node X-Displ Y-Displ Slope 
1 1 1 a 
9 0 0 0 


ES) © 


G) 


Solution Vector: 


Node 


O@M~ INN & WHF 


Lal 


OF 
-Q. 
-0. 
~0. 


-0 


~0. 
.383113E+00 


-0 


-0. 
-Q. 


X-Displ 

OO0000E+00 
170799E-O1 
657226E-01 
138527E+00 


.224410E+00 


310298E+00 


431768E+00 
448854E+00 


-000000E+00 
.109773E-01 
.428386E-01 
~.924699E-01 
~-155017E+00 
-224357E+00 
-293705E+00 
.356268E+00 
-405923E+00 
.437804E+00 
.448790E+00 


Ooooo oaoo°]en 


pan EE an i we Tl we Bl oe OE o> i oe i we Et we oe i we) 


Y~-Displ 


-000000E+00 
~174758E-02 
-165686E-01 
-555485E-01 


.126097E+00 
.230817E+00 
-367110E+00 
-527570E+00 
.701109E+00 


-000000E+00 
-916237E-03 
.861779E-02 
~292281E-01 
.676091E-01 
.126884E+00 
.208132E+00 
.310279E+00 
~430208E+00 
-563056E+00 
-702697E+00 
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oOoooeo oa0o0c a 


oOoooooooo0°0 Oo 


Slope 


-0O00000E+00 
- 388006E-02 
-761158E-02 
~-110510E-01 


-140661E-01 
-165408E-01 
.183798E-01 
~195124E-01 
0. 


198948E-01 


.000000E+00 
sll S6SEs02 
~-615142E-02 
-903846E-02 
~117035£E-01 
-L140808E-01 
peel SE=01 
-177456E-01 
~-189426E-01 
~196729E-01 
-199184E-01 


Validation #4 


OPTIMIZATION SOLUTION 


Arch Angle 90.000 Youngs Modulus: 30000000.0 
Arch Radius: 45.000 Yield Strength: 5200020 
No of Design Var: 26 No of Elements: ali 
B) Derived Constants: 
No of System Nodal Points... 13 
No of Degrees of Freedom.... 39 
Lemogthwper Element.......... 5.8863 
Mempe@r of Iterations........ 
C) Structure Loading: 
RETIRE , otns 6 5 ce ee ec eee 0.0000 
So 5 racer erernr eran 1000.0000 
re 0.0000 a 
9 oe 0.0000 
D) Elemental Dimensions and Stress Distribution: 
Node Height Base Length Area 
i a. 00000 1. S0000 5.88628 4.50000 
2 3.00000 1.50000 3.88628 4.50000 
3 3.00000 1.50000 5.88628 4.50000 
4 3.00000 1720000 5.88628 4.50000 
5 3.00000 1. 50000 5.88628 4.50000 
6 3.00000 10 O00 5.88628 4.50000 
i 3.00000 1.59000 5.88628 4.50000 
8 3.00000 ie. 000 2.00028 4.50000 
9 3.00000 1; 50000 5.88628 4.50000 
10 3.00000 1. SO0IG 5.88628 4.50000 
tee 3.00000 1.50000 >ooe025 4.50000 
a2 3.00000 1.50000 Sogo 25 4.50000 
13 3.00000 1.50000 5.88628 4.50000 
E) Objective Function: 
Total structure Volume: 317.859253 
Node Normal Stress Bending Stress Total 
1 7a Pid Ra 997 3.5 20195. 2 
2 219-8 19805.9 200:25.8 
3 214.2 19299.1 Gis, 3 
4 204.9 18461.6 ECO 6..S 
5 Po2 1 17308.0 7 002 
6 17-6.. 0 15857.4 16033.4 
y 6, 9 ley 3523 a2 2.3 
8 Mes >< 2 2 7 Le0 e306). 2 
9 eel. 1 9998.4 O19. 5 
10 S50 TOS 3 co 113841 
a1 57.4 D1) Oe 3755.7 
12 742 a ZO1027 2639.8 
AS a. O21 4.7 
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F) Boundary Conditions: 


Node X-Displ 
ak 1 
13 0 
G) Solution Vector: : 
Node X-Displ 
1 0.000000E+00 
2 -0.764866E-02 
3 -0.300750E-01 
4 -0.657529E-01 
5S =~0.112252E+00 
6 -0.166406E+00 
7 -0.224523E+00 
8 -0.282643E+00 
9 -0.336806E+00 
10 -0.383318E+00 
11 -0.419010E+00 
12 -0.441447E+00 
13 -0.449100E+00 


oOoo°ooo°o0o°o°o°;co 0° 


Y-Displ 
i 


0 


Y-Displ 


-000000E+00 
-544914E-03 
-504941E-02 
.172040E-01 
.401787E-01 
.764064E-01 
~127417E+00 
-193735E+00 
.274838E+00 
.369200E+00 
~-474387E+00 
.587230E+00 
.704036E+00 


oOoooooeooo°ooo eo 


Slope 


Slope 


.000000E+00 
-260170E-02 
-515929E-02 
-762896E-02 
~996841E-02 
~121375E-01 
~140991E-01 
HL Oe 7 bod 
»-172696E-01 
~-184241E-01 
~192632E-01 
~197725E-01 
~199432E-01 
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Validation #5 


OPTIMIZATION SOLUTION 





Arch Angle 90.000 Youngs Modulus: 30000000.0 
Arch Radius: 32.000 Yield Strength: 52000.0 
No of Design Var: 10 No of Elements: 4 
B) Derived Constants: 
No of System Nodal Points... 5 
No of Degrees of Freedonm.... iS 
eength per Element........-.. 12.4059 
Number of Iterations........ a 
C) Structure Loading: 
MEEPS) els 2 ose s 6 6 tie es es : 0.0000 
<2 0 6 oes err -5000.0000 
a Meee getsia et cuneate aes 0.0000 s 
MMM) 5s) «oe 4 9) ose 6 ose Woe 8 0 ee 0.0000 A 
D) Elemental Dimensions and Stress Distribution: 
Node Height Base Length Area 
it 3.00000 i750000 12.48578 4.50000 
2 3.00000 1.50000 12.48578 4.50000 
3 3.00000 1.50000 12.48578 4.50000 
4 3.00000 1.50000 12.48578 4.50000 
5 3.00000 i 50000 12.48578 4.50000 
E) Objective Function: 
Total structure Volume: 224.744095 
Node Normal Stress Bending Stress Total 
tk e231. 2 O70 M2 ai i2 
2 1278.8 23350 3616.9 
3 ey 3. 2 aT Ee i Gas 13245.0 
4 NO73.8 1043.2 2116.9 
5 g27.6 24/25)..3 25652.9 
F) Boundary Conditions: 
Node X-Displ Y-Displ Slope 
i 1 al 0 
5 I 0 i 
G) Solution Vector: 
Node X-Displ Y-Displ Slope 
1 0.000000E+00 0.000000E+00 0.302544E-02 
2 -0.301621E-01 0.547718E-02 0.131377E-02 
3 -0.265139E-01 0.237552E-02 -0.205875E-02 
4 -0.530910E-02 -0.302736E-01 -0.357488E-02 
5 0.000000E+00 ~-0.589432E-01 0.000000E+00 
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Validation #5 


OPTIMIZATION SOLUTION 


—_= = «= op oe ap oe oe oe om oe ae w= om oe oe oe oe oe we oe ew om oe oe om oe © © oe © oe © om oo oe ow © oo oe oe ew oe oe oe oe oe oe oe oe oe oe oe ee ee ee oe 


A) Problem Parameters: 
Arch Angle 90.000 Youngs Modulus: 30000000.0 
Arch Radius: 32,000 Yield Strength: 52000.0 
No of Design Var: 14 No of Elements: 6 
B) Derived Constants: 
No of System Nodal Points... 7 
No of Degrees of Freedon.... Pash 
Length per Element.......... 6.3537 
Number of It@rations........ i: 
C) Structure Loading: 
) PEP ee re 6 5 6 bad o 0.0000 
i Sara er I 5 RN -5000.0000 
PM 5 ws 6S 5 % sre raepee ate ene ge eee 0.0000 Z 
PA 0.) 0 cohol’ oie coneouatce see) «cb emei en nen 0.0000 z 
D) Elemental Dimensions and Stress Distribution: 
Node Height Base Length Area 
1 3.00000 Ie Ota, 8.35368 4.50000 
2 3.00000 50006 8.35368 4.50000 
3 3.00000 1.50000 8.35368 4.50000 
4 3.00000 20000 8.35368 4.50000 
5 3.00000 1.50006 a725560 4.50000 
6 3.00000 1.50000 8.35368 4.50000 
7 3.00000 1.50000 8.35368 4.50000 
E) Objective Function: 
Total structure Volume: §225.549301 
Node Normal Stress Bending Stress Total 
uh 1194.9 0.0 1194.9 
2 1247.4 9411.2 10658 .6 
3 PS 0OR 2 13334.9 14643.1 
4 gratis By. 11503.9 LZ 
5 1164-2 4042.8 $207 
6 969.3 8539.9 SLU) 5 7 
7 653.4 25300.c 26240.2 
F) Boundary Conditions: 
Node X-Displ Y-Displ Slope 
nf 1 nl 0 
ui ih 0 i 
G) Solution Vector: 
Node X-Displ Y-Displ Slope 
1 0.000000E+00 0.000000E+00 0.316688E-02 
2 -0.238606E-01 0.280572E-02 0.229335E-02 
3 -0.340200E-01 0.662210E-02 0.182085E-03 
4 -0.276224E-01 0.125100E-02 -0.212342E-02 
5 -0.128425E-01 -0.185793E-01 -0.356644E-02 
6 -0.176538E-02 -0.461115E-01 -0.314903E-02 
7 Q0.000000E+00 -0.613415E-01 0.000000E+00 
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wa~veation #5 


OPTIMIZATION SOLUTION 





A) Problem Parameters: 
Arch Angle 90.000 Youngs Modulus: 30000000.0 
Arch Radius: 32.9000 Yield Strength: 52000720 
No of Design Var: 18 No of Elements: 8 
B) Derived Constants: 
No of System Nodal Points... 9 
No of Degrees of Freedon.... aa 
Length per Element.......... G.27 31 
Number of Iterations........ 1 
C) Structure Loading: 
5 6 GSoC nsee mae 0.0000 
So 6 rr ; -5000.0000 
so 5 Se 5 0.0000 7 
6 A ere 0.0000 2 
D) Elemental Dimensions and Stress Distribution: 
Node Height Base Length Area 
1 3.00000 1 90000 6227 50 4.50000 
2 3.00000 i 90000 6227510 4.50000 
3 3.00000 i. 50000 627310 4.50000 
4 3.00000 i. 509000 6.27310 4.50000 
5 3.00000 i. 50000 6.27310 4.50000 
6 3.00000 1.50000 6227310 4.50000 
u 3.00000 50000 627310 4.50000 
8 3.00000 1.50000 27310 4.50000 
9 3.00000 1.50000 6727 310 4.50000 
E) Objective Function: 
Totaimesccucture Volume: 225.831528 
Node Normal Stress Bending Stress Total 
1 1175.4 GO. 1175.4 
z 1 PACAP A 1509'.3 8731.8 
3 eZ o2 3 11997 :.3 13289 .6 
4 1312.4 1329126 14604.0 
=) 262.1 sac. 3 26245 
6 1202.5 6224.5 1427.9 
7 1076.8 1865.4 2942.1 
8 909.6 L266)... 3 173525 °.8 
9 816.4 256155 .1 26431.5 
F) Boundary Conditions: 
Node X-Displ Y-Displ Slope 
1 1 aL 0 
9 iL 0 il 
G) Solution Vector: 
Node X-Displ Y-Displ Slope 
1 0.000000E+00 0.000000E+00 0.321682E-02 
2 -0.190172E-01 0.162605E-02 0.269341E-02 
3 -0.314948E-01 0.513368E-02 0.133378E-02 
4 -0.342107E-01 0.627359E-02 -0.428882E-03 
5S -0.280319E-01 0.848441E-03 -0.214589E-02 
6 -0.170219E-01 -0.129808E-01 -0.337032E-02 
7 -0.6540S51E-02 -0.331005E-01 -0.367416E-02 
8 -0.742209E-03 -0.529373E-01 -0.266477E-02 
9 0.000000E+00 -0.622147E-01 0.000000E+00 
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Validation #5 


OPTIMIZATION SOLUTION 


Arch Angle 90.000 Youngs Modulus: 30000000.0 
Arch Radius: 32.000 Yield Strength: §2000.0 
No of Design Var: 22 No of Elements: 10 
B) Derived Constants: 
No of System Nodal Points... aI 
No of Degrees of Freedonm.... 33 
Length per Element.......... 5.0214 
Number of Iterations........ 
C) Structure Loading: 
IX oc ee ee + Wleene ¢. sone memeneein cece 0.0000 
Grrr ee ee ee -5000.0000 
BM oe 0 6% ereste eueiel eleilee) otc eer 0.0000 ~ 
ft er no ee co OO oe 0.0000 - 
D) Elemental Dimensions and Stress Distribution: 
Node Height Base Length Area 
1 3.00000 i-30000 a202138 4.50000 
Z 3.00000 1.50000 S02 138 4.50000 
3 3.00000 1.50000 5.02136 4.50000 
4 3.00000 1.50000 5202136 4.50000 
5 3.00000 1.50000 S. 02138 4.50000 
6 3.00000 1.50000 5202138 4.50000 
7 3.00000 1.50000 3.02138 4.50000 
8 3.00000 1.50000 5.02136 4.50000 
9 3.00000 1.50000 02 oe 4.50000 
10 3.00000 1.50000 502136 4.50000 
tae 3.00000 1.50000 >.0Ziag 4.50000 
E) Objective Function: 
Total structure Volume: 225.962219 
Node Normal Stress Bending Stress Total 
ik ino s". 3 0.0 os. 3 
2 1204.7 6225.1 7429.8 
3 We7 2.0 10546.0 ee? . 9 
4 1307.9 TZ8s6.2 14164.1 
5 5 ic) oe ae | 13098.8 14410.6 
6 P2605 22 11268:.0 eZ oie 2 
7 223 7408.7 8631.5 
8 PisZac 1616.0 2748.9 
9 1014.8 5967.4 6982.2 
10 871.6 1$155.0 16026.7 
11 794.2 25720.4 26914. 7 
F) Boundary Conditions: 
Node X-Displ Y-Displ Slope 
1 ah ah 0 
11 1 0 2 
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G) 


Solution Vector: 


Node 


Me 


POW OAINMN AWN he 


X-Displ 


-000000E+00 
~156551E-01 
-277397E-01 
~339727E-01 
-337863E-01 
»282253E-01 
-195378E-01 
-105138E-01 
-366515E-02 
~355416E-03 
- O0O0000E+00 


oooo°oeo 


=o 


Y¥-Displ 


.000000E+00 
.103676E-02 
ao tes 5402 
~607021E-02 
~569723E-02 
-660388E-03 
-0. 
-0. 
=0. 
-0. 
-626255E-01 


983642E-02 
249420E-01 
419486E-01 
564151E-01 
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(= a 3 @ Xa) 


Slope 


~324004E-02 
~289272E-02 
-195701E-02 
-651332E-03 
-/96782E-03 
~215629E-02 
~319832E-02 
~-370184E-02 
-345906E-02 
~ 228057E-02 
-000000E+00 


Validation #5 


OPTIMIZATION SOLUTION 


Arch Angle 90.000 Youngs Modulus: 30000000.0 
Arch Radius: 32.000 Yield Strength: 5200020 
No of Design Var: 26 No of Elements: TZ 
B) Derived Constants: 
No of System Nodal Points... 13 
No of Degrees of Freedom.... 39 
Length per Element.......... 4.1858 
Number of Iterations........ 
C) Structure Loading: 
Pere ar Ape Oe C. 8.0 6.6 bc 0.0000 
earner a ee a er 58 -5000.0000 
BM eb ecsce- es vel'ss ease? cc neue on ene eee 0.0000 - 
BA esos o sceeie Gece pasion ensue ee tome teen 0.0000 - 
D) Elemental Dimensions and Stress Distribution: 
Node Height Base Length Area 
1 3.00000 1.59000 4.18580 4.50000 
2 3.00000 1.50060 4.18580 4.50000 
3 3 .C0GGC 1.59000 4.18580 4.50000 
4 3.00000 1.50000 4.18580 4.50000 
=) 3.00000 1.50000 4.18580 4.50000 
6 3.00000 1.50000 4.18580 4.50000 
i 3.00000 1.59000 4.18580 4.50000 
8 3.90000 1.50000 4.18580 4.50000 
9 3.00000 1.50008 4.18580 4.50000 
10 3.00000 1.50000 4.18580 4.50000 
Le 3.00000 1.50000 4.18580 4.50000 
2 3.00000 py ay 10N0 10) 4.18580 4.50000 
is 3.00000 iS C000 4.18580 4.50000 
E) Objective Function: 
Total structure Volume: 226.033264 
Node Normal Stress Bending Stress Total 
1 115550 OF Liss) 
Z 11915 5308.7 6500.2 
3 125352 9309.8 1056327 
4 1294.8 1193'570 13229.8 
=) 131356 1313983 14452.9 
6 1309.9 1290222 14212571 
7 126320 112275 12523 
8 P2352. 8144.1 9379.9 
9 11665 3704.7 4871.2 
10 uO Pee ae ZO L510 309223 
ps 969.7 8917.1 9886.8 
iz 845.5 16683 .3 LI7Z266 
3 119.5 25717%.3 26556.9 
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F) Boundary Conditions: 


Node X-Displ 
ik 1 
13 1 
G) Solution Vector: 
Node X-Displ 
1 0.000000E+00 
2 -0.132523E-01 
3 -0.243569E-01 
4 -0.317555E£-01 
5 -0.347043E-01 
6 -0.332895E-01 
7 -0.283308E-01 
8 -0.211878E-01 
9 -0.134844E-01 
10 -0.678007E-02 
11 -0.221779E-02 
12 -0.181929E-03 
13 0.000000E+00 


Y-Displ 


-000000E+00 
~-707096E-03 
-274126E-02 
-506415E-02 
-631460E-02 
~514815E-02 
-557844E-03 
~-785492E-02 
~-196869E-01 
~336374E-01 
.475236E-01 
~584105E-01 
-628492E-01 


ar 


Slope 


-325259E-02 
~-300569E-02 
~232580E-02 
2L33772E=-02 
~171546E-03 
~103961E-02 
»- 216186E-02 
-306281E-02 
~- 361389E-02 
~-369247E-02 
~-318404E-02 
~-198409E-02 
-000000E+00 


Validation #5 


OPTIMIZATION SOLUTION 


Arch Angle 90.000 Youngs Modulus: 30000000.0 
Arch Radius: 22.2000 Yield Strength: 52000.0 
No of Design Var: 30 No of Elements: 14 
B) Derived Constants: 
No of System Nodal Points... i> 
No of Degrees of Freedon... 45 
Length per Element. .7 1.2.2... 3.5885 
Number of Iterations........ 1 
C) Structure Loading: 
| > rr ee ar erry st ob 6 0.0000 
BY ous 0 0 sree 5 she eieneh a en -5000.0000 
ty Pr re kc sc 0.0000 z 
PA . « wis oa 0 eae 00s) ere eco) eaten 0.0000 . 
D) Elemental Dimensions and Stress Distribution: 
Node Height Base Length Area 
il 3.00000 1.50000 3.56051 4.50000 
az 3.00000 1.50000 3.50051 4.50000 
3 3.00000 1.50000 3.58851 4.50000 
4 3.00000 1.50000 sSea0do1 4.50000 
5 3.00000 1.50000 3256651 4.50000 
6 3.00000 1.50000 3.506 OL 4.50000 
i 3.00000 1.50000 3750851 4.50000 
8 3.00000 1.50000 Ei clelepk 4.50000 
9 3.00000 1.50000 3256851 4.50000 
10 3.00000 1.50000 Soe OL 4.50000 
Esl 3.00000 1.50000 3250051 4.50000 
a2 3.00000 1.50000 S.50601 4.50000 
a3 3.00000 1.50000 3.56651 4.50000 
14 3.00000 1.50000 S566 OL 4.50000 
5 3.00000 1.50000 3. 906 511 4.50000 
E) Objective Function: 
Total structure Volume: 226.076035 
Node Normal Stress Bending Stress Total 
1 1149.0 0.0 1149.0 
2 61.5 4624.5 5805.9 
x Zsi6. 6 6296.5 gos5.3 
4 1280.5 10970.1 2250.6 
5 7306.1 EZ611.. 5 3917.5 
6 13.5% 3 3 2.010). 1. 14515.4 
7 1307.9 T2726. 6 14036.5 
8 1284.1 i027 12486.8 
9 1244.2 8641.9 9886.1 
10 1188.6 5076.2 6266.8 
ae i 16:.'0 556.4 1674:.5 
12 1033.5 4866.7 5900.2 
13 935.9 LEiZ3e 1 12058.9 
14 826.5 e133. 9 18960.3 
15 16950 2610.9 26579.9 


S10 


Bimeoundary Conditions: 


G) 


Node 
1 
LS) 


X-Displ 
i} 


ak 


Solution Vector: 


Node 


OF 
-0. 
-0. 
-0. 
-Q0. 
-0. 
-0. 
-0. 
-0. 
m5 7056E=01 
-949819E-02 


-0 
-0 


=0), 
-0. 
~-946979E-04 


-0 


0. 


X-Displ 

O00000E+00 
114679E-0O1 
215484E-01 
291563E-01 
336353E-01 
347835E-01 
328342E-01 
283939E-01 
223430E-01 


456627E-02 
142229E-02 


000000E+00 


oOO9O0O000 0 


Y-Displ 


Y-Displ 


.O00000E+00 
.506369E-03 
-207183E-02 
-410633E-02 
-579362E-02 
-624852E-02 
~467242E-02 
-495743E-03 
-650287E-02 
.161095E-01 
9271627 2E=02 
~398712E-01 
soLZ1IOE=S=O2 
-596559E-01 
~-629827E-01 


aes)! 


Slope 


© 


Slope 


-326008E-02 
-307569E-02 
-256050E-02 
~179229E-02 
-852037E-03 
eld 31 ee 03 
~121096E-02 
-216516E-02 
-295641E-02 
-350346E-02 
»372813E-02 
-355627E-02 
~-291872E-02 
-175218E-02 
-000000E+00 


Validation (#5 


OPTIMIZATION SOLUTION 


A) Problem Parameters: 
Arch Angle 90.000 Youngs Modulus: 
Arch Radius: 32000 Yield Strength: 
No of Design Var: 34 No of Elements: 
B) Derived Constants: 
No of System Nodal Points... 17 
No of Degrees of Freedom.... Syal 
Length per Element.......5... 3.1803 
Number of IGerations........ i 
C) Structure Loading: 
) CP ee oe co 0.0000 
CR MED © cox -~5000.0000 
OMe. 0 ee he wes 6 60s 5 0 6 eee 0.0000 
FAs be ose 0 esses 6% sceunitels. acs 0.0000 
D) Elemental Dimensions and Stress Distribution: 
Node Height Base Length 
al 3.00000 1.50000 3.14033 
2 3.00000 1.50000 3.14033 
3 3.00000 1.50000 3.14033 
é 3.00000 1.50000 3.14033 
5 3.00000 1.59000 3.14033 
6 3.00000 1.50000 3.14033 
7 3.00000 1.50000 3.14033 
8 3.00000 1.50000 3.14033 
9 3.00000 1.50000 3.14033 
10 3.00000 1.50000 3.14033 
net 3.00000 1.50000 3.14033 
M7 3.00000 1.50000 3.14033 
i! 3.00000 1.50000 3.14033 
14 3.00000 1.50000 3.14033 
15 3.00000 1.50000 3.14033 
16 3.00000 1.50000 3.14033 
17 3.00000 1.50000 3.14033 
E) Objective Function: 
Total structure Volume: 226.103821 
Node Normal Stress Bending Stress 
1 1144.4 Or 
Zz SS 4094.6 
3 226.1 7465.0 
4 1266.9 10078.9 
5 29 5.6 tF911..1 
6 earl 7 12943.9 
7 W152 1367.7 
8 1306.0 l2579..8 
9 28a 2 11186.0 
10 PASO) a ak 8999.8 
st 1204.0 6042.2 
2 1146.2 2341.6 
3 OR. 2 Z060.5 
14 998.3 Ws 9! 5 
15 909.5 2826.6 
16 812.0 19079.4 
lee ice net 25831.4 


LZ 


30000 
32 


Tota 


1144. 
5268. 
B69 1. 
TiE345.. 
Ta 06 . 
14255. 
14482. 
136 O)27. 
12470. 
10250. 
7246. 
3487. 
3143. 


Sas 7. 
en sar 
19891. 


26592 


000.0 
000.0 
16 


Area 
4.50000 
4.50000 
4.50000 
4.50000 
4.50000 
4.50000 
4.50000 
4.50000 
4.50000 
4.50000 
4.50000 
4.50000 
4.50000 
4.50000 
4.50000 
4.50000 
4.50000 


1 


mbm WOrFrOW IO AVADONE- A 
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F) Boundary Conditions: 


Node X-Displ 
ak of 
ey a 
Goesolution Vector: 
Node X-Displ 
1 0.000000E+00 
2 -0.100961E-01 
3 -0.192461E-01 
4 -0.266726E-01 
S -0.318426E-01 
6 -0.344882E-01 
7 -0.346060E-01 
8 -0.324379E-01 
9 -0.284326E-01 
10 -0.231915E-01 
11 -0.174003E-01 
12 -0.117521E-01 
13 -0.686387E-02 
14 -0.319243E-02 
15 -0.953191E-03 
16 -0.472050E-04 
17 0Q.000000E+00 


tot 2. 0 te 
ooooo0o0oceoo0o0o0 900 000 


Y-Displ 


.000000E+00 
-376056E-03 
-160606E-02 
~333145E-02 
-503854E-02 
-613847E-02 
-604840E-02 
~-426921E-02 
-455661E-03 
-552506E-02 
-135497E-01 
223 21350E=-01 
-338212E-01 
-444054E-01 
-537569E-01 
-604801E-01 
-630648E-01 


Lo 


Slope 


- © 


Slope 


-326459E-02 
olen? ZE-02 
-271838E-02 
-210622E-02 
~133893E-02 
-471678E-03 
-439423E-03 
sl SS782E-02 
-216707E-02 
- 287140E-02 
-339626E-02 
-368879E-02 
- 3698 38E-02 
EUG) 210) 72 
- 268042E-02 
-156706E-02 
-Q00000E+00 


APPENDIX E 
CASE STUDIES 


OPTIMIZATION #1 


OPTIMIZATION SOLUTION 


A) Problem Parameters: 


Arch Angle : 0.002 Youngs Modulus: 30000000.0 
Arch Radius: 1000000.000 Yield Strength: 52000.0 
No of Design Var: 18 No of Elements: 8 
B) Derived Constants: 
No of System Nodal Points... 9 
No of Degrees of Freedom.... PAT 
Length per Element... .--fee 4.0000 
Number of Iterations........ al 
C) Structure Loading: 
FX. 0.6 Seu es 6.6 6 See Eee concn 2000.0000 
pb Pare rere so 6 ee 0.0000 
|, Pe EPPEr Eye cc 5 cece c 0.0000 = 
ey re rrr Go co nec 0.0000 : 
D) Elemental Dimensions and Stress Distribution: 
Node Height Base Length Area 
1 4.19530 0.41953 4.00000 1.76005 
2 4.01266 0.40127 4.00000 1.61015 
3 3.8169 0.381ig 4.00000 1.45290 
4 3.58695 0.35869 4.00000 i. 28662 
5 3.32962 0.33296 4.00000 1.10877 
6 3.02540 0.30254 4.00000 0.91530 
7 2.64292 0.26429 4.00000 0.69850 
8 2.09722 0.20977 4.00000 0.44004 
9 0.03041 0.03000 4.00000 0.00091 
E) Objective Function: 
Total structure Volume: 33.1263862 
Node Normal Stress Bending Stress Total 
0.0 52002.4 52002.4 
2 0.0 52002.4 5§2002.4 
3 0.0 5200255 52002358 
4 O70 52002.4 5§2002.4 
5 0.0 5200322 §2003.1 
6 0.0 52000.9 52000.9 
7 0.0 5§2001.6 5§2001.6 
8 0.0 51999.4 51999.4 
9 O.1 844.7 844.8 
F) Boundary Conditions: 
Node X-Displ Y-Displ Slope 
1 1 1 1 
9 0 0 0 
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Se 


Pesolution Vector: 


Node 
al 


O@D~IHMN Hh WwW po 


oRokokolokolo Nolo) 


X-Displ 


.0O00000E+00 
.691840E-02 


moe 7399E-01 
.632283E-01 
.114358E+00 
.182449E+00 
-269490E+00 
.379305E+00 
-618095E+00 


oookokokolo More, 


Y=Dispil 


.000000E+00 
»295492E-09 


.119801E-08 
-274074E-08 
-496595E-08 
.793075E-08 
-117209E-07 
-165004E-07 
-268468E-07 


isis: 


Slope 


.O0O0000E+00 
~338401E-02 


-693567E-02 
~-L06947E-01 
~L147211E-01 
-L9LISZE=-01 
-240748E-01 
~301571E-01 
~- 744676E-01 


OPTIMIZATION #la 


OPTIMIZATION SOLUTION 


A) Problem Parameters: 
Arch Angle : 0-002 Youngs Modulus: 
Arch Radius: 1000000.000 Yield Strength: 
No of Design Var: 18 No of Elements: 
B) Derived Constants: 
No of System Nodal Points... 9 
No of Degrees of Freedom.... 27 
Length per Element.......... 4.0000 
Numbér of Iterations......:.. a 
C) Structure Loading: 
PX uo: Stace Piel <ua 10: erento a Re meee eis sts 2000.0000 
BY on 's eceriar eden peuaurasomeue te emma 0.0000 
bs Pa En i Tia EY Or pe 0.0000 
PA. 5 ccs a, eee ee 0.0000 
D) Elemental Dimensions and Stress Distribution: 
Node Height Base Length 
i 4.19499 0.41961 4.00000 
2 4.01244 0.40133 4.00000 
3 3.81144 0.38124 4.00000 
4 3.58705 On 3s 0 4.00000 
5 3.32942 0. 33.306 4.00000 
6 3.02541 0.30254 4.00000 
7 2.64171 0.26456 4.00000 
8 2.09811 0.20981 4.00000 
9 0.10080 0.03000 4.00000 
E) Objective Function: 
Total structure Volume: 33.148262 
Node Normal Stress Bending Stress 
1 0.0 Sz001.2 
2 0.0 SZe0O1 22 
6 0.0 52000.0 
4 0.0 51998.6 
5 0.0 52003.4 
6 0.0 52000.0 
u 0.0 51997.4 
8 0.0 5197055 
9 0.0 00 
F) Boundary Conditions: 
Node X-Displ Y-Displ Slope 
1 1 1 iL 
9 0 0 0 
G) Solution Vector: 
Node X-Displ Y-Displ Slope 
1 0.Q000000E+00 0.000000E+00 0.000000E+00 
2 0.691870E-02 0.295505E-09 -0.338415E-02 
3 0.277409E-01 0.119806E-08 -0.693591E-02 
4 0.632300E-01 0.274082E-08 -0.106947E-01 
5 0.114360E+00 0.496604E-08 -0.147212E-01 
6 0.182452E+00 0.793089E-08 -0.191155E-01 
7 0.269496E+00 0.117212E-07 -0.240760E-01 
8 0.379315E+00 0.165008E-07 -0.301576E-01 
9 0.607042E+00 0.263666E-07 -0.703189E-01 


1 516 


30000000.0 
52000.0 


CODORPRRREH 


Omh Ol HDONW 
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- OPTIMIZATION #2 


OPTIMIZATION SOLUTION 


Arch Angle 90.000 Youngs Modulus: 30000000.0 
Arch Radius: 32.000 Yield Strength: 52000.0 
No of Design Var: 18 No of Elements: 8 
B) Derived Constants: 
No of System Nodal Points... 9 
No of Degrees of Freedonm.... PAT 
Pemeeen per Element.......... 622731 
Bimmoercr Of Iterations......e- B 
C) Structure Loading: 
6 eer 0.0000 
5 6 6 -2000.0000 
1-505 5 er 0.9000 = 
4 6 6 5 Se 0.0000 
D) Elemental Dimensions and Stress Distribution: 
Node Height Base Length Area 
1 1.95274 Peg5611 6.27310 3 78L977 
2 1.94555 1.93241 Gine 310 3.7 3960 
3 1.89934 1.90956 Geran 310 3.62691 
4 1.32039 179:2025 6527310 3.68763 
5 YSZ 03 ets 27 6.21310 3.68698 
6 1.98325 1.92011 62727310 3.668326 
7 1.91849 1.92476 6.2/310 3.69264 
8 Lagze2 2 J. 31786 62273510 3.69422 
9 1.96640 1.96000 6.27310 3sp413 
E) Objective Function: 
fecal structure Volume: 186.151276 
Node Normal Stress Bending Stress Total 
ik 52a 51491.0 S202. 2 
2 5 19.3 51498.4 520 Mie 7 
3 50720 52507228 720146 
4 448.8 45091.9 45540.7 
=) ciel oy 38340.4 3872260 
6 300.0 Z0197.7 30497.7 
zu zZ0>29 20743 .2 20949.1 
8 104.9 NOS2Z7.46 1 O62 27 
9 307.9 0.3 op Lay 
F) Boundary Conditions: 
Node X-Displ Y-Displ Slope 
1 1 1 ak 
9 0 0 0 
G) Solution Vector: 
Node X-Displ Y-Displ Slope 
1 0.Q000000E+00 0.000000E+00 0.000000E+00 
2 0.345895E-01 -0.351714E-02 -0.110489E-01 
3 0.134858E+00 -~0.340464E-01 -0.222560E-01 
4 0.287683E+00 -0.115848E+00 -0.328233E-01 
5 0.469409E+00 -0.265100E+00 -0.419063E-01 
6 0.651576E+00 -0.487184E+00 -0.493731E-01 
7 0.806243E+00 -0.776659E+00 -0.549243E-01 
8 0.909657E+00 -0.111768E+01 -0.583234E-01 
9 0.945957E+00 -0.148635E+01 -0.594191E-01 


sR s) 


OPTIMIZATION 42a 


OPTIMIZATION SOLUTION 


Arch Angle 90.000 Youngs Modulus: 30000000.0 
Arch Radius: 32.000 Yield Strength: 52000.0 
No of Design Var: 18 No of Elements: 8 
B) Derived Constants: 
No of System Nodal Points... 9 
No of Degrees of Freedonm.... og | 
Length per Element.......... Oe 31 
Number of Iterations........ 1 
C) Structure Loading: 
)) Pe rere Se Se yO Oe 0.0000 
PPP ens OO eee -2000.0000 
it PP re ra a oe 0.0000 = 
FA osinerene  clleteeees « o-< 2 . « seen 0.0000 7 
D) Elemental Dimensions and Stress Distribution: 
Node Height Base Length Area 
a 3.91339 0.49197 6.27310 1392527 
2 4.10780 0.43837 632727310 1.80073 
3 4.11506 Oa 125 1 6.27310 1 69237 
4 3.97207 0.39721 6.27310 1 Se 3 
5 3623 0.39033 6.27310 1.44130 
6 2.42308 0.70760 6.27310 1.71458 
7 2.32561 O%53330 6.27310 lv22503 
8 1.30739 0.84851 627) 3.10 1.10933 
9 0.80670 1.49956 6.27310 1.20969 
E) Objective Function: 
Total structure Volume: 77.179106 
Node Normal Stress Bending Stress Total 
i 1033.8 50968.7 52002.6 
2 1084.1 50916.9 52000.9 
3 1085.9 $0913.5 §$1999.4 
4 1048.9 50949.6 51998.5 
5 976.5 51020:..7 $1997.2 
6 645.0 51350.-6 51995 .6 
7 616.7 oy Hac en | $1995.9 
8 3501.0 5$1653.2 SZ200322 
9 or 0 0.2 162.2 
F) Boundary Conditions: 
Node X-Displ Y-Displ Slope 
u 1 1 af 
9 0 0 0 
G) Solution Vector: 
Node X-Displ Y-Displ Slope 
1 0.Q000000E+00 0.000000E+00 0.000000E+00 
2 0.165765E-01 -0.185682E-02 -0.530039E-02 
3 0.640268E-01 -0.164901E-01 -0.104791E-01 
4 0.136735E+00 -0.556094E-01 -0.157515E-01 
S 0.226861E+00 -0.129851E+00 -0.213260E-01 
6 0.324976E+00 -0.249654E+00 -0.277857E-01 
7 QO.421348E+00 -0.430238E+00 -0.369135E-01 
8 0.499870E+00 -0.689424E+00 -0.482073E-01 
9 0.534109E+00 -0.103740E+01 -0.595043E-01 
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OPTIMIZATION #3 


OPTIMIZATION SOLUTION 


Arch Angle : 90.000 Youngs Modulus: 30000000.0 

Arch Radius: 32.000 Yiete Strength: $2000.0 

No of Design Var: 18 No of Elements: 8 
B) Derived Constants: 

No of System Nodal Points... 9 

No of Degrees of Freedom.... 27 

bength per Element.......... CG. 273L 

Number of Iterations........ 1 


C) Structure Loading: 


oo 5 ORS Seer 2000.0000 
0 34 ee 0.0000 
v5 9 6 Siar eer 0.0000 a 
5 6 errr 0.0000 = 
D) Elemental Dimensions and Stress Distribution: 
Node Height Base Length Area 
Bf 3.66465 O72 100 6@27310 2.01972 3 
2 3.83344 0.40649 6.2/310 lee b2S 
3 3.34200 0.41258 6.27310 1337663 
4 2.46546 0.54860 6227310 1.35256 
5 1.99742 Ga 2o7 2 6227310 ee al lollolak 
6 1.48940 0.58247 6.27310 0.86752 
7 lqei374 0.66160 6.27310 OF 8Cs01 
8 0.58386 0.48072 6.27310 0.28068 
9 7) 7 Sal 7603606 6.27310 1.43612 
E) Objective Function: 
Total structure Volume: 56.657707 
Node Normal Stress Bending Stress Total 
il 918 Soo 5.5 $1992.6 
2 249.2 517 4522 51994.4 
3 552.4 §1445.1 51997.4 
4 817-5 51180.4 51997.9 
5 126759 507 306 S1996r5 
6 19177 50087.6 5199503 
7 2290.0 2999178 32201 6 
8 6955-2 450251 51980. 3 
9 1386.0 O50 1386.0 
F) Boundary Conditions: 
Node X-Displ Y-Displ Slope 
at il 1 
9 0 0 0 


G) Solution Vector: 

Node X-Displ Y-Displ Slope 
~O0O0000E+00 0.000000E+00 0.000000E+00 
~-185814E-01 -0.180717E-02 -0.574511E-02 
~-720148E-01 -0.179297E-01 -0.117837E-01 
-158405E+00 -0.639459E-01 -0.190296E-01 
-275821E+00 -0.160028E+00 -0.286801E-01 
-416853E+00 -0.331363E+00 -0.409667E-01 
-559677E+00 -0.597638E+00 -0.537583E-01 
-678467E+00 -0.986551E+00 -0.722173E-01 
~725130E+00 -0.145498E+01 -0.764554E-01 


WOMnA WSL WNP 
oOooOoOOoO0O000 0 


Lo 


OPTIMIZATION #3a 


OPTIMIZATION SOLUTION 


Arch Angle 90.000 Youngs Modulus: 
Arch Radius: 32.000 Yield Strength: 
No of Design Var: 18 No of Elements: 
B) Derived Constants: 
No of System Nodal Points... 9 
No of Degrees of Freedonm.... 27 
Length per Element.......... Gric7 SL 
Number of Iterations........ 2 
C) Structure Loading: 
EX.i sles ss «0s spas ene cn On momen oe 2000.0000 
1 Perera herr Te ety) ae 0.0000 
BM. acd 0 ai-aca eueiegece cu smedcmera) aeons 0.0000 
fr rr rea 6 Sa A 0.0000 
D) Elemental Dimensions and Stress Distribution: 
Node Height Base Length 
iL 4.19827 0.41983 622/310 
2 3.90887 0.39089 6.27310 
3 358501 0.35858 6.27310 
4 3.250 1 Ons2237 6.27310 
5 2.81849 Oze165 6227310 
6 2.36310 0.23631 6.27310 
zl 1.84434 0.18443 6.27310 
8 0.99509 0.18101 6.27310 
9 0.60305 0.06348 6223310 
E) Objective Function: 
Total structure Volume: 41..3410122 
Node Normal Stress Bending Stress 
1 lit §1891.6 
2 2541 §1749.1 
3 892.4 §1411.0 
4 1064.0 50938.9 
5 LAG S02317.°5 
6 2903.0 49040.2 
7 §406.0 46591.3 
8 10838.0 41165.8 
9 §1995.7 0.0 
F) Boundary Conditions: 
Node X-Displ Y-Displ Slope 
1 aL uk aL 
9 0 0 0 
G) Solution Vector: 
Node X-Displ Y-Displ Slope 
1 0.000000E+00 0.000000E+00 0.000000E+00 
2 0.17360S5E-01 -0.168480E-02 -0.536752E-02 
3 0.677523E-01 -0.168806E-01 -0.111584E-01 
4 0.148066E+00 -0.596161E-01 -0.175022E-01 
§ 0.251532E+00 -0.144153E+00 -0.246090E-01 
6 0.367672E+00 -0.284910E+00 -0.328323E-01 
7 0.482192E+00 -0.497395E+00 -0.428618E-01 
8 0.577701E+00 -0.806931E+00 -0.575020E-01 
9 0.627437E+00 -0.126845E+01 -0.822389E-01 
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| OPTIMIZATION #4 


OPTIMIZATION SOLUTION 


A) Problem Parameters: 
3 


Arch Angle 0.000 Youngs Modulus: 30000000.0 
Arch Radius: 32,000 Yield Strength: 52000.0 
No of Design Var: 18 No of Elements: 8 
| B) Derived Constants: 
No of System Nodal Points... 9 
: No of Degrees of Freedom.... 27 
Length per Element.......... G2. 3:k 
Number of Iterations........ 2 
C) Structure Loading: 
6k SSS earner 0.0000 
00 9 6 eee ee ee ae -2000.0000 
5 6 Acerca 1000.0000 : 
5 6 ieee cern 0.0000 = 
D) Elemental Dimensions and Stress Distribution: 
Node Height Base Length Area 
i 4.20406 0.42041 Gree 7 310 1.76741 
Zz 4.17639 0.41764 6.27310 1.74422 
3 4.09209 0.40921 Gaze 310 1.67452 
4 3.94749 0.39475 6.27310 aooG2 7 
5 3.73436 0.37344 627310 e455 
6 3943729 0.34373 627310 feel 47 
7 3..02033 0730203 Gaz? 310 0.91224 
8 ie) ee 0.35108 6427310 0.68592 
9 0.79249 0.21196 6.27310 0.16798 
E) Objective Function: 
Total structure Volume: 63.252686 
Node Normal Stress Bending Stress Total 
1 li ZiGez 5$0868.9 51995..1 
2 Pio 50874.8 1994.1 
3 1098.2 5$0896.4 51994.6 
4 1062.0 $0929.1 51991 .1 
5 2009%3 50986.3 51395..6 
6 936.0 B05 5.6 51991.6 
7 835.0 1156.6 51991.8 
8 566.1 $1424.8 51990.9 
9 FICC. > 45071.7 46238.2 
F) Boundary Conditions: 
Node X-Displ Y-Displ Slope 
ah 1 iL ih 
9 0 0 0 
G) Solution Vector: 
Node X-Displ Y-Displ Slope 
1 0.000000E+00 0.000000E+00 0.000000E+00 
2 0.158777E-01 -0.180201E-02 -0.507744E-02 
3 0.619007E-01 -0.160076E-01 -0.102264E-01 
4 0.133298E+00 -0.544288E-01 -0.155283E-01 
5S 0.222312E+00 -0.127764E+00 -0.210892E-01 
6 0.318467E+00 -0.245256E+00 -0.270698E-01 
7 O.408886E+00 -0.414819E+00 -0.337708E-01 
8 0.479103E+00 -0.646809E+00 -0.425065E-01 
9 Q0.512886E+00 -0.990889E+00 -0.605568E-01 


ow 


OPTIMIZATION #4a 


OPTIMIZATION SOLUTION 


A) Problem Parameters: 
Arch Angle 90.000 Youngs Modulus: 30000000.0 
Arch Radius: 32.000 Yield Strength: 52000.0 
No of Design Var: 18 No of Elements: 8 
B) Derived Constants: 
No of System Nodal Points... 9 
No of Degrees of Freedom.... Za 
Length per Blement.......... 622731 
Number of Iterations........ 2 
C) Structure Loading: 
BX 6 sca ee: 0 ie 0 cuca eemmeaetete de: emeae evens 0.0000 
BY. 6 fee 5 0-4 elias eter n -2000.0000 
sy rarer ror i reo eld chee 10000.0000 Bs 
FA. 6 0 elses << iss sss) s) etemee ort enemene 0.0000 - 
D) Elemental Dimensions and Stress Distribution: 
Node Height Base Length Area 
M 3.99646 0.39965 6727310 1 S977 
2 3.96579 0.39658 6.27310 LeSiz5 
3 3.67202 0.38720 6.27310 1.49925 
4 3.70940 0.37094 6.27310 1.37596 
5 3.46539 0.34654 6.27310 1.20089 
6 3 pee 3 8 0.31124 6.27310 0.96869 
7 2.57622 0225763 6.23310 Osos 71 
8 1.43829 0.14383 Oe 310 0.20687 
9 2.259350 0.22769 6.27510 0.51444 
E) Objective Function: 
Total structure Volume: 53.206200 
Node Normal Stress Bending Stress Total 
ih 1246.2 5079772 52003.4 
2 1241.2 507617, 6 52002.8 
3 1226.5 50776. 2 52002.8 
4 1202.8 50799 .4 52002.1 
5 Pi 220 50828 .7 52000.7 
6 114156 50859.8 52001.4 
7 1147.6 50851.9 51999.5 
8 Ley 7a SOPz7 20 52004.0 
9 361770 Sl6Z2Ze2 5 ZO sret 
F) Boundary Conditions: 
Node X-Displ Y¥-Displ Slope 
1 1 J ik 
9 0 0 0 
G) Solution Vector: 
Node X-Displ Y-Displ Slope 
1 0.000000E+00 0.000000E+00 0.000000E+00 
2 0.166829E-01 -0.190698E-02 -0.533238E-02 
3 0.650796E-01 -0.168603E-01 -0.107525E-01 
4 0.140297E8+00 -0.573S560E-01 -0.163633E-01 
S 0.234385E+00 -0.134897E+00 -0.223074E-01 
6 0.336635E+00 -0.259876E+00 -0.288281E-01 
7 0.434064E+00 -0.442670E+00 -0.365050E-01 
8 0.515013E+00 -0.710561E+00 -0.496258E-01 
9 0.544307E+00 -0.100920E+01 -0.415951E-01 
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OPTIMIZATION #5 


OPTIMIZATION SOLUTION 


Arch Angle 90.000 Youngs Modulus: 30000000.0 
Arch Radius: 82-000 Yield Strength: 52000.0 
No of Design Var: 18 No of Elements: 8 
B) Derived Constants: 
No of System Nodal Points... 9 
No of Degrees of Freedon.... 754) 
Length per Element.......... 6.24 31 
Number of Iterations........ 2 
C) Structure Loading: 
000 60 0 eee 0.0000 
0 66 6 Sa 0.0000 
4 4 6 6a er 0.0000 a 
oo 6 Ae -100.0000 i 
D) Elemental Dimensions and Stress Distribution: 
Node Height Base Length Area 
jh 4.94233 0.49484 6927310 2.44566 
2 4.59883 0.46036 6.273510 2. tia 
3 4.20711 0.42104 6.27450 1.7 7b 
4 3.76635 0.37698 652750 1.42060 
5 3.20609 0.32787 6227.20 1.07416 
6 2.72509 Oi 2aezsd 6.27510 0.74152 
7 2.06293 020953 6227 210 0.43644 
8 1.31418 Ovis2s5 6.27310 0.17394 
9 0.14293 0.03000 62273 rh 0.00429 
E) Objective Function: 
Total structure Volume: §5.704273 
Node Normal Stress Bending Stress Total 
a lij3.3 $0828.3 5 290252 
2 1210.8 5079126 52002.6 
3 1109%3 5063s-1 §2002-9 
4 996.4 51006.4 SZ200247 
5 868.5 S1135.7 52004.2 
6 7124,.0 512797 5200357 
7 555.5 51446.8 S 2002-03 
8 35176 51646.)5 51998.3 
9 1.8 4.8 60 
F) Boundary Conditions: 
Node X-Displ Y-Displ Slope 
i 1 sh 1 
9 0 0 0 
G) Solution Vector: 
Node X-Displ Y-Displ Slope 
1 0.000000E+00 0.000000E+00 0.000000E+00 
2 0O.144353E-01 -0.168651E-02 -0.447172E-02 
3 0.564283£-01 -0.146790E-01 -0.932896E-02 
4 0.1236576+00 -0.508635E-01 -0.147199E-01 
S 0.2108486+00 -0.122671E+00 -0.208707E-01 
6 0.309650E+00 -0.243321E+00 -0.281898E-01 
7 0.408546E+00 -0.428622E+00 -0.375918E-01 
8 0.492877E+00 -0.706927E+00 -0.522128E-01 
9 0.557215E+00 -0.136017E+01 -0.130850E+00 
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OPTIMIZATION #6 


OPTIMIZATION SOLUTION 


A) Problem Parameters: 
Arch Angle 90.000 Youngs Modulus: 30000000.0 
Arch Radius: 32.000 Yield Strength: 52000.0 
No of Design Var: 18 No of Elements: 8 
B) Derived Constants: 
No of System Nodal Points... 9 
No of Degrees of Freedom.... rae | 
Length per Element.......... 6.2731 
Number of ~fterations........ i 
C) Structure Loading: 
FXG: 5 so. wile: evtey Rear Cth ne ele eo ners 0.0000 
BOL so 5 Sou ene sues one c cee a retoneuens -2000.0000 
PM esc Siue wks orehe tans terete ee esmeua omen ae 0.0000 b 
1 \ Pree ere erry yaa eo -100.0000 - 
D) Elemental Dimensions and Stress Distribution: 
Node Height Base Length Area 
1 5 e2s26 0.58232 62273510 3.39091 
2 5.56534 0.551653 6.27310 3.09731 
3 5 25560 0.52550 6.27310 2.16025 
4 4.88828 0.48883 67273510 2.36953 
= 4.46201 0.44620 62273510 1.99096 
6 3.96199 0.39620 627310 1.56974 
7 3.35911 0233591 6.27310 1 . D286 
8 2.27093 0.32719 6.273510 0.74303 
9 0.40176 0.51850 6.27310 0.20631 
E) Objective Function: 
Total structure Volume: 97.474487 
Node Normal Stress Bending Stress Total 
at 1433.6 50561.4 31995... 
?, 1457.9 50537 .3 51995.. 2 
3 137854 50617.0 513995.-5 
4 1284.9 $0710.93 51995.8 
= 117538 50821.4 91996.9 
6 1046.4 §0951.4 51997.8 
ii 889.9 S1 10972 51999.1 
8 605.0 513937.6 51998 .6 
9 940.6 0.6 941.2 
F) Boundary Conditions: 
Node X-Displ Y-Displ Slope 
iL ih 1 1 
9 0 0 0 
G) Solution Vector: 
Node X-Displ Y-Displ Slope 
1 0.000000E+00 0.000000E+00 0.000000E+00 
2 0.118390E-01 -0.148107E-02 -0.371830E-02 
3 0.461727E-01 -0.122074E-01 -0.763781E-02 
4 0.100311E+00 -0.414621E-01 -0.118319E-01 
5S 0.169098E+00 -0.982482E-01 -0.164007E-01 
6 0.244936E+00 -0.191026E+00 -0.215055E-01 
7 QO.317850E+00 -0.327870E+00 -0.274519E-01 
8 0.375826E+00 -0.519533E+00 -0.353806E-01 
9 0.404874E+00 -0.815207E+00 -0.533505E-01 
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OPTIMIZATION #6a 


OPTIMIZATION SOLUTION 


} 


— naan 


er 


A) Problem Parameters: 
Arch Angle 90.000 Youngs Modulus: 30000000.0 
Arch Radius: 32-600 Yield Strength: 52000.0 
No of Design Var: 18 No of Elements: 8 
B) Derived Constants: 
No of System Nodal Points... 9 
No of Degrees of Freedon.... 74] 
Length per Element.......... Soest 
Number of Iterations........ 1 
C) Structure Loading: 
ss 0.0000 
6 6 a -2000.0000 
"5 3 SSS 0.0000 7 
0 esa -100.0000 : 
D) Elemental Dimensions and Stress Distribution: 
Node Height Base Length Area 
1 5.82323 OFoo232 6427310 3.239201 
2 5.56545 0.55654 6.27310 3.09742 
3 §.25394 0552539 6.27310 2.76039 
4 4.88843 0.48884 623 10 2.38968 
5 4.46221 0.44622 6.27310 1.99114 
6 3.96197 029623 6527310 1.56986 
7 2.96078 0.43150 6427310 Laz s9 
8 2. Oiled 7 i 0.41514 6227310 0.83639 
9 O7Ss831 iwa0792 6227310 0.65023 
eyeorpyective Function: 
Total structure Volume: 101.764938 
Node Normal Stress Bending Stress Total 
L 143306 50 Si516..2 51989.7 
vi 1457925 50531.4 51989.2 
3 131 613 50609.9 51988.2 
4 1284.9 50703 .3 S198e..2 
5 11753 50812.4 51 987'28 
6 1046.3 5$0946.8 $1993.21 
7 TBOae Sl 21 eee 51997.3 
8 53720 51463.1 SZ000m6 
9 301.4 O=5 302.9 
F) Boundary Conditions: 
Node X-Displ Y-Displ Slope 
1 Hi it 1 
9 0 0 0 
G) Solution Vector: 
Node X-Displ Y-Displ Slope 
1 O.Q000000E+00 O0.000000E+00 0.000000E+00 
2 0.118375E-01 -0.148090E-02 -0.371783E-02 
3 0.461666E-01 -0.122058E-01 -0.763674E-02 
4 0.100297E+00 -0.414562E-01 -0.118301E-01 
5S 0.169072E+00 -0.982332E-01 -0.163980E-01 
6 0.244899E+00 -0.190995E+00 -0.215019E-01 
7 0O.318251E+00 -0.328633E+00 -0.277224E-01 
8 0.377833E+00 -0.525529E+00 -0.367195E-01 
9 0.404775E+00 -0.799478E+00 -0.474619E-01 


dS) 


OPTIMIZATION #7 


OPTIMIZATION SOLUTION 


Arch Angle 90.000 Youngs Modulus: 30000000.0 
Arch Radius: a2,.000 Yield Strength: 52000.0 
No of Design Var: iS No of Elements: 8 
B) Derived Constants: 
No of System Nodal Points... 9 
No of Degrees of Freedon.... Pa 
Length per Element.......... 6.2731 
Number of Iterations........ 1 
C) Structure Loading: 
Bere “e's 6 -actdiws le eee nots ae nenememeneetetteme pt omremes 0.0000 
BY «os wp Wak acs ene ote eat. omen mem aes -8000.0000 
FM < «1s aleve. 69 ouemessuenenets seem m ene memenen 0.0000 ~ 
BA cjichs «0 5 3» 0 eee ts ete 3 0) eee 0.0000 - 
D) Elemental Dimensions and Stress Distribution: 
Node Height Base Length Area 
1 0.49766 O32385 6227210 0.16.7 
Z 2.49460 0.47039 6.27310 1.17348 
3 3.05960 0.46175 6.20320 1.41278 
4 2.96505 0.48983 6.273510 1245238 
5 2.41709 0.49892 6.27310 1.20593 
6 0.96211 0.15996 6.273510 0.15390 
7 3.08575 0.43109 6.27510 1.33023 
8 4.00244 0.54499 6.27310 2.16329 
9 4.21890 0.79543 6.27310 3.35565 
E) Objective Function: 
Total structure Volume: 64550076 
Node Normal Stress Bending Stress Total 
i S199 86 0.2 S199 eg 
2 1364066 44637.4 59 Siz 
3 6359. 0 45640.2 S1999-3 
4 61657. 7 4$613..3 5199920 
5 TiGsee7 44836.1 5199927 
6 51734.2 267.2 52001.4 
7 $246.1 46199.3 $1445.4 
8 2625.0 49374.5 5 2000s 
9 1501.8 50497.5 5$1999.4 
F) Boundary Conditions: 
Node X-Displ Y-Displ Slope 
i i) i 0 
9 i) 0 i) 
G) Solution Vector: 
Node X-Displ Y-Displ Slope 
1 0.000000E+00 0.000000E+00 0.307764E-01 
2 -0.149671E+00 0.117776E-01 0.102322E-01 
3 -0.192288E+00 0.232035E-01 0.336248E-02 
4 -0.193966E+00 0.225997E-01 -0.298150E-02 
S -0.161968E+00 -0.546921E-02 -0.100998E-01 
6 -0.785744E-01 -0.112048E+00 -0.273030E-01 
7 -0.16108S5E-01 -0.234507E+00 -0.111216E-01 
8 -0.133252E-02 -0.285889E+00 -0.514754E-02 
9 0.000000E+00 -0.303321E+00 0.000000E+00 
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OPTIMIZATION #7a 


OPTIMIZATION SOLUTION 


ow a a @ oe oO oO oe Oo oe oe Oo ow oe oe oe © 6 6 oe oe oe we 6 ee 6 ee oe oe oe ee oe oe oe ee ee oe oe oe ee ee oe oe Oo oe oe 6 ee ee ee ow ow oe ee oe ow oe oe oe ow oe oe oe oo oo = 


Arch Angle 90.000 Youngs Modulus: 30000000.0 
Arch Radius: 32.000 Yield Strength: a 0. 0 
No of Design Var: 18 No of Elements: 8 
B) Derived Constants: 
No of System Nodal Points... 9 
No of Degrees of Freedonm.... 2] 
Length per Element.......... 6227 31 
Number of Iterations........ 2 
C) Structure Loading: 
PES os 5 5 oe Se es ce ee ee 0.0000 
Son 65) rere -8000.0000 
so 6 5 ARIS Ae 0.0000 ~ 
oO 5 6 6 ee 0.0000 - 
D) Elemental Dimensions and Stress Distribution: 
Node Height Base Length Area 
af 0267.57 0 0.23920 6.27310 0216115 
Z 2.61798 0.42995 6a 7510 lez oo 1 
3 3.27094 0.40726 Geez 75 10 1 Siz 2 
4 Sel 5266 0.43628 6.27310 137546 
5 2 9 31 052326 6.27310 lezen 
6 0.77194 Onl9791 6.27310 0.15277 
q| 30g 183 2 0.44320 6227510 1.34348 
8 4.32486 0.46885 Gaz 310 2 Oar Z 
9 4.62704 Ue 55.0 6.27310 3.069450 
E) Objective Function: 
Total structure Volume: 61.674786 


Node Normal Stress Bending Stress Total 
af 20056 Or 52005..8 
2 7674.6 44327.7 520023 
3 6743.9 45258 .6 52002235 
4 6531.34 45470.1 §2001.4 
5 07 3.5 44927.7 3200 la 
6 Size 10 wl 11218 .4 532365 
7 5193 a7 46613.6 5180723 
8 2824.0 4917733 320 Gli 
9 1641.8 50360.2 5200271 
F) Boundary Conditions: 
Node X-Displ Y-Displ Slope 
sh 1 1 0 
9 i 0 yh 
G) Solution Vector: 
Node X-Displ Y-Displ Slope 


0.000000E+00 
-~-0.140675E+00 


0.000000E+00 
0.106575E-01 
-0.184847E+00 0.224797E-01 
-0.190542E+00 0.239361E-01 


ih 0.285876E-01 
2 

3 

4 

5 -0.163888E+00 0.229189E-03 

6 

7 

8 

9 


0.102757E-01 
0.386082E-02 
-0.204050E-02 
~0.886111E-02 


~0.790007E-01 .108136E+00 -0.284664E-01 


~0.148789E-01 ~-0.233588E+00 -0.104461E-01 
~0.115786E-02 -0.281574E+00 -0.472253E-02 
0.000000E+00 -0.297573E+00 0.000000E+00 
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OPTIMIZATION #8 


OPTIMIZATION SOLUTION 


Arch Angle 90.000 Youngs Modulus: 30000000.0 
Arch Radius: 32.000 Yield Strength: 52000.0 
No of Design Var: 18 No of Elements: 8 
B) Derived Constants: 
No of System Nodal Points... 9 
No of Degrees of Freedom.... 27 
Length pér Element...:...... 6.2731 
Number of Iterations........ Z 
C) Structure Loading: 
2) RP PE ea so Se yee ee 0.0000 
¢) QE are a eras coo 0 cee -8000.0000 
BAM ss. 6 ss vei: oleic 61s ore rer eure ememamemen tis 0.0000 4 
PAs ce © eielele wie eee sseneneieeeset se tenemos 0.0000 - 
D) Elemental Dimensions and Stress Distribution: 
Node Height Base Length Area 
ak 3.96647 0.39665 6.27310 1.59329 
2 2.66054 0.26605 6.27310 0.70784 
3 2.69799 0.26980 6427310 O.7 2791 
4 3.29169 0.32917 6.27310 1.08352 
5 3229261 0.32928 622310 1.08426 
6 2.70298 0.27030 6.27310 0.73061 
7 2.65240 0.26524 6.27310 0270352 
8 3.95012 0.391511 65277510 1.56114 
9 4.87425 0.48694 6.27310 2.31347 
E) Objective Function: 
Total structure Volume: 52.992058 
Node Normal Stress Bending Stress Total 
i 5469.7 46526.9 51996.6 
Z 12833.6 39153.1 51986.7 
3 13542.6 38456.9 51999".5 
4 9462.4 42532.7 51995..1 
5 9456.6 42753 3..2 51989.9 
6 13496.0 38480.4 51976.4 
7 L129 Tbe 39108.3 52026.4 
8 S103 2 46906.0 52009.2 
9 3085., 5 48816.0 SL 90.5 
F) Boundary Conditions: 
Node X-Displ Y-Displ Slope 
i il ak 1 
9 sl 0 ne 
G) Solution Vector: 
Node X-Displ Y-Displ Slope 
1 0.000000E+00 0.000000E+00 0.000000E+00 
2 -0.237838E-01 0.695637E-03 0.631554E-02 
3 -0.683462E-01 0.113025E-01 0.624288E-02 
4 -0.894876E-01 0.199189E-01 0.336519E-03 
5 -0.792890E-01 0.896548E-02 -0.506634E-02 
6 -0.479138E-01 -0.329894E-01 -0.109648E-01 
7 -0.147824E-01 -0.100892E+00 -0.110874E-01 
8 -0.821552E-03 -0.152604E+00 -0.471284E-02 
9 0.000000E+00 -0.168974E+00 0.000000E+00 
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| OPTIMIZATION #8a 


OPTIMIZATION SOLUTION 


Arch Angle 90.000 Youngs Modulus: 30000000.0 
Arch Radius: 32.000 Yield Strength: §2000.0 
Arch Height: 2.000 No of Elements: 12 
B) Derived Constants: 
No of System Nodal Points... 13 
No of Degrees of Freedon.... 39 
Length per Element.......... 4.1858 
Number of Iterations........ 
C) Structure Loading: 
og 5 0.0000 
0 oS -8000.0000 
5 6 RS Cee O20 000 
3 g Ser 0.0000 : 
D) Elemental Dimensions and Stress Distribution: 
Element Height Base Length Volume 
sl 2.00000 LotZao7 4.18580 14.45838 
2 2.00000 0.94577 4.18580 1egl62 
3 2.00000 0.27194 4.18580 2.27661 
4 2.00000 0.65040 4.18580 5.44488 
5 2.00000 1.09075 4.18580 9.133 0 
6 2.00000 0.78986 4.18580 6.61256 
zi 2.00000 Ve0622 3 4.18580 8.89255 
8 2.00000 0.75368 4.18580 6-3 1116 
9 2.00000 0.25849 4.18580 226 4100 
10 2.00000 On9327 5 4.18580 7.80861 
hat 2.00000 e792 03 4.18580 14.67236 
Wd 2.00000 2.65143 4.18580 Pap ees oie Ea) 
E) Objective Function: 
Total structure Volume: 107.886574 
Node Stress 
1 S991). 36 
2 49986.21 
3 46231.32 
4 51991.33 
=) 49404.85 
6 51991.04 
uu 51841.79 
8 42080.21 
9 $1997.34 
10 5279.07 
lz $2434.15 
eZ 51908.94 
13 54,935.26 
F) Boundary Conditions: 
Node X-Displ Y-Displ Slope 
q 1 af 
1s - 0 1 


Weg 


G) 


Solution Vector: 


Node 


0. 
=n 
20m 
-0. 
moe 
-0. 
.110454E+00 
.850803E-01 


-0 
-0 


-0. 
-0. 
-0. 
=Or 

0. 


X-Displ 

OO0O0000E+00 
121854E-01 
442765E-01 
861533E-01 
112986E+00 
121058E+00 


$37900E-01 
252856E-01 
744947E-02 
654979E-03 
000000E+00 


Y-Displ 
0.000000E+00 
0.457452E-03 
0.615066E-02 
0.177060E-01 
0.297132E-01 
0.342979E-01 
0.237526E-02 

-0.622361E-02 
-0.547404E-01 
-0.118485E+00 
-0.173190E+00 
-0.209063E+00 
-0.222069E+00 


ube (C, 


Uh 8 to Ua 
OoOooo0o0oeo0o09c0ec°oeo 


Slope 


-000000F+00 
.527930E-02 
-947556E-02 
.907372E-02 
-435047E-02 
.109393E-03 
-701848E-02 
~-115516E-01 
-LS5SOQ000E-01 
~155S5S01E-01 
~-110976E-01 
-580765E-02 
-O00000E+00 


OPTIMIZATION #9 


OPTIMIZATION SOLUTION 


A) Problem Parameters: 


Arch Angle 180.000 Youngs Modulus: 30000000.0 
Arch Radius: 32.000 Yield Strength: 52000.0 
No of Design Var: 18 No of Elements: 8 
B) Derived Constants: 
No of System Nodal Points... 9 
No of Degrees of Freedon.... 27 
Gength per Element.......... 12.4858 
Number of Iterations........ 2 
C) Structure Loading: 
5 6 6 ISIS AS Sareea 16000.0000 
MEI 5 ct et tt eee eee 0.0000 
err 0.0000 . 
a EE cn 0.0000 - 
D) Elemental Dimensions and Stress Distribution: 
Node Height Base Length Area 
af 3387225 2.6 18130 12.48578 10723610 
2 2.99788 1.94284 12.48578 5.82439 
3 0.78829 0732181 12.48578 0225566 
4 2.77540 1.34286 12.48578 3572096 
5 Zeoze9l 10s S748 abe 12.48578 3274337 
6 2.47737 eOlce 2 1 12.48578 2.54492 
7 P6017 2 1.17975 12.48578 1.88963 
8 0.98586 0.84212 12.48578 0.83021 
9 Deel 3 2 1.14827 12.48578 1.46670 
E) Objective Function: 
Total structure Volume: 287.147583 
Node Normal Stress Bending Stress Total 
Jt 5.1 35 51426.4 §1999.9 
2 1458.0 50542.0 §2000-0 
3 $1241.47 770.7 §2012.4 
4 4166.3 47833.6 §2000.0 
5 22392 49761.4 §2000.5 
6 404.8 §1595..6 5200055 
ri 100771 S0Ssi3 51998.4 
8 2995.1 49002.5 51997 .6 
9 183522 0.1 16S 5 2 
F) Boundary Conditions: 
Node X-Displ Y-Displ Slope 
ue ie 3 1 
9 0 1 0 
G) Solution Vector: 
Node X-Displ Y-Displ Slope 
1 0.000000E+00 0.000000E+00 0.000000E+00 
2 0.904933E-01 -0.176853E-01 -0.130647E-01 
3 0.558262E+00 -0.327630E+00 -0.608843E-01 
4 0.883004E+00 -0.806151E+00 -0.171659E-01 
5 0.909118E+00 -0.928170E+00 -0.266783E-02 
6 0.923784E+00 -0.854074E+00 0.134924E-01 
7 0.109988E+01 -0.590018E+00 0.348143E-01 
8 0.169034E+01 -0.194612E+00 0.717773E-01 
9 0.266363E+01 0.000000E+00 0.833529E-01 


ey. 


OPTIMIZATION #9a 


OPTIMIZATION SOLUTION 


A) Problem Parameters 


Arch Angle 180.000 Youngs Modulus: 
Arch Radius: 32.000 Yield Strength: 
No of Design Var: 18 No of Elements: 
B) Derived Constants: 
No of System Nodal Points... 9 
No of Degrees of Freedonm.... 27 
Length per Element.......... 12.4858 
Number of Iterations........ 2 
C) Structure Loading: 
1) CParererer owe Rea Gra. Gunso cutis ic 16000.0000 
sp PEE er re eae oA eG ee 0.0000 
OPM sos shaiecelgWeleteleds ie acelsm ee mem eiemt Ee 0.0000 
DAs Sue cuseerer: ME : 0.0000 
D) Elemental Dimensions and Stress Distribution: 
Node Height Base Length 
1 lL u0rs55 0.19099 12.48578 
2 5.49600 0.69864 12.48578 
3 5.93999 0.96376 12.48578 
4 5 Sela 7 1.04035 12.48578 
5 6.00000 0.54928 WZ .46578 
6 5.573 0.58244 12.48578 
7 4.24805 0.49093 12.48578 
8 257 0.52420 12.48578 
9 0.24362 0.66203 12.48578 
E) Objective Function: 
Total structure Volume: 344.338989 
Node Normal Stress Bending Stress 
1 53383 .0 ees 
2 3451.4 S0lcie! 
3 2907.2 50644.1 
4 28 L250 50764.5 
5 2854.3 77667.4 
6 9257 52414.6 
7 200072 50774.2 
8 6355.6 4711726 
9 48642.3 0.0 
F) Boundary Conditions: 
Node X-Displ Y-Displ Slope 
1 1 ih 0 
9 0 i 0 
G) Solution Vector: 
Node X-Displ Y¥-Displ Slope 
1 0.Q000000E+00 0.000000E+00 -0.532524E-01 
2 0.538764E+00 -0.103983E+00 -0.253326E-01 
3 0.766928E+00 -0.254801E+00 -0.178851E-01 
4 0.867834E+00 -0.403592E+00 -0.108023E-01 
5S 0.885116E+00 -0.482743E+00 -0.246003E-02 
6 0.891319E+00 -0.450538E+00 0.696965E-02 
7 0.974611E+00 -0.324620E+00 0.161255E-01 
8 0.122814E£+01 -0.153178E+00 0.301513E-01 
9 0.197489E+01 0.000000E+00 0.765027E-01 


30000000.0 
52000.0 


Tota 


sh 


53384. 
336i. 
535517 
535708 
80521. 
53333 
53434. 
53473. 
48642. 


OFrPNWWHUW © 


eZ 


WN & J Uw UW 


8 


OPTIMIZATION #10 


OPTIMIZATION SOLUTION 


A) Problem Parameters: 


Arch Angle 180.000 Youngs Modulus: 30000000.0 
Arch Radius: 32.9000 Yield Strength: 52000.0 
No of Design Var: 18 No of Elements: 8 
B) Derived Constants: 
No of System Nodal Points... 9 
No of Degrees of Freedom.... 27 
Length per Element.......... 12.4858 
Number of Iterations........ 2 
C) Structure Loading: 
Oa ae 16000.0000 
Os 6 a 0.0000 
1 55 SS ee 0.0000 - 
6 6 oa 0.0000 - 
D) Elemental Dimensions and Stress Distribution: 
Node Height Base Length Area 
ah 3.87807 e257 526 12.48578 9.98702 
2 3.14748 D269 3 12.48578 5.43547 
3 0.88930 0729131 12.48578 0.25906 
4 2.92467 114525 12.48578 3.34919 
5 296633 1705599 12.48578 3e 13335 
6 220310 1 03918 12.48578 2.10318 
a 0.77820 0.60572 12.48578 0.47137 
8 1.64419 0.87326 12.48578 1.43580 
9 Wes 70 1 1.15400 iZeaoo 76 1.94680 
E) Objective Function: 
Total structure Volume: 256.608276 
Node Normal Stress Bending Stress Total 
ii 63177 51368.8 52000..5 
2 1646.5 $0354 .2 52000.7 
3 5153056 479.7 52010n3 
4 4692.9 47307 .6 5200025 
S 2706.4 49294.1 52000.5 
6 5798 §1421.,7 52001. > 
7 4780.6 47222.8 52003.4 
8 2050.8 49951.2 5200220 
9 l637-2 50364.7 52001.8 
F) Boundary Conditions: 
Node X-Displ Y-Displ Slope 
1 ah 1 ar 
9 0 i! 
G) Solution Vector: 
Node X-Displ Y-Displ Slope 
1 0.000000E+00 0.000000E+00 0.000000E+00 
2 0.882631E-01 -0.172024E-01 -0.127257E-01 
3 0.520760E+00 -0.303339E+00 -0.559400E-01 
4 0.823131E+00 -0.747606E+00 -0.172292E-01 
5 0.850524E+00 -0.874578E+00 -0.359896E-02 
6 0.865174E+00 -0.800410E+00 0.136811E-01 
7 QO.115273E+01 -0.368882E+00 0.570417E-01 
8 0.165293E+01 -0.331472E-01 0.251769E-01 
9 0.181555E+01 0.000000E+00 0.000000E+00 


PS 


OPTIMIZATION #11 


OPTIMIZATION SOLUTION 


A) Problem Parameters: 


Arch Angle 180.000 Youngs Modulus: 30000000.0 
Arch Radius: S2,000 Yield Strength: §2000.0 
No of Design Var: 18 No of Elements: 8 
B) Derived Constants: 
No of System Nodal Points... 9 
No of Degrees of Freedon.... y Hae 
Length per Element.......... 12.4858 
Number of Iterations........ 
C) Structure Loading: 
Pre ERG 5c do Oe 0.0000 
FA so oo: o'slave se Siemeuery eaeeteteneueeene sists -12000.0000 
i Pr rrr or es go So 1000.0000 Z 
2, rere ere ieee 4 0.0000 : 
D) Elemental Dimensions and Stress Distribution: 
Node Height Base Length Area 
aL 1.15824 IV02z2565 12.48578 1.19142 
Z 2.00490 iz 219 12.48578 2.249387 
3 2.04000 iO e818 22.465/8 2.21989 
4 0.57320 0.21096 12.48578 0.12092 
5 2.42160 1.36445 12.48578 3.30425 
6 1.23847 Ce 7 On22 12.48578 0.87587 
y 1.54335 1.08778 12.48578 1.67882 
8 1.43388 0.95977 12.48578 1.37619 
9 1.49848 1.26079 12.48578 1.88926 
E) Objective Function: 
Total structure Volume: 153.076752 
Node Normal Stress Bending Stress Total 
iL 5944.7 0.0 5944.7 
iz 3270 48720.7 51998.4 
3 33zZ179 48527.9 51849 .8 
4 51700me 2535.5 5195324 
5 1618.5 $0381.52 51999.8 
6 681720 45180.4 §1997.4 
7 4083.4 47917.5 52000.9 
8 4865.8 47133.8 51999 .6 
9 3360.4 48642.0 5200223 
F) Boundary Conditions: 
Node X-Displ Y-Displ Slope 
ak alk al 0 
9 a at ah 
G) Solution Vector: 
Node X-Displ Y-Displ Slope 
1 0.000000E+00 0.000000E+00 0.456664E-01 
2 -0.471806E+00 0.920810E-01 0.241675E-01 
3 -0.619634E+00 0.189140E+00 0.415893E-02 
4 -0.359468E+00 -0.206478E+00 -0.589599E-01 
5S -0.269622E+00 -0.667966E+00 0.495288E-02 
6 -0.213488E+00 -0.379826E+00 0.343967E-01 
7 0.128494E-01 -0.370866E-01 0.214467E-01 
8 0.888700E-01 0.160224E-01 -0.531602E-02 
9 0.000000E+00 0.000000E+00 0.000000E+00 
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OPTIMIZATION #lla 


OPTIMIZATION SOLUTION 


- oe oe oe ee oe ee ae ee es i i a ss a es a es i a i a es a es es es es ee a a es es a a 


A) Problem Parameters: 


Arch Angle : 180.000 Youngs Modulus: 30000000.0 

Arch Radius: 32.000 Yreld Strength: 200070 

No of Design Var: 18 No of Elements: 8 
B) Derived Constants: 

No of System Nodal Points... 9 

No of Degrees of Freedom.... ai 

Length per Element.......... 12.4858 

Number of Iterations........ Z 


C) Structure Loading: 


6 o 6 eae 0.0000 

4 5 -24000.0000 

7. 4 ivoke momen tenes: 1000.0000 _ 

SS 0.0000 7 

D) Elemental Dimensions and Stress Distribution: 

Node Height Base Length Area 
1 1.25459 eee 70 12.48578 1.39472 
Z 2.56424 1.39293 12.48578 3.57180 
3 2.61780 1.33869 12.48578 3.50443 
4 0.76364 Om3i775 12.48578 0.24265 
5 3.03972 1.74821 02246578 5231405 
6 1.59975 0.89444 12.48578 1.43088 
7 1.98431 1.34883 12.48578 2.67651 
8 Wn GS 7S i8 1.19860 12.48578 Fins PA PAU EO) 
9 1.96692 1.46917 12.48578 2.88973 


E) Objective Function: 


Total structure Volume: 241.778809 


Node Normal Stress Bending Stress Total 
ih Or-3 Sir6 0-0 O23 Siw 
2 4120.0 47880.2 5200052 
3 4199.2 47797.9 51997 1 
4 5 471 .3r 5 otis $2025,.0 
5 2009.1 49990.9 5200070 
6 8340.7 48657.1 51997.9 
7 $124.9 46875.8 52000.8 
8 6039.8 45960.0 51999.8 
9 4402.2 47599.3 5 2.0 0.5 

F) Boundary Conditions: 

Node X-Displ Y-Displ Slope 
af i i 0 
9 1 1 1 


Gy solution Vector: 
Node X-Displ Y-Displ Slope 
0.Q000000E+00 0.000000E+00 0.377959E-01 
-0.377866E+00 0.726540E-01 0.168604E-01 
-0.474151E+00 0.134826E+00 0.149644E-02 
~0.277207E+00 -0.167482E+00 -0.440847E-01 
~0.210637E+00 -0.513900E+00 0.340668E-02 
-0.168729E+00 -0.295879E+00 0.262287E-01 
0.35986S5E-02 -0.330045E-01 0.166200E-01 
0.647848E-01 0.107673E-01 -0.365121E-02 
0.000000E+00 0.000000E+00 0.000000E+00 


w@mr~yInnh WwW NF 
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OPTIMIZATION #12 


OPTIMIZATION SOLUTION 


A) Problem Parameters: 


Arch Angle 180.000 Youngs Modulus: 30000000.0 
Arch Radius: 32.000 Yield Strength: 52000.0 
No of Design Var: 18 No of Elements: 8 
B) Derived Constants: 
No of System Nodal Points... 9 
No of Degrees of Freedom.... 27 
Length per Element....-..... 12.4858 
Number of Iterations........ 2 
C) Structure Loading: 
a er 5 or -9000.0000 
hE rE Hr GeO ere -17000.0000 
PSAP EPEC eres Fo Co Go 1000.0000 - 
7 Serre eer oe 0.0000 = 
D) Elemental Dimensions and Stress Distribution: 
Node Height Base Length Area 
ak 2.07490 1764551 12.48578 3.41427 
2 0.81408 0.49069 12.48578 0.39946 
3 2.06842 0.85360 12..46578 1.76560 
4 2.08744 We 85802 12.48578 1.78898 
> 1 o66 2 0. 1166 12.48578 0.12468 
6 2.82593 0.84646 12.48578 2.39204 
7 3.62188 1218265 12.48578 4.28341 
8 0.81741 Oma oo 3 1 12.48578 0.20869 
9 2.18518 1.96045 12.48578 4.28395 
E) Objective Function: 
Total structure Volume: 156.554611 
Node Normal Stress Bending Stress Total 
iL 133122 50668.2 §51999.4 
2 131592 38839.6 51998.8 
3 3525.9 48473.0 51998 .9 
4 3491.6 48507.9 §1999.4 
5 42640.2 6458.5 49098 .6 
6 1495.5 5050326 51999.1 
i 136622 $0631.9 5200022 
8 5193555 45.6 §1981.4 
9 2914.9 49085.0 520000 
F) Boundary Conditions: 
Node X-Displ Y-Displ Slope 
1 1 1 ah 
9 ah dh uh 
G) Solution Vector: 
Node X-Displ Y-Displ Slope 
1 0.000000E+00 0.000000E+00 0.000000E+00 
2 -0.385756E+00 0.754815E-01 0.480384E-01 
3 -0.777271E+00 0.334000E+00 0.140370E-01 
4 -0.808788E+00 0.378437E+00 -0.538832E-02 
5 -0.737499E+00 0.351551E-02 -0.430591E-01 
6 -0.812780E+00 -0.364503E+00 -0.423870E-02 
7 -0.800801E+00 -0.345997E+00 0.959411E-02 
8 -0.415861E+00 -0.859055E-01 0.510140E-01 
9 Q.000000E+00 0.000000E+00 0.000000E+00 
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OPTIMIZATION #13 


OPTIMIZATION SOLUTION 


A) Problem Parameters: 
Arch Angle 180.000 Youngs Modulus: 30000000.0 
Arch Radius: 3272000 Yield Strength: 52000.0 
No of Design Var: 18 No of Elements: 8 
B) Derived Constants: 
No of System Nodal Points... 9 
No of Degrees of Freedom.... 27 
Bengensper Element.......... 12.4858 
Number of Iterations........ Zz 
Coeseaue ture Loading: 
oo een 16000.0000 
- sR 0.0000 
Es tt ccc ccc ee 1000.0000 z 
Te sc st te eee 0.0000 : 
D) Elemental Dimensions and Stress Distribution: 
Node Height Base Length Area 
1 3.66535 2.61929 12.48578 10.18466 
iz 3.06710 1.86559 12.48578 5.18326 
3 0.92232 0.26590 12.48578 0.24524 
4 2.98689 1.21410 12.48578 3.62638 
5 3.00467 1.32694 12.48578 3.98700 
6 1.88382 0.97124 12548578 1.82964 
V 0.67362 0.55665 12.48578 0.37497 
8 1.49367 0.78025 12.48578 1.16544 
9 1.679697 1.18504 12.48578 1.99048 
E) Objective Function: 
Total structure Volume: 265.960205 
Node Normal Stress Bending Stress Total 
i 236.3 51465.6 §2001.9 
2 141223 $0589.8 5200220 
3 §1994.3 96.5 §2090.8 
4 4244.8 47 7S $2002 
5 3935.9 48065.6 5200-2 
6 4523.6 47478.1 52001.8 
7 4413.4 4757220 $1986.1 
8 18 552 $0144.8 §2000.0 
9 agli coy 7/ 3653845 3771142 
F) Boundary Conditions: 
Node X-Displ Y-Displ Slope 
1 ah i uh 
9 0 1 sh 
Gymsewucion Vector: 
Node X-Displ Y-Displ Slope 
1 0.000000E+00 0.000000E+00 0.000000E+00 
2 Q0.892425E-01 -0.174556E-01 -0.128954E-01 
3 0.527041E+00 -0.307447E+00 -0.566398E-01 
4 0.833359E+00 -0.758307E+00 -0.176816E-01 
5S 0.862231E+00 -0.894399E+00 -0.436349E-02 
6 0.880363E+00 -0.814806E+00 0.140516E-01 
7 0.117817E+01 -0.368085E+00 0.590398E-01 
8 0.168336E+01 -0.291622E-01 0.222087E-01 
9 0.182676E+01 0.000000E+00 0.000000E+00 


Ii 7) 


OPTIMIZATION #13a 


OPTIMIZATION SOLUTION 


A) Problem Parameters: 
Arch Angle 180.000 Youngs Modulus: 30000000.0 
Arch Radius: 32.000 Yield Strength: 52000.0 
No of Design Var: 18 No of Elements: 8 
B) Derived Constants: 
No of System Nodal Points... 9 
No of Degrees of Freedonm.... 27 
Length per Element.......... T2248 56 
Numbet of Iterations........ 2 
C)_ Structure Coadine: 
FN eae 5 ene os «20 cs 0; erent ee meee 16000.0000 
brn a re OC ett 6. 0.0000 
Pee aceon aside di coke os (ey teens 1000.0000 = 
of WEEE ra en A Sa Ga OG 0.0000 - 
D) Elemental Dimensions and Stress Distribution: 
Node Height Base Length Area 
i 5 . 58igion 1 2as7 12.48578 7.11864 
Z S.10ia6 0.69418 12.48578 3.54200 
3 1. S6ou7 0.1652 12.48578 0.24497 
4 4.38397 0.58589 12.48578 2.56854 
5 4.63376 0.58228 12.48578 2.69813 
6 3.27335 0.34568 12.48578 1 sissies 
7 1.61894 0.16189 12.48578 0.26210 
8 2.57960 0.25796 12.48578 0.66544 
9 2.91537 0.29154 12.48578 0.84994 
E) Objective Function: 
Total structure Volume: 175.647415 
Node Normal Stress Bending Stress Total 
1 169.0 5123471 51999.8 
2 2303.6 49696.1 51999.8 
3 52019 .3 211-4 52230.6 
4 5991.5 46008.2 51999.7 
5 S6l6.¢0 46183.7 52999..7 
6 130161 44698 .3 51999.4 
7 6284.3 35672.9 41957.2 
8 3234.0 48765.9 51999.9 
9 2740.5 47790.4 5057029 
F) Boundary Conditions: 
Node X-Displ Y-Displ Slope 
1 1 iL a 
9 0 1 i 
G) Solution Vector: 
Node X-Displ Y-Displ Slope 
1 0.000000E+00 0.000000E+00 0.000000E+00 
2 0.563288E-01 -0.107647E-01 -0.813259E-02 
3 0.306688E+00 -0.174214E+00 -0.318503E-01 
4 0.481161E+00 -0.425410E+00 -0.977707E-02 
5 0.497296E+00 -0.493441E+00 -0.126135E-02 
6 0.513184E+00 -0.431268E+00 0.955349E-02 
7 0.666860E+00 -0.199686E+00 0.294584E-01 
8 0.932019E+00 -0.202692E-01 0.147632E-01 
9 0.102734E+01 0.000000E+00 0.000000E+00 


Lvs 


OPTIMIZATION #14 


OPTIMIZATION SOLUTION 


A) Problem Parameters: 


Arch Angle 180.000 Youngs Modulus: 30000000.0 
Arch Radius: 320 00 Yield Strength: 52000.0 
No of Design Var: 18 No of Elements: 8 
B) Derived Constants: 
No of System Nodal Points... 9 
No of Degrees of Freedonm.... 77 
Length perc Element...... Aree 12.4858 
Number of Iterations....... : Z 
C) Structure Loading: 
a eo : 9000.0000 
Cerra eae are ~ -5000.0000 
0566 8 eer 1000.0000 2 
Ai. . eeeris< seeds 5 oe ee 0.0000 
D) Elemental Dimensions and Stress Distribution: 
Node Height Base Length Area 
iL 2.30012 1.89649 12.48578 4.49111 
Pa 0.84905 0.08490 12.48578 0.07209 
3 3.71044 0.71004 12.48578 2.63458 
4 3705 120 0.50168 12.48578 53373 
5 1750758 0.15986 12.40978 0225059 
6 ZeoO8 0.60867 ae 78 1.28654 
7 2acoo 70 0.73028 2.46978 1267139 
8 2.09919 0.65663 12.48578 17838 
9 0.85662 0.28601 B2.48578 0.24500 
E) Objective Function: 
Total structure Volume: 121.283012 
Node Normal Stress Bending Stress Total 
1 2513 487221 48973.4 
Z 18883.4 $3152.8 §2036..2 
3 422.4 S157 51999.5 
& 637 21 S16 8 51993). 1 
5 9637.3 42361.6 51998 .9 
6 2478.0 49521.0 51998 .9 
i] 2079.4 49920.1 5193975 
8 2346.2 49653.4 51999 .6 
9 1216521 a), 72 12165005 
F) Boundary Conditions: 
Node X-Displ Y-Displ Slope 
il 1 1 1 
9 1 1 0 
G) Solution Vector: 
Node X-Displ Y-Displ Slope 
1 Q.000000E+00 0.000000E+00 0.000000E+00 
2 0.427757E+00 -0.853866E-01 -0.525055E-01 
3 0.819383E£+00 -0.346183E+00 -0.786212E-02 
4 0.822885E+00 -0.351655E+00 0.527751E-02 
S 0.767984E+00 -0.809962E-01 0.313859E-01 
6 0.828458E+00 0.231809E+00 0.108871E-01 
7 O.841128E+00 0.252532E+00 -0.809329E-02 
8 0.653208E+00 0.1281152+00 -0.271114E-01 
9 0.000000E+00 0.000000E+00 -0.664133E-01 
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OPTIMIZATION #15a 


OPTIMIZATION SOLUTION 


A) Problem Parameters: 


Arch Angle : 180.000 Youngs Modulus: 30000000.0 
Arch Radius: 32.000 Yield Strength: §2000.0 
No of Design Var: Lis No of Elements: 8 
B) Derived Constants: 
No of System Nodal Points... 9 
No of Degrees of Freedon.... an 
Length per Slement.......... 12.4858 
Number of Iterations........ 2 
C) Structure Loading: 
ip PEP a Rare Sy hen 8 SR 0.0000 
eb PP ere rs oe no one 32000.0000 
BM. «. ocekeveuselacerete ares ee oo 1000.0000 = 
PA i ..< 4s 06) 4: bieueeiaeasnis ee e 100.0000 - 
D) Elemental Dimensions and Stress Distribution: 
Node Height Base Length Area 
1 27725280 2.04372 12.48578 6.64780 
zZ 3.20982 1.57865 12.48578 5.06717 
3 P0558 5 On 59193 12.48578 0.56164 
4 S4oli2 1.93749 12.48578 6.68769 
5 4.53938 Zeal 7 12.48578 10.90026 
6 3.87991 2.04928 t2.48578 7.95104 
7 Ph fie ew e206! 12.48578 5225150 
8 1.78269 Las O38 12.48578 Ze S655. 
9 1.43268 1.44887 12.48578 2.07577 
E) Objective Function: 
Total structure Volume: 516.579224 
Node Normal Stress Bending Stress Total 
i 323547 48783.9 S20Le.> 
2 3969.0 48027.9 §1996.9 
3 28718.2 2630.5 31348.8 
4 1520 56 50476.2 §1997.0 
5 574.6 Sasa 51998.7 
6 1027.3 50976.6 51997.8 
7 2359.0 49638.0 §1997.0 
8 6361.0 45636.9 §1998.0 
9 7761.0 One. welt 
F) Boundary Conditions: 
Node X-Displ Y-Displ Slope 
1 uu 1 1 
9 0 i 0 
G) Solution Vector: 
Node X-Displ Y-Displ Slope 
1 0.000G00E+00 0.000000E+00 0.000000E+00 
2 -0.800741E-01 0.174874E-01 0.125050E-01 
3 -0.427959E+00 0.254098E+00 0.443218E-01 
4 -0.654608E+00 0.596950E+00 0.100603E-01 
5 -0.666896E+00 0.660397E+00 -0.108817E-02 
6 -0.682870E+00 0.578787E+00 -0.114177E-01 
7 -0.811784E+00 0.384649E+00 -0.243721E-01 
8 -0.118582E+01 0.132794E+00 -0.440123E-01 
9 -0.183804E+01 0.000000E+00 -0.579555E-01 
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